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Dr.  Ronald  W.  Armstrong,  a  member  of  the  scientific  and  professional  cadre  of  senior 
executives,  is  Senior  Scientist  in  Energetic  Materials,  Munitions  Directorate,  Air  Force 
Research  Laboratory,  Eglin  Air  Force  Base,  Fla.  He  is  responsible  for  contributing 
technical  leadership  to  the  Air  Force  basic  research  mission  of  sustaining  superiority  in 
the  science  and  engineering  of  energetic  and  related  materials  properties  and 
performances.  In  this  activity,  he  provides  technical  reviews,  assesses  programs, 
establishes  research  strategies  for  new  study  areas,  and  performs  theoretical  and 
experimental  research.  He  interacts  with  other  Air  Force  and  government  agencies 
relating  to  energetic  material  technologies,  and  promotes  collaboration  with  foreign 
governments  and  organizations,  other  services  and  agencies,  universities  and  industry, 
as  relevant  to  the  directorate's  thrusts  and  long-range  goals  involving  energetic 
materiais  and  their  uses.  Research  topics  include  fundamental  aspects  of  the 
mechanics,  physics  and  chemistry  of  high  explosives:  the  mechanical  and  initiation 
behavior  of  nanoparficle  systems;  and  the  high-deformation-rate  and  shock  behavior  of 
condensed  matter. 

Dr.  Armstrong  has  spent  sabbatical  leave,  summer  and  consulting  appointments  at  U.S. 

Steel  Research  Laboratory.  Monroeville,  Pa.;  Oak  Ridge  National  Laboratory,  Tenn.; 

Lawrence  Livermore  National  Laboratory,  Calif.;  Institute  for  Defense  Analyses. 

Washington.  D.C.;  and  the  United  Technologies  Research  Center,  East  Hartford,  Conn, 

Additional  appointments  include  the  U.S.  National  Science  Foundation,  Science 
Education  Directorate,  Bethesda,  Md.;  Naval  Surface  Weapons  Center,  Silver  Spring,  Md.;  National  Institute  of  Standards  and 
Technology,  Gaithersburg,  Md.;  Naval  Research  Laboratory,  Washington,  D  C.;  U.S.  Office  of  Naval  Research,  London,  England; 
University  of  Cambridge,  England;  and  the  University  of  Strathclyde,  Glasgow.  Scotland.  Other  appointments  include  the  Naval 
Ship  Research  and  Development  Center,  Annapolis,  Md.;  Laboratory  for  Physical  Sciences,  College  Park,  Md.;  Johns  Hopkins 
University.  Baltimore,  Md.;  and  the  Naval  Surface  Warfare  Center,  Indian  Head,  Md. 

During  his  career.  Dr.  Armstrong  published  approximately  275  research  periodical  and  book  articles,  including  48  invited  articles. 
He  has  edited  several  books,  has  been  awarded  two  patents,  and  has  organized  a  number  of  national  and  international 
conferences  In  association  with  colleagues  and  students. 
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in  2003  Nano  Materials  for  Aerospace  Applications  Symposium:  Expanding  the 
Enyelope,  (Corpus  Christ!,  TX,  2003),  CD-ROM  (Some  charts  contain  data  from 
publications  in  the  open  literature  prior  to  MNME  appointment  and  others  are  cleared  yia 
AAC/PA  01-504,  -02-404,  -02-287,  -01-502,  -02-291,  -02-301,  and  -01-098). 

24.  R.W.  Armstrong,  B.  Baschung,  D.W.  Booth,  and  M.  Samirant,  “Enhanced 
propellant  combustion  with  nanoparticles”.  Nano  Letters,  3,  253  (2003)  (AAC/PA  02- 
401). 

III.  Articles  in  process  of  publication 

25.  *R.W.  Armstrong,  W.  Arnold,  and  F.J.  Zerilli,  “Grain  size  dependence  of  shock- 
induced  twinning  stresses  for  Armco  iron”.  Chemical  Physics  (Russian),  2003  (AAC/PA 
02-291). 

26.  *R.W.  Armstrong,  “Dislocation  mechanics  aspects  of  impact/shock  initiation  of 
detonation”.  Chemical  Physics  (Russian),  2003,  in  print  (AAC/PA  02-291). 

27.  *R.W.  Armstrong  and  W.L.  Elban,  “Dislocations  in  energetic  crystals”,  in 
Dislocations  in  Solids,  edited  by  E.R.N.  Nabarro  and  J.P.  Hirth,  (Elseyier  Science 
Limited,  London,  2003),  12,  in  print  (AAC/PA  01-502). 

28.  R.B.  Clough,  S.C.  Webb,  and  R.W.  Armstrong,  “Dynamic  hardness  measurements 
using  a  dropped  ball:  with  application  to  1018  steel”.  Materials  Science  and  Engineering, 
2003,  in  print. 

29.  R.W.  Armstrong,  “Crystal  imperfections  seen  by  x-ray  diffraction  topography”,  in 
Nano  and  Microstructural  Design  of  Adyanced  Materials,  edited  by  M.A.  Meyers,  R.O. 
Ritchie,  and  M.  Sarikaya,  (Elseyier  Science  Publications,  N.Y.,  2003)  (AAC/PA  02-551). 
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30.  H.W.  Haslach  ,  Jr.,  and  R.W.  Armstrong,  “Deformable  Bodies  and  Their 
Deformation  Behaviors”,  (John  Wiley  and  Sons,  Publishers,  N.Y.,  2003),  in  print. 


IV,  Articles/Presentations  submitted/to  be  submitted 

31.  M.Y.D.  Lanzerotti,  J.  Sharma,  and  R.W.  Armstrong,  “AFM  studies  of  fraeture 
surfaees  of  eomposition  B  energetic  materials”,  article  submitted  to  Metallurgical  and 
Materials  Transactions,  2003. 

32.  J.G.  Glenn  and  R.W.  Armstrong,  “The  seeds  of  tornado  prevention”,  article 
submitted  to  Atmospheric  Science  Letters,  2003  (AAC/PA  03-187). 

33.  P.  Rodriguez,  R.W.  Armstrong,  and  S.L.  Mannan,  “The  dependence  of  activation 
area  on  grain  size  in  cadmium”,  article  submitted  to  Transactions  of  the  Indian  Institute 
of  Metals,  2003. 

34.  K.L.  Kline,  R.W.  Armstrong,  D.W.  Richards,  and  M.P.  Kramer,  “Comparative 
model  x-ray  diffraction  characteristics  of  RDX  and  aluminum”,  article  submitted  to 
Shock  Compression  of  Condensed  Matter  -  2003,  American  Institute  of  Physics  Meeting, 
Portland,  OR,  July  20-25,  2003  (AAC/PA  03-272). 

35.  M.P.  Kramer,  A.E.D.M.  van  der  Heijden,  D.W.  Richards,  W.H.  Wilson,  R.W. 
Armstrong,  and  A.C.  van  der  Steen,  “Coating  and  recrystallization  of  energetic/reactive 
nanomaterials”,  abstract  submitted  for  the  Sixth  International  Workshop  on  Crystal 
Growth  of  Organic  Materials  (6-ICGQM),  Glasgow,  U.K.,  August  17-21,  2003  (AAC/PA 
03-361). 

36.  *R.W.  Armstrong,  M.P.  Kramer,  and  W.H.  Wilson,  “The  Power  and  Strength  of 
Energetic/Reactive  Nanomaterials”,  abstract  submitted  for  the  International  Workshop  on 
Novel  Energetic  Materials  and  Applications  (9-IWCP),  Eerica,  Ea  Spezia,  IT,  September 
14-18,  2003  (AAC/PA  03-269). 

37.  W.R.  Grise  and  R.W.  Armstrong,  “I-D  to  3-D  dislocation  pile-ups  adiabatically 
released”,  abstract  submitted  for  The  Minerals,  Metals,  and  Materials  Society  (TMS)  3- 
Dimensional  Materials  Science  Symposium,  Chicago,  IE,  November  9-12,  2003. 

38.  *R.W.  Armstrong,  H.  Conrad,  and  E.R.N.  Nabarro,  “Meso-to-nano-scopic 
polycrystal/composite  strengthening”,  abstract  submitted  for  the  Materials  Research 
Society  (MRS)  Symposium  Q:  Mechanical  Properties  of  Nanostructured  Materials  and 
Nanocomposites,  Boston,  MA,  December  1-5,  2003. 

39.  M.Y.D.  Lanzerotti,  J.  Sharma,  and  R.W.  Armstrong,  ’’Nanofractography  of 
composition  B  fracture  surfaces  with  AEM”,  abstract  submitted  for  the  Materials 
Research  Society  (MRS)  Symposium  AA:  Synthesis,  Characterization,  and  Properties  of 
Energetic/Reactive  Nanomaterials,  Boston,  MA,  December  1-5,  2003. 
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40.  W.R.  Grise  and  R.W.  Armstrong,  “Hot  spots  from  crystal  size  dependent  disloeation 
pile-up  avalanehes  in  energetie  materials”,  abstraet  submitted  for  the  Materials  Researeh 
Soeiety  (MRS)  Symposium  AA:  Synthesis,  Charaeterization,  and  Properties  of 
Energetie/Reactiye  Nanomaterials,  Boston,  MA,  Deeember  1-5,  2003. 

41.  L.  Ferranti,  Jr.,  R.W.  Armstrong,  andNareshN.  Thadhani,  “Elastie/plastie 
deformation  behayior  in  a  eontinuous  ball  indentation  test”,  draft  artiele  prepared  for 
submittal  to  Materials  Seienee  and  Engineering  A,  August,  2003. 

B.  Lectures/Briefings 

1 .  “Energetie  erystal  properties  and  model  analyses”;  Inyitational  Workshop  on  Nano-, 
energeties,  Uniyersity  of  Florida  Graduate  Engineering  and  Researeh  Center,  Shalimar, 
EE;  Noyember  16,  2000. 

2.  *“Disloeation  meehanics/mierostruetural  deseriptions  of  strength  properties  of  metals 
and  alloys”;  THERMEC’2000,  International  Conferenee  on  Proeessing  and 
Manufaeturing  of  Adyaneed  Materials,  Eas  Vegas,  NV;  Deeember  5,  2000  (see  A,II,14). 

3.  “Nanoseienees  teehnology”;  AFOSR  Program  Offieers’  Visit,  Eglin  Air  Foree  Base, 
EE;  February  8,  2001. 

4.  *  “Charaeterization  of  eleotronie  deyiee  and  energetie  erystal  properties”  (presented 
by  M.  Dudley);  International  Workshop  on  Properties  and  Charaeterizations  of 
Teehnologieally  Important  Single  Crystals  (PCSC  -  2001);  New  Delhi,  India;  February 
27,  2001. 

5.  “Grain  size  dependent  alumina  fracture  meehanies  stress  intensity”;  International 
Conferenee  on  the  Seienee  of  Hard  Materials  (ICSHM7);  Ixtapa,  Mexieo;  Mareh  5,  2001 
(see  A.L2)  (AAC/PA  01-098). 

6.  *“Air  Foree  research  in  nano-energetics  -  potential  adyantages  and  ehallenges” 
(presented  by  W.H.  Wilson);  22  Ameriean  Chemieal  Soeiety  National  Meeting,  San 
Diego,  CA;  April  3,  2001  (see  A.IL15)  (AAC/PA  01-130). 

7.  “Nanoenergeties”  (with  W.H.  Wilson  and  M.P.  Kramer);  AFOSR  yisit  of  Eyle 
Sehwartz,  AFRL/MNME,  Eglin  Air  Foree  Base,  EE;  April  10,  2001. 

8.  “Nanoseienee  teehnology”,  AFOSR  Presentation,  Arlington,  VA;  May  29,  2001. 

9.  *“AFRL  strategic  workshop  on  nanoscience  and  teehnology”,  Uniyersal  Technology 
Corporation,  Dayton,  OH;  May  30,  2001  (Proprietary). 

10.  *“Energetie  erystal  researeh  actiyity”;  Institut  Saint-Eouis  (ISE),  France;  June  18, 

2001. 
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1 1 .  *“Constitutive  equation  research  activity”;  Fraunhofer  Institut  fur  Kurzzeitdynamik, 
Ernst-Mach-Institut  (EMI),  Ereiburg,  Germany;  June  19,  2001. 

12.  (Poster)  “Nanoscience  mechanics”;  AEOSR  Joint  Metallic  and  Ceramic  Materials 
Programs  Review,  Snowbird,  UT;  August  21,  2001  (see  A.ILIV)  (AAC/PA  01-367). 

1 3 .  *“lmpact  of  nanoenergetics”;  National  Academy  of  Science/National  Research 
Council  (NAS/NRC)  Committee  on  Implications  of  Energetic  Micro-  and  Nano- 
Technologies;  Albuquerque,  NM;  November  8,  2001;  listed  in  “Implications  of  Emerging 
Micro-  and  Nanotechnologies”,  National  Research  council  of  The  National  Academies, 
Washington,  D.C.,  2002,  p.241. 

14.  “6.1  nanoenergetics”;  Scientific  Advisory  Board  (SAB)  Review,  High  Explosives 
Research  and  Development  (HERD),  Eglin  Air  Eorce  Base,  EE;  November  29,  2001 
(Proprietary). 

15.  “Advanced  energetics  briefing”;  General  Chedister,  Air  Armament  Command 
(AAC),  High  Explosives  Research  and  Development  (HERD),  Eglin  Air  Eorce  Base,  EE; 
January  3,  2002  (Proprietary). 

16.  *”Nanometric  aspects  of  energetic  crystal  properties/testing”;  Mardi  Gras: 
Biotechnology  Conference,  Eouisiana  State  University,  Baton  Rouge,  EA;  Eebruary  8, 
2002  (see  A.IL18). 

17.  *“Armco  iron  twinning/grain  size  dependence”;  Shock  Wave  Processes  in 
Condensed  Matter,  Edinburgh,  Scotland;  May  20,  2002  (see  A,III,25)  (AAC/PA  02-291). 

18.  *”Model  studies  of  explosive  crystal  impacts/shocks”;  Shock  Wave  Processes  in 
Condensed  Matter,  Edinburgh,  Scotland;  May  21,  2002  (see  A.III,26)  (AAC/PA  02-291). 

19.  *“Dislocation  mechanics  activities”;  Physics  and  Chemistry  of  Solids  Colloquium 
Series,  Cavendish  Eaboratory,  University  of  Cambridge.  U.K.,  May  30,  2002  (AAC/PA 
02-178  and  -02-287). 

20.  *”Discussant  overview;  energetic  nanomaterials”;  Gordon  Research  Conference: 
Energetic  Materials,  Tilton,  NH;  June  16,  2002. 

21 .  *“Dislocation  mechanics  for  dynamic  plasticity/fracturing”;  14^*^  National 
Theoretical  and  Applied  Mechanics  Congress,  Virginia  Polytechnic  Institute  and  State 
University,  Blacksburg,  VA;  June  24,  2002  (see  A,II,19). 

22.  (Poster)  “The  power  of  nano-energetics”  (presented  by  M.P.  Kramer);  29^*^ 
International  Pyrotechnics  Seminar,  Westminster,  CO;  July  16,  2002  (see  A,II,20). 
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23.  (Poster)  “Nanoscience  mechanics  of  materials”  (presented  by  S.L.  Semiatin); 
AFOSR  Joint  Materials  Programs  Reviews,  Bar  Harbor,  ME;  August  10,  2002  (see 

A.IL21). 

24.  “Crystal  imperfections  seen  by  x-ray  topography”;  The  Minerals,  Metals,  and 
Materials  Society  (TMS)  Fall  Meeting,  Columbus,  OH;  October  9,  2002  (see  A,III.29) 
(AAC/PA  02-551). 

25.  *(Poster)  “Dislocation  mechanics  based  constitutive  equations”,  US  -  Korea 
Workshop  on  Advances  in  Metallic  Structural  Materials,  AFOSR  Advanced  Computing 
Center,  Kehei,  Maui,  HI;  January  21,  2003  (see  A.II.22)  (AAC/PA  01-361,  -02-287,  and 
-02-291). 

26.  *“  The  power  of  energetic/reactive  nano-materials”;  2003  Nano  Materials  for 
Aerospace  Applications;  Expanding  the  Envelope,  Corpus  Christi,  TX;  January  29,  2003 
(see  A.IL23)  (AAC/PA  01-504,  -02-404,  -02-287,  -01-502,  -02-291,  -02-301,  and  -01- 
098). 

27.  *”The  wonders  of  high  strain  rate  dislocation  mechanics?”;  Chemistry  Division,  Eos 
Alamos  National  Eaboratory,  Los  Alamos,  NM;  February  19,  2003. 

*  invited  participations 


C.  Degree  Associations 


1.  External  examiner;  Ph.D.  thesis  of  Dr.  John  Corley,  “Characterization  and  modeling 
of  the  dynamic  mechanical  properties  of  a  particulate  composite  material”,  Universitat 
der  Bundeswehr  Munchen,  Germany,  June  28,  2002. 

2.  Advisory  Committtee:  M.Sc.  thesis  of  Kinnan  L.  Kline,  “Comparative  Model 
Powder  X-ray  Diffraction  Characteristics  of  Aluminum  and  Cyclotetramethylene- 
tetranitramine  (RDX)”,  University  of  Florida,  Graduate  Engineering  and  Research 
Center,  November  26,  2002. 


D.  Patent  Disclosure 


J.G.  Glenn  and  R.W.  Armstrong,  “Tornado  disruption  through  nanometric  sized  cloud 
seeding”,  OMB-No.  900  -  0095,  May,  2002. 

E.  Associated  Symposia 

1 .  M.P.  Kramer,  W.H.  Wilson  and  R.W.  Armstrong,  “Invitational  Workshop  on  Nano¬ 
energetics”,  University  of  Florida,  Graduate  Engineering  and  Research  Center,  Shalimar, 
EL,  November  15-16,  2000. 
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2.  G.T.  Gray,  M.A.  Meyers,  N.N.  Thadhani,  and  K.S.  Veeehio,  “Symposium  in  honor 
of  Professor  Ronald  W.  Armstrong,  Dynamic  deformation;  constitutive  modeling,  grain 
size,  and  other  effects”.  The  Minerals,  Metals,  and  Materials  Society  (TMS)  Annual 
Meeting,  San  Diego,  CA,  March  3-5,  2003. 

3.  N.N.  Thadhani,  R.W.  Armstrong,  W.H.  Wilson,  J.J.  Gilman,  Z.  Munir,  and  R.L. 
Simpson,  co-organizers  of  Materials  Research  Society  (MRS)  Symposium  AA; 

Synthesis,  Characterization,  and  Properties  of  Energetic/Reactive  Nanomaterials,  Boston, 
MA,  December  1-5,  2003. 
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Abstract 

Various  grain  size  dependent  hardness  test  results  obtained  on  alumina  materials  are  brought  together  on 
a  Hall-Peteh  (H-P)  inverse  square  root  of  grain  diameter  basis  to  illustrate  the  normal  expeetation  of  greater 
strengthening  for  finer  grained  material,  not  unrelated  to  the  observation  that  struetural  eeramie  materials 
are  employed  often  in  a  very  mueh  finer  grain  size  eondition  than  eounterpart  struetural  metals  and  alloys. 
Here,  the  H-P  dependenee  is  shown  to  earry  over  to  deseribing  reeent,  mutually  eonfirming,  hardness-based 
measurements  of  the  fraeture  meehanies  stress  intensity  of  alumina  materials  at  small  grain  sizes.  The 
somewhat  surprising  dependenee  relates  to  the  frequent  observation  that  an  inerease  in  plastie  flow  strength 
of  a  struetural  material  is  assoeiated  with  a  deerease  in  its  fraeture  toughness.  For  the  ease  of 
strengthening  by  grain  size  refinement,  however,  the  fraeture  stress  of  a  material  is  raised  more  than  the 
plastie  flow  stress  and  so  greater  loeal  plastie  work  is  required  for  fraeturing. 

1.  Hardness/grain  size/cracking  dependencies 

Considerable  results  are  available  to  show  that  the  hardnesses  of  both  cubie  and  non-cubie  eeramie 
materials  are  greater  for  materials  of  finer  polyerystalline  grain  sizes,  for  example,  as  established  from 
various  measurements  made  on  eubie  magnesia  and  rhombohedral  alumina  polyerystals.  The  dependenee 
is  often  demonstrated  on  an  inverse  square  root  of  grain  diameter  basis  as  shown  for  polyerystalline 
alumina  materials  in  Figure  1  that  has  been  adapted  from  measurements  reported  by  Riee,  Wu  and 
Borehelt  [1].  The  adapted  Figure  eontains  from  the  original  work  only  the  Viekers  hardness  measurements 
reported  at  two  applied  load  values  for  mainly  hot  pressed  alumina  material,  as  open  and  elosed  triangle 
points,  but  to  whieh  have  been  added  newer  measurements  reported  by  Krell  [2]  for  pressureless  sintered 
material  (open  eirele  points)  and  by  Franeo,  Roberts  and  Warren  [3]  for  three  grain  sizes  (elosed  eirele 
points).  Franeo  et  al.  employed  pressureless  sintering  for  their  two  largest  grain  size  materials  and  hot 
pressing  during  sintering  for  their  smallest  grain  size  material. 
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The  lines  in  Figure  1  are  those  drawn  by  Rice  et  al.  [1]  for  their  measurements  at  the  indicated  applied 
loads  of  100  and  500  g.  Rice  et  al.  associated  the  lowest  hardness  values  obtained  at  500  g,  shown  at  the 
break  in  the  hardness  lines,  with  cracking  and  observed  crack  sizes  being  comparable  to  the  material  grain 
size.  The  upper  closed  circle  points  of  Franco  et  al.  were  obtained  at  a  smallest  applied  load  of  500  g  and 
the  lower  closed  circle  points  were  obtained  as  essentially  constant  hardness  values  over  a  range  of  loads 
from  1.0  to  25  kg.  Cracking  was  reported  for  all  of  the  measurements  reported  by  Franco  et  al.,  with 
comparable  crack  and  hardness  indentation  diagonal  lengths  occurring  at  the  smallest  500g  load.  The 
results  shown  from  Krell,,  for  ground  surfaces  on  pure  alumina  microstructures,  were  obtained  at  a  load  of 
10  kg.  Krell  obtained  slightly  lower  hardnesses  at  the  same  applied  load  for  polished  materials. 


Grain  Size  fim 
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Fig.  1.  The  Vickers  microhardness  dependence  on  grain  size,  adapted  from  results  reported  by  Rice,  Wu 
and  Borchelt  [1],  and  incorporating  additional  measurements  from  Krell  [2]  and  Franco,  Roberts  and 
Warren  [3]. 

The  new  results  added  to  Figure  1  appear  to  be  in  agreement  with  the  emphasis  given  by  Rice  et  al.  to  the 
importance  of  cracking  in  affecting  the  determined  hardness  values.  A  similar  effect  of  cracking  on 
reducing  the  hardness  of  MgO  single  crystals  was  demonstrated  for  reported  single  crystal  measurements 
[4].  Flowever,  such  cracking  is  only  one  consideration  relating  to  the  increase  in  hardness  that  occurs  with 
reduced  applied  load.  Detailed  measurements  made  on  residual  indentations  in  the  {100}  cleavage  surface 
of  an  MgO  crystal  [5]  and  computation  of  the  expected  Flertzian  elastic  contribution  to  the  load  bearing 
capability  [6]  have  provided  evidence  for  an  important  elastic  influence  on  this  so-called  indentation  size 
effect  (ISE).  Cracking  reduces  the  contribution  of  elastic  recovery  to  the  ISE.  Otherwise,  the  occurrence 
of  controlled  indentation  cracking  experiments  has  led  to  the  method  of  indentation  fracture  mechanics, 
either  by  employing  Vickers  hardness  measurements  or  Flertzian  ball-type  ring  cracking  methods  for 
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determining  the  fracture  toughness  of  alumina  polycrystals  [3,7].  In  the  present  report,  recent  indentation 
fracture  mechanics  stress  intensities  determined  for  both  testing  methods  are  shown  to  follow  predictions  of 
an  H-P  inverse  square  root  of  grain  diameter  dependence  [8,9]. 

2.  A  Hall-Petch  dependence  for  the  fractnre  mechanics  stress  intensity. 

2. 1 .  Collected  experimental  results 

Franco  et  al.  [3]  measured  the  fracture  mechanics  stress  intensities  for  their  materials  by  the  two  methods 
of  Hertzian  ring  cracking  and  Vickers  indentation  tests.  For  the  Hertzian  case,  well-defined  ring  cracks 
were  obtained  in  the  minimum  load  range  of  30  to  50  kg  applied  at  a  rate  of  0. 14  mm/min  to  a  2.5  mm 
radius  alumina  sphere.  The  investigators  gave  emphasis  to  smaller  individual  crack  sizes  being  obtained  in 
the  Hertzian  tests,  say,  at  about  5  micrometer  crack  lengths  or  smaller,  as  compared  to  larger  crack  sizes 
occurring  in  the  Vickers  test  results  obtained  at  applied  loads  of  0.5  to  25  kg.  The  Vickers  measurements 
showed  a  decrease  in  the  grain  size  dependence  of  the  fracture  mechanics  stress  intensity  at  greater  loads 
until  reaching  a  negligible  dependence  for  crack  sizes  of  about  40  micrometers  and  thereafter  showing  a 
negative  H-P  dependence. 
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Fig.  2.  A  Hall-Petch  description  of  the  Hertzian-determined  indentation  fracture  mechanics  stress  intensity 
dependence  on  grain  size  for  fine  grained  alumina  as  measured  by  Franco,  Roberts  and  Warren  [3]  and,  for 
low-load  Vickers  measurements,  comparable  measurements  reported  by  Muchtar  and  Lim  [7],  also 
including  in  the  Figure  a  range  of  single  crystal  values  reported  by  Rice  [10]  along  with  indication  of 
thermal  strain  influence  complications  at  large  grain  sizes. 
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The  Hertzian  fraeture  meehanies  stress  intensity  results  obtained  by  Franeo  et  al.  shown  in  Figure  2  are 
eompared  with  other  Viekers  stress  intensity  determinations  reported  by  Muehtar  and  Lim  [7]  for  indenter 
load  values  of  0.3,  0.5,  1.0  and  2.0  kg.  The  latter  measurements  give  essentially  parallel  H-P  dependeneies 
that  are  raised  at  inereasing  applied  loads.  Also  shown  in  Figure  2  along  the  vertieal  ordinate  axis  is  a 
braeketed  range  of  single  erystal  hardness  values  reported  in  a  reeent  eompilation  by  Riee  [10].  In  the  same 
eompilation,  widely  varying  stress  intensity  values  were  indieated  at  larger  values  of  polyerystalline  grain 
diameters,  say,  greater  than  ten  mierometers,  as  is  meant  to  be  indieated  in  the  Figure  by  the  vertieal 
arrows,  and  appears  attributable  to  thermal  strain  influenees.  The  total  range  in  stress  intensity  values 
shown  in  Figure  2  may  be  eompared  with  a  value  of  3.7  MPa.m'^'quoted  by  Congleton  [11]  for  dynamie 
fraeturing  studies  eondueted  on  pre-eraeked  thin  plates. 

2.2  Model  Hall-Peteh  interpretation 

Returning  to  Figure  1,  it  is  important  that  an  inerease  in  the  alumina  yield  stress  with  reduetion  in  grain 
size  is  indieated  for  the  H-P  type  hardness  dependenee.  Generally,  an  inerease  in  the  yield  stress  of  a 
material  is  assoeiated  with  a  reduetion  in  the  fraeture  meehanies  stress  intensity.  The  reason  is  that  sueh 
ehanges  often  do  not  affeet  the  material  fraeture  stress  and  so,  during  loading  of  a  pre-eraeked  tensile 
speeimen,  less  strain  hardening  is  neeessary  at  the  eraek  tip  in  order  to  reaeh  the  material  fraeture  stress  [9]. 
For  strengthening  produeed  by  grain  size  refinement,  however,  the  yield  stress  is  raised  but  the  fraeture 
stress  is  raised  more  on  an  H-P  basis  than  the  yield  stress  and,  so,  the  eounter  result  of  greater  strain 
hardening  must  oeeur  at  the  eraek  tip  before  the  fraeture  stress  is  reaehed.  The  basis  for  a  greater  H-P 
dependenee  for  the  fraeture  stress,  eompared  to  the  yield  stress,  rests  on  the  theoretieal  assoeiation  of 
disloeation  eoaleseenee  with  formation  of  a  eraek  nueleus. 

For  alumina,  as  in  the  brittle  eleavage  of  steel,  mieroyielding  is  presumed  to  oeeur  at  the  eraek  tip  in 
order  to  promote  further  eraek  growth  and/or  to  nueleate  new  eraeks  at  about  the  same  size  seale  as  the 
grains  themselves.  Evidenee  for  a  substantial  plastie  work  eontribution  to  the  experimentally  measured 
fraeture  meehanies  stress  intensity  eomes  from  evaluation  of  the  theoretieal  surfaee  energy  term. 

Congleton  pointed  out  that  the  stress  intensity  determined  value  of  approximately  18  J/m^  was  far  larger 
than  indieated  for  the  theoretieal  value.  Reeently,  a  value  of  2.15  J/m^has  been  reported  for  alumina  [12], 
in  agreement  with  the  previous  eonjeeture. 

For  eraek  sizes  very  mueh  larger  than  the  material  grain  size,  the  plastie  zone  size  at  the  tip  of  the  larger 
eraek  should  eonneet  with  the  material  grain  size  and,  thus,  a  reversal  from  the  H-P  dependenee  may  oeeur 
[8].  Thus,  the  reversed  H-P  dependenee  shown  by  applieation  of  the  higher  loads  in  the  Viekers  tests  of 
Franeo  et  al.  may  be  explained  on  this  basis. 

3.  (Grain/plastic  zone)  size  influences  on  the  flexural  strength. 

To  eonneet  with  the  grain  size  dependent  stress  intensity  measurements,  several  features  of  the  flexural 
(bend)  strength  dependenee  on  grain  size  for  variously  maehined  alumina  materials  are  shown  on  an  H-P 
basis  in  Figure  3,  also  adapted  from  a  eompilation  of  measurements  reported  by  Riee  [13];  see  his  Figure  3 
that  eompares  a  great  amount  of  data  and  builds  onto  earlier  work.  First,  in  the  present  Figure  3,  the  steep 
line  through  the  origin  with  slope  of  3.5  MPa  m'^  was  shown  to  deseribe  the  fraeture  stress  dependenee  on 
grain  size  at  relatively  large  grain  sizes.  The  pair  of  lines  with  relatively  lower  slopes  and  turned  down 
segments  at  their  large  grain  size  ends  were  ehosen  to  deseribe  two  flexural  strength  features:  (1)  in  the 
ease  of  the  lower  line,  a  lower  limiting  boundary  that  was  established  for  smaller  grain  size  fraeture 
strength  measurements  from  a  previous  survey;  and,  (2)  for  the  higher  segmented  line,  a  ehosen 
dependenee  for  numerous  but  seattered  fraeture  strength  measurements  when  the  plotted  grain  diameters 
were  determined  at  the  fraeture  origins.  It  should  be  noted  that  the  eolleeted  measurements  ineluded 
fraeture  strength  values  following  the  line  with  steepest  slope  to  a  flexural  strength  value  of  approximately 
600  MPa  at  a  grain  size  of  approximately  35  mierometers. 

It  seems  reasonable  to  take  the  steep  line  in  Figure  3  as  the  fraeture  stress  dependenee  on  Riee’s 
suggestion  of  the  material  eraek  size  being  determined  in  depth  by  one-half  of  the  material  grain  size,  and, 
on  this  basis,  the  slope  value  of  3.5  MPa  m'^  is  in  reasonable  agreement  with  the  stress  intensity  values 
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shown  in  Figure  2.  The  dashed  line  of  lesser  slope  through  the  origin  in  Figure  3  has  been  newly  eomputed 
for  a  fracture  stress  dependence  on  crack  size,  also  taken  in  depth  equal  to  one-half  of  the  material  grain 
diameter,  and  with  a  theoretical  value  of  the  stress  intensity  employing  a  Young’s  modulus  of  420  GPa  and 
theoretical  fracture  surface  energy  of  2.15  J/m^;  see  the  Figure  legend.  As  mentioned  earlier,  Congleton 
had  pointed  out  that  the  measured  stress  intensity  even  for  rapid  fracturing  in  his  macroscopic  dynamic 
fracture  tests  gave  indication  of  a  significant  plastic  flow  contribution  to  the  alumina  material  fracture 
toughness,  as  apparently  confirmed  on  the  Fl-P  basis  of  Figure  3. 


t,  M-m 


10  1.0  0.3 


Fig.  3.  Linear  indications  from  collected  measurements  of  the  flexural  (bend)  strength  dependence  on  grain 
size,  as  adapted  from  Rice  [13],  in  comparison  with:  (1)  a  lower  limiting  (dashed)  theoretical  fracture  stress 
dependence  on  grain  size  when  determining  the  material  crack  size;  and,  (2)  a  model  description  fitted  to 
the  observed  fracture  stress  dependencies  and  incorporating  influence  of  a  plastic  zone  size  correction  at  the 
crack  tip. 

Also  shown  in  Figure  3,  via  closed  triangular  points,  are  pre-cracked  fracture  stress  values  computed  in 
terms  of  a  reference  crack-free  fracture  stress  and  the  material  crack  radius,  c,  with  an  associated  plastic 
zone  extension,  s,  at  the  crack  tip  [8,14];  see  the  legend  in  Figure  3.  For  this  model  description,  a  limiting 
small  value  of  s  applies  for  the  theoretical  surface  energy  and  s  becomes  increasingly  larger  to  match  the 
greater  stress  intensity  values  shown  in  Figure  2.  The  higher  computed  fracture  stress  (filled  triangular) 
points  in  Figure  3  were  arrived  at  by  a  two-step  procedure.  First,  a  relevant  value  of  the  fracture  stress 
intensity,  K,  was  selected  from  Figure  2,  along  with  an  estimated  crack-free  fracture  stress  taken  a  large 
grain  size  from  the  upper  line  in  Figure  3,  to  provide  an  estimated  plastic  zone  size,  s.  Then,  an 
appropriate  pre-cracked  fracture  stress  was  computed  from  the  selected  crack-free  fracture  stress  and  s. 
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Such  pre-cracked  plastic-  zone-size-determined  fracture  stresses  rapidly  converged  to  the  crack-free 
fracture  stress  values.  The  values  of  s  decreased  from  0.026  mm  at  a  grain  diameter  of  0.25  mm  to  a 
smallest  size  of  0.012  mm  at  the  smallest  grain  size  limit  of  just  less  than  0.0003  mm  in  the  Figure,  thus, 
indicative  for  the  smallest  grain  size  material  of  significant  microyielding.  For  the  lower  triangular  points 
in  Figure  3,  match  was  made,  first,  at  the  largest  grain  size,  with  the  dashed  theoretical  stress  intensity  line 
for  K  =  1.34  MPa  m''^  and  an  estimated  crack- free  fracture  stress  of  130  MPa  to  give  a  plastic  zone  size,  s 
=  0.19  mm  for  the  0.25mm  grain  size,  and  so  on. 

For  the  foregoing  model  description  then,  the  (crack/plastic  zone)  sizes  are  shown  to  tie  stress  intensity  and 
fracture  stresses  together,  so  as  to  give  at  large  crack  sizes,  taken  approximately  equal  to  the  material  grain 
size,  a  conventional  fracture  mechanics  type  dependence  on  the  material  grain  size;  and,  at  small  grain 
sizes,  a  plastic  zone  size  substantially  larger  than  the  material  grain  size  and  so  giving  an  apparent  Fl-P 
dependence  of  the  fracture  stress  measurements. 
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ABSTRACT 


The  experimental  observation  of  an  inverse  square 
root  of  grain  size  influence  on  the  strength 
properties  of  reasonably  single  phase 
polycrystalline  materials  has  been  explained  by 
model  calculations  of  grain  boundaries  blocking 
temporarily  the  egress  of  dislocation  pile-ups  in 
embryonic  slip  bands  or  deformation  twins  [I]. 

Metallurgical  composite  connection 

Extension  of  the  inverse  square  root  of  grain  size 
dependence  to  metallurgical  composite  materials 
has  been  arguably  successful  on  the  basis  of 
determining  an  effective  grain  size  for  the  overall 
composite  through  measurement  of  the  mean  free 
intercept  line  lengtli  involving  both  interphase 
boundaries  and  grain  boundaries  [2,3].  Significant 
applications  have  been  to  very  small  effective  grain 
sizes  for  patented  eutcctoid  steel  (piano)  wire  [4] 
and  unidirectionally-solidified  eutectic,  superlattice 
(candidate  turbine  blade),  materials  [5]. 

Nanoscale  strength  connection 

Extrapolation  of  the  effective  grain  size  analysis  to 
nanoscale  dimensions  has  produced  a  number  of 
interesting  consequences:  (i)  the  achievement  of 
experimental  strength  levels  near  to  tlie  theoretical 
limiting  strength  of  defect-free  material  [6]  as 
shown  in  Figure  I;  (ii)  theoretical  prediction  of 
discontinuous  strength  level  changes  at  limiting 
small  grain  sizes  [7];  and,  (iii)  suggestion  of  a 
reverse  trend  to  effective  grain  size  weakening 
caused  by  activation  of  grain  boundary  deformation 
mechanisms  at  the  smallest  grain  sizes  [8].  The 
last  mentioned  consideration  relates  to  the 
observation  of  effective  high  temperature  creep 
deformation  nieclianisms  [9j. 

Thermally-activated  flow 

Distinction  between  low-temperature,  dislocation- 
based,  model  descriptions  of  viscopiastic  behavior 
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Figure  1.  Inverse  square  root  of  grain  size 
dependence  for  iron  and  steel  materials  [6]. 

and  higher  temperature  diffusion-based  flow 
apparently  rests  on  direct  consideration  in  the 
former  case  of  an  appropriate  stress  influence  on  the 
required  Gibbs  free  energy  for  dislocation 
activation,  G  [10],  Such  consideration  enters 
generally  into  the  so-called  friction  resistance  force 
that  operates  on  each  dislocation  within  an 
embryonic  slip  band,  rather  than  being  an  influence 
on  the  microstructural  stress  intensity  coefficient  of 
the  reciprocal  square  root  of  grain  size  dependence. 

The  thermally-activated  dislocation  model 
description  provides  for  significant  differences  in 
crystal-structure-dependent  flow,  for  example, 
between  strain-independent  thermal  stress  terms  for 
iron  and  other  bcc  metal  deformation  as  compared 
with  a  strain-dependent  thermal  stress  for  coppwr 
and  related  fee  metals  [11].  Such  thermally- 
activated  description  has  been  shown  to  be 
applicable,  as  well,  to  certain  polymer  material 
deformation  properties,  as  shown  below  in  Figure  2 
for  reported  compression  results  on  PMMA  [12]. 


Shear  banding 
41 
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Beyond  consideration  of  the  stress-assisted  thermal 
energy  influence  on  individual  dislocation  or  flow 


G/k  (K) 

Figure  2.  Yield  stress/shear  modulus  ratio  versus 
G/k,  [Boltzmann’s  constant],  for  PMMA  [12]. 


unit  creation  and  movement,  there  is  the  afore¬ 
mentioned  important  added  influence  of  local  stress 
elevations  caused  by  the  collective  action  of  the 
multiple  flow  units  in  the  embryonic  shear  zones 
blocked  at  internal  obstacles.  The  latter  influence 
is  a  minor  consideration  for  aluminum  and  a  major 


Crystal  Size,  t 


Figure  3.  Crystal  size  dependence  of  the  drop- 
weight  impact  initiation  of  explosives  [13], 

one  for  steel,  and  is  even  more  important  for 
understanding  the  deformation  properties  of 
energetic  (explosive)  materials,  such  as  RDX 
(cyclotrimethylenetrinitramine),  for  which  the 


mechanical  initiation  of  chemical  decomposition  in 
drop-weight  impact  tests  is  shown  in  Figure  3  to  be 
dependent  on  crystal  size  [13]. 

Ductility  connection 

Although  both  reduction  in  crystal  or  grain  size  and 
addition  of  particles  add  to  material  strengthening, 
the  influence  on  ductility  is  potentially  increased  in 
the  former  case  and  clearly  reduced  in  the  latter 
case  [14],  thus  composite  material  strengthening 
with  particles,  even  at  nanocrystalline  grain  sizes,  is 
obtained  at  the  price  of  ductility  reduction.  One 
model  consideration  being  further  investigated  is 
based  on  the  proposal  that  the  true  fracture  strain  is 
detennined  by  tlie  particles  being  brought  by  plastic 
deformation  to  a  smaller  critical  spacing  [15]. 
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The  limits  of  conventional  CHNO  chemistry  are  in  sight 
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Nanometric  Fuels 


2 

O 

P  O 


P 

<0 

Q> 

c 

TO 


o  t: 

:  tX. 


to 

<i> 

%m0 

q: 

Q> 

to 

<C 

-2 

q: 

I 

ti) 

£ 

k. 

Q> 


Q> 

to 

TO 

TO 

TO 

QC 

g> 

TO 

£ 

P  o' 

§!  C 
P 


2  P 

P 

•S’  C  s2 
a:  £u4 

I  I 


c 

TO 

'*X 

C 

o 

o 

rk, 

g>g 
TO  o 
C  *2- 

£  TO 

TO  g 

■c  F 
•2>o 
p  £ 


o 

JO 

(0  ^ 

Q)  C 

O  m 
C 

TO  P 
C  O 

.S> 

Jfc  Q. 

Q  0) 

gw 

CO  o 

3  P 
00  TO 
.  0) 


^-  TO 
O  TO 
^  TO 

c  +- 

1° 

o 

2  to  ^ 
TO  ><5 
•O^x 
C  0.  = 

3  -o  ^ 

TO  TO  C 
c  J>*TO 

2  F  TO 
E  .to  TO 
TO  'E  'C 

^  C  0) 
c  TO 
SJZ  TO 

u.  U  ” 

I 


O) 

TO 

•a 

re 

TO 

L. 

0) 

Q 

>. 

(Q 

**- 

_l 

O 

C 

TO  .2 

S’  TO 

-  ■% 
C  TO 

TO  2 

e  ^ 

TO  -r 
0.  o 
X  = 
O) 


TO 

TO  X 

<5 

TO  m 
o  p 
(0  0- 

=  2. 

w  ^ 

o  TO 
it=  O 

O  (A 

TO  > 
Q. 

(O  •*= 
u> 

II 
•*■  0) 
> 


S  X 


TO 

> 

I 


re 

TO 


TO 

o 

re 

P 

3 

(O 


TO 

(A 

TO 

TO 

c 

TO  i’ 
S  TO 


<0 

Q> 

O 

c 

2 

s£ 


re 

o 

CO 


u 

TO 

Q. 

CO 

£ 

O) 

X 

TO 

> 


™  -C 

U>  4-> 

TO  re 
S  0. 

°  i 

(A  •— 
TO  ^ 

•o  ^ 

55 

4  ^ 

o 


(O 

I 


c 

§ 

o 

Q. 


c 

o 

2 

re 

o 

c 

re 


O  Uj 


p 

re 

re 

c 

re 


re 

o 

c 

re 

g 

Q 


o  t: 

.  re 

:  0. 


to 

TO 

to  o 
TO  •— 

0^  S. 

TO  — 
to  TO 
re  3 
TO 


>> 

O) 

TO  to 


TO  LL  c  TO 
<1)  ^ 
O 

>'  c 

-I 
2 1 
Lu  m 

TO  (U 
m  TO 
o  “■ 

=  I 

”  O 


UJ  re 
<i)  X 
^  c 

t  3 
3  CO 

-o 

c  ii 
re  re 
o  ^ 

.Z  (J) 
>*-  _ 

E  TO 

O)  g 

CO  < 


UJ 


"D 

_TO 

TO 

3 


C 

O 

Q-  'S 
</)  re 


■o 

TO 

o 

c 

re 

JZ 

c 

UJ 


c 

O 

4-> 

TO 

o 

“re 

TO 

■O 


I  I  I  I  I  I 


TO 

i_ 

3 

o 

(A 

o 

(A 

TO 

TO 

> 

0. 

c 

c 

o 

o 

'■3 

TO 

TO 

c 

c 

o 

o 

■*-> 

TO 

TO 

o 

Q 

TO 

TO 

O) 

O) 

X 

X 

• 

• 

24 


Handling,  Processing  Difficulties 
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ADVANCED  ENERGETICS 
Current  Materials 
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Encapsulated  Reactive  Materials 
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Submicron  Coated  Aluminum 


Nanoscale  Explosive  Composites 


(/) 

3 

(/) 

0) 


E 

3 

C 

E 

_3 

(U 

I 

o 

c 

(C 

c 


E 

c 

o 


o 

m 


E 

3 

C 

1 

_3 

I 

O 

c 

ra 

c 

X 

a> 


0) 

re 

E 

x 

o 

a 

Q. 

re 

0) 

N 

'w 

re 

o 

t 

re 

Q. 


re 

> 

'w 

o 

CL 

X 

re 

X 

Q 

0^ 


"O 

re 

(/) 


re  c 
=  °  = 

Q.  re  jp 
TJ  ^ 

c  w  ® 

<0  ^  55 

£  .2  o 

•g,t  0, 

I  a? 

>  re  o 

>%  N  w 

^  «  2 

”  S  - 

^  E  E 

T3  O  _ 

re  '3  re 
■D  E  o 


re  ^ 
o  3  a 
-I  (A  o 


re 

re 

.y 

'3 

re 

re 

■4-f 

re 

o 

re 

£ 

re 

o 

o 

§: 

o 

re 

c 

re 

re 

■a 

Ui 

re 

o 

E 

E 

4-« 

o 

re 

> 

i- 

o 

Q. 

E 

re 

D) 

re 

re 

ti 

re 

re 

o 

> 

o 

z 

o 

32 


Nanoscale  Explosive  Composites 


(/) 

^-1 

3 

(/) 

0) 


■o 

0) 

E 

o 

Q-  E 

(0  3 

c  .E 

i| 

^  (0 
(D 

a«? 
X  -r 
0) 

__  ^ 
(U  O 

5  ^ 

o  ' 


T3  £ 
■o  .ti 

<  ^ 


E 

3 

C 

E 

re 


I 


o 

CM 


I 


0} 

> 

w 

o 

a 

X 

0) 

X 

Q 

q: 


X 

o 
m  > 

I  I 

gl 


lO 


TJ  ^ 

Q>  .E 
E 

re  = 
o  = 
-I  re 


■o 

a> 

re 

o 

o  a, 

E  w 

II 

E 

3  re 
<  w 

re  u 
- 

re  O 
■E  0^ 


w  « 

re  ’«-• 
'•-*  _ 
(/)  o 

C  c 
o  o 


E  - 


(/> 

a> 

3 

a 

O 

o 

0) 

o 

c. 

(/) 

0) 

o 

G) 

O 

o 

E 

E 

o 

re 

c 

_o 

‘43 

■D 

a 

<D 

o 

re 

> 

<i> 

a 

Q. 

o 

re 

4-> 

o 

re 

re 

k. 

re 

a> 

re 

o 

c 

■  Hi 

o 

(/> 

T3 

O 


TJ 

C 

re 

re 

E 

4-» 

c 

c 

<D 

re 

> 

> 

O 

o 

(/) 

(A 

0) 

">3 

Ui 

c 

C 

re 

^3 

c 

■o 

re 

(U 

_> 

(/) 

o 

D 

(/) 

I  I 


33 


crystallization 


Nanoscale  Explosive  Composites 


(Z 

L. 

0) 

> 

0) 

V) 

■O 

o 

E 


0) 

u 

<u 


O 

(U 


Ui 

c 

w 

3 

(/) 

c 

0) 

E 

a> 

Q. 

X 

0) 

(D 

C 


E 

3 

C 

E 


3 

re 

X 

o 

re 


o 

X 

Q 

0^ 

*♦- 

O 


re 

> 

o 


E 


(D 


C 

0) 

> 

o 

(A 

■D 

re 

■O 

o 


c. 

re 

(/) 

(A 

C 

o 

o 

■D 

re 

c. 

're 

o 


(A 

_re 

u 

r 

re 

a 

■D 

re 

N 

'55 

c 

o 

t. 

o 

E 

I 

n 

3 

(A 


X 

re 

re 


o 

<N 


I 


o 


T3 

re 

re 

o 

o 

o 

jA 

< 


E 

3 

C 

E 

_3 

re 


X 

_re 

re 

o 

re 

re 

re 

re 

o 

re 

3 

(A 

TJ 

re 


T3 

re 

re 

o 

u 

o 

04 


I 


o 


T3 

C 

re 

E 

3 

C 

E 

_3 

re 


H 

LU 

CQ 

O) 

£ 

‘(A 

3 

X 

re 

re 


•  • 


34 


Nanoscale  Explosive  Composites 


35 


750  nm  500  nm  450  nm  400  nm  350  nm  300  nm  250  nm  200  nm  150  nm  100  nm  50  nm 


Nanoscale  Explosive  Composites 


(/) 

3 

(/) 

a> 

(Xi 


o 

■O 

TO 

■ 

■ 

o 

TO 

c 

Q. 

TO 

E 

1 

c 

o 

o 

c 

o 

o 

■O 

c 

TO 

C 

■D 

TO 

c 

TO 

■O 

O 

TO 

TO 

TO 

'(/) 

o 

TO 

O 

> 

O 

Q. 

(/) 

o 

(/) 

X 

TO 

TO 

cL 

5 

TO 

'4-» 

■ 

5 

X 

TO 

■D 

TO 

TO 

"O 

"O 

X 

O) 

TO 

■  ■■■ 

E 

C 

TO 

o 

>» 

O 

c 

TO  _ 

TO 

E 

2  E 

3 

C  3 

C 

TO  C 

^  ■ 

TO 

■  ^■1 

E 

o  E 

(H 

_3 

TO  3 

CO 

TO 

E  03 

(U 


03 

0) 

"O 


0) 
(U 

o 

(/) 

I 

O 

o  ^ 
a  c 

Q.  a> 
3  £ 
(/)  ♦- 

O  TO 

<2  ® 

(5 

^  ® 
TO  (fl 

iS  TO 

TO  TO 

^  </) 

2  TO 
TO  > 

TO  2. 

Q. 

'c  ^ 

^  o 


^  (/) 
^  TO 

^  5 


■D 

O 

TO 

O 

O 

c 

TO 


TO 

■  MBI 

TO 

TO 

E 

(/> 


</) 

Vi 

TO 


(0  L"  5 
TO  0  ^ 

5 

c  5 

c  o 

3  Q-  a 

“0)0 

E  c  c. 

3  -5  " 

^  OH 

E  P 

TO  ^ 

•i-ij  0  (/ 

S  ^  c 

V  I  T 

.12  I 

(/)  u.  Q 

O  TO 

0.2  c 

HI  !5  0 


36 


Nanoscale  Explosive  Composites 


(/} 

3 

(/) 

0) 

QC 

•4-* 

(U 

CO 

a> 

(U 

o 

CO 


(U 

E 

CO 


37 


Surface  Area  Advantages  and 


(0 

0) 

G) 

C 

o 

75 

O 


38 


Surface  area  =  O.lm^/g  Surface  area  =  74  m^/g 


Bomb  Calorimeter 


>• 

j’ 


lU 


■D 

0) 

« 

3 

O 

re 

o 

E 

-S’* 

re 

u 

0) 

3 

(A 

re 

0) 

£ 

E 


0) 

re 

■D 


o 

o 


00 

CO 

<D 


N. 

M 

CO 

(O 


O 

o 

CO 

00 

•6 

U) 

;d 

'5 

re 

u 

o 

N 

c 

0) 


05 

CM 

CO  00 

to 

CO 

CM  CM 

to 

00 

h-' 

r». 

N. 

CO  to 

to 

to 

05 

05 

CM  <M 

in 

in 

^  Tf 

h-  1^ 

00  1^ 
h- 

T-’  fsj 
00  00 


OJ  cs 


x»  O 
CN  tf) 
CNJ  Ol 


^  in 

0>  O) 
o>  O) 


N-  O 
^  O 


O 

(O 


(O 


o 

o 


O) 


re 

c 

o 

0) 


o 

o 

« 

a> 


E  c 

2  o 

C  (A 

E  5 

3  £ 
re  re 

^  P 
O  .fc 


o 

o 

in 

O) 


:=> 

E 

3 

C 

E 

3 

re 

X 

0) 

re 


00  in 
00  in 

in  (d 
o>  o> 


CM  CM 


00  oo 
T-  to 

T"  T— 

h- 


o  o 
o  o 


05  05 


E 

3 

C 

E 

3 

re 

in 

I 

X 


00  in 

O  CO 

to  to 

05  05 


CM  CM 
^  "O’ 
1^ 


CO  CO 

CO  in 


1^ 


o 

o 

CO 


05 

E 

3 

C 

£ 

_3 

re 

u 

re 

w 

0) 

U 


c 

o 


re 

JO 

re 

> 

< 


39 


Aveka/PVDF 
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Aveka  Coated  Nano-Particle 
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AVEKA:  Polymer-coated  Nano- 
Aluminum  Particles 
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feature  that  adds  energy 


X-ray  diffractometer 
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Future  challenges 
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Is  energy  transfer  to  the  working  fluid  efficient? 


NANOSCIENCE  MECHANICS  OF  MATERIALS 
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Abstract 

Predietive  eonstitutive  relations  are  to  be  obtained  for  the  deformation  and  fraeturing 
properties  of  materials  eontaining  nanometrie-sized  eonstituents.  Sueh  relations  are  of 
partieular  interest  for  nanoenergetie  materials.  The  approaeh  is  to  aeeomplish  this  goal 
by  extending  eurrent  understanding,  test  methods,  and  disloeation  meehanies  based 
models  that  have  been  established  for  predieting  the  behaviors  of  eounterpart  mieroseale 
materials. 


Introduction 

1 .  Dislocation  pile-up  based  Hall-Petch  (H-P)  considerations. 

The  dependence  of  the  plastic  yield,  flow  and  fracture  stresses  of  polycrystalline 
materials  on  the  reciprocal  square  root  of  their  average  grain  diameters,  as  interpreted  on 
a  dislocation  pile-up  (Hall-Petch)  basis,  was  extended  to  provide  a  model  understanding 
of  localized  adiabatic  (hot  spot)  heating  when  associated  with  discontinuous  load-drop 
behavior  in  impacted  materials,  including  energetic  materials  [1].  In  further  work,  it  was 
proposed  that  the  susceptibility  of  materials  to  such  localized  heating  in  shear  bands  was 
directly  proportional  to  the  stress  versus  inverse  square  root  of  grain  diameter  slope, 
known  as  the  H-P  “k”  value,  or  microstructural  stress  intensity  [2,  3].  The  stress  intensity 
label  was  applied  by  analogy  to  the  macroscopic  fracture  mechanics  stress  intensity  in 
that  the  mathematical  description  of  a  dislocation  pile-up  has  direct  connection  with  the 
description  of  a  continuum  mechanics  modeled  shear  crack  [4]. 

2.  Constitutive  equation  development. 

Recent  grain  size  dependent  H-P  type  results  and  model  analyses  have  been  based  on 
structural  material  properties  and  extended  to  consideration  of  the  influence  of 
temperature  and  strain  rate  effects  that  mostly  enter  into  the  average  friction  stress  for 
dislocation  movement  within  the  pile-ups  [5].  Maximum  load  condition,  dynamic 
recovery,  dynamic  strain  aging,  deformation  twinning,  and  true  fracture  strain 
considerations  have  been  described,  in  the  latter  case,  relating  to  particle  size  influences 
on  the  ductility  of  materials.  The  thermally  activated  equation  analysis  for  temperature 
and  strain  rate  effects  was  combined  with  the  dislocation  pile-up  model  to  describe 
mechanically-induced  “hot  spot”  behavior  in  energetic  material  decompositions.  An  H-P 
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type  prediction  was  obtained  for  initiation  of  explosion  of  RDX  (cyclotrimethylene- 
trinitramine)  crystals  in  drop-weight  impact  tests  [6]. 

3.  Elastic/plastic/cracking  hardness  stress-strain  analysis. 

Conventionally,  energetic  crystal  particles  of  sub-millimeter  to  micrometer  size  are 
employed  in  explosive  and  propellant  formulations.  The  individual  crystals  are  known 
to  be  relatively  brittle,  hence,  microindentation  hardness  testing  has  been  utilized  as  a 
convenient,  though  complicated,  quasi-static  method  of  evaluating  the  material 
deformation  properties,  much  as  employed  for  ceramic  materials.  Such  energetic  crystals 
are  found  to  be  relatively  hard  in  comparison  with  their  relatively  weak  molecular  or 
hydrogen  bonding.  The  occurrence  of  cracking  at  microhardness  impressions  in  RDX 
crystals  has  been  evaluated  on  an  indentation  fracture  mechanics  basis  and  compared 
with  parallel  measurements  on  MgO  crystals,  also  including  a  comparison  of  dislocation 
characteristics  [7].  More  recently,  these  material  hardness  properties  have  been  assessed 
on  an  indentation  hardness  stress-strain  basis  [8]  that  was  utilized  for  comparing  the 
relative  elastic,  plastic  and  cracking  behaviors  of  explosive  and  related  inert  crystal 
materials  [9].  Microhardness  testing  is  extended  recently  to  very  low  load 
“nanoindentation”  testing  for  which  the  proposed  hardness  stress-strain  method  of 
analysis  appears  to  be  useful  in  deciphering  the  meaning  of  results  [10].  Such 
nanoindentation  testing  is  being  qualitatively  applied  at  present  to  RDX  crystals 
examined  at  nanoscale  dimensions  with  the  Atomic  Force  Microscope  [11]. 

4.  Molecular  dislocation/lattice  modeling. 

An  early  investigation  of  comparative  dislocation  etch  pitting  at  residual 
microindentation  impressions  in  RDX  and  LiF  crystals  established  that  dislocation 
generation  and  movement  were  particularly  difficult  in  RDX  [12].  The  results  were 
explained  on  the  basis  of  a  special  resistance  to  shear  displacements,  even  within  the 
dislocation  displacement  field,  because  of  the  interlocked  nature  of  adjacent  molecules 
within  the  structure  of  the  RDX  orthorhombic  Pbca  unit  cell  [13].  Extension  to  modeling 
the  deformation  of  related  HMX  (cyclotetramethylenetetranitramine)  explosive,  with  a 
monoclinic  unit  cell,  led  to  explanation  of  twinning  being  preferred  in  HMX  because  of 
the  flexibility  of  the  larger  molecule  [14].  The  occurrence  of  nitroso  compounds  as  an 
early  decomposition  product  in  drop-weight  impact,  mechanical  sensitivity,  testing  of 
RDX  crystals  was  shown  to  be  an  expectation  of  forced  slip  deformation  within  the 
corresponding  lattice  structure  [15]. 


Current  research 

The  preceding  research  topics  are  being  pursued  for  application  to  understanding 
energetic  material  properties  and  the  influence  of  explosive  loading  on  dynamic 
behaviors  of  structural  materials.  The  topics  are  summarized  as:  (1)  dislocation  pile-up 
prediction  of  crystal  and/or  polycrystal  size  dependent  mechanical  behaviors;  (2) 
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thermally-activated  defect  rate  equation  analysis;  (3)  modeled  hardness  stress-strain 
description;  and  (4)  mechanical  property  interpretation  at  the  atomic/molecular  scale. 

1 .  Extension  to  nanometric-sized  energetic  materials. 

Evaluation  of  the  dislocation  self  energies  in  RDX  and  related  energetic  crystals,  relative 
to  surface  energy  considerations  and  by  comparison  with  a  dislocation-assisted  vapor 
growth  mechanism,  gives  indication  that  the  dislocation  cores  might  contain  included 
solvent  carried  over  from  conventional  crystal  growth  from  supersaturated  solutions  [16]. 
A  fuller  review  of  experimental  and  modeled  dislocation  properties  in  energetic  crystals 
is  in  progress  [17].  Of  particular  interest  is  the  extension  of  these  considerations  to 
nanoscale  crystal  sizes  where  the  surface  energy  plays  an  increasingly  important  role  in 
material  processing  behavior  [18].  Also  of  interest  in  connection  with  the  generation  of 
thermomechanical  stresses  are  the  relative  surface  energies  of  energetic  crystals  and 
related  polymer  binder  materials  that  are  employed  in  explosive  and  propellant 
formulations  [19]. 

2.  Extension  to  nano-indentation  and  -impact  testing. 

Collaboration  with  J.  Sharma  and  colleagues  [11]  has  involved  analysis  of  very  recent 
results  obtained  on  nanoindentations  put  into  an  RDX  crystal  using  an  atomic  force 
microscope  (AEM)  and  relating,  apparently,  to  microscale  indentation  effects  described 
for  MgO,  RDX  [7],  and  ammonium  perchlorate  (AP)  crystals  [20].  At  this  scale, 
individual  dislocation  generation  and  movements  should  be  able  to  be  evaluated  from 
AEM  surface  profile  measurements  that  are  in  progress.  Also,  the  measurements  are  to 
be  correlated  with  the  hardness  stress-strain  characteristics  modeled  on  the  microscale. 

In  addition,  collaboration  with  R.  Clough  has  been  recently  undertaken  to  consider 
development  of  a  dynamic  hardness  test  method  achieved  by  ball  impacts,  relating  to 
previous  analyses  for  usefully  assessing  the  influence  of  deformation  rate  on  material 
properties  [21].  Development  of  a  nano  indentation  test  facility  at  AERL/MNME  is 
planned. 

3.  Extension  to  composite  explosive/propellant  formulations. 

As  mentioned  above,  the  extension  of  dislocation  pile-up  based  H-P  model  considerations 
to  predict  a  material  sensitivity  to  shear  banding  behavior  shows  the  consideration  to  be 
important  for  energetic  materials.  Recent  application  of  the  same  type  of  pile-up  model 
consideration  has  been  to  predict  a  positive  H-P  dependence  for  the  indentation  fracture 
mechanics  stress  intensity  of  fine  grained  alumina  materials,  that  is,  grain  size  refinement 
provides  in  general  for  both  strengthening  a  material  and  improving  its  fracture  toughness 
[22].  The  question  of  how  far  such  considerations  may  be  extended  in  the  direction  of 
nanoscale  materials  continues  to  be  of  interest  [23]  and,  also,  the  extent  to  which  such 
considerations  may  be  applied  to  composite  materials  is  being  investigated  [24].  In  this 
regard  encouraging  results  have  been  obtained  on  extending  to  modeled  “flow  unit” 
displacements  in  polymeric  materials  the  thermal-activation  type  description  employed 
for  temperature  and  strain  rate  influences  in  polycrystalline  viscoplastic  materials  [25]. 
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An  interesting  eomparison  between  grain  size  versus  partiele  size  strengthening  of 
materials  is  that  grain  size  reduetion  potentially  eontributes  to  inereasing  material 
ductility  whereas  particle  strengthening  reduces  material  ductility  [26].  The  intention  is 
to  connect  research  activity  on  this  topic  with  an  award  that  has  been  received  recently 
for  an  AFOSR  International  Research  Initiative  (IRI)  effort  in  collaboration  with  the 
University  of  Cambridge  (UK)  on  the  dynamic  deformation  properties  of  composite 
materials  [27]. 
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ABSTRACT 


Previous  deseription  has  been  given  of  potential  advantages  and  ehallenges  assoeiated 
with  the  development  of  nano-energetic  materials  [1].  Such  nano  metric-scale  energetic 
materials,  as  well  as  counterpart  combustible  metals,  and  reactive  materials,  are  all  of 
current  research  interest,  for  one  reason,  because  surface-control  of  their  decomposition 
properties  provides  for  the  possibility  of  orders  of  magnitude  increase  in  energy  release 
rates  at  ultrafine  particle  sizes  [2].  Recent  report  of  enhanced  burn  rate  results  for 
aluminum  particle  combustion  [3]  appears  to  substantiate  the  prediction. 

The  surface  energy  property  is  of  paramount  importance  at  nanometric  particle  sizes,  as 
evidenced  by  nano-particle  processing  difficulties  with  agglomeration.  Particle  coating 
can  alleviate  the  problem,  especially  relating  to  coating  energetic  materials  with  polymers 
that  have  relatively  lower  surface  energies.  Otherwise,  the  added  surface  areas  of  smaller 
energetic  particles  provide  for  energy  enhancement  and  this  can  be  substantial,  say, 
compared  to  heats  of  fusion  among  different  materials.  At  the  smallest  particle  sizes,  or 
at  the  level  of  atomic  or  molecular  clusters,  issues  of  strained  bonds  at  the  material 
surfaces  and  size-dependent  surface  or  interfacial  energies  are  to  be  considered  [4].  The 
topic  lends  itself  to  molecular  dynamics  (MD)  modeling  as  described  for  the  simulated 
oxidation  of  aluminum  nano-clusters  [5]. 

Generally,  more  severe  conditions  of  phase  transformation  or  material  deformation  are 
required  to  produce  smaller  material  particles  and,  consequently,  added  internal  defects, 
say,  point  or  line  (dislocation)  defects,  might  be  expected  to  be  contained  within  the 
particles.  Eor  energetic  materials,  the  large  lattice-structure-determined  Burgers  vectors 
of  dislocations  provide  internal  strain  energies  that  so  outweigh  the  relative  surface 
energies  that  hollow  dislocation  cores  of  nm  dimensions  are  predicted  to  be  energetically 
feasible,  based  on  a  pioneering  model  description,  particularly,  for  vapor  grown  crystals 
[6]. 

The  mechanical  strength  and  fracturing  behaviors  of  nanometric-scale  material  particles 
and  structures  are  challenging  areas  of  research  so  far  as  concerns  a  dearth  of  material 
testing  methods  and  also  prediction  of  size-dependent  material  de formation/ decompo¬ 
sition  properties.  Extrapolations  to  nm  crystal  sizes  of  significantly  increased  drop- 
weight  height  measurements  required  for  mechanical  initiation  of  energetic  crystals 
remain  to  be  confirmed.  Recent  atomic  force  microscopy  (AFM)  observations  on  RDX 
(cyclotrimethylenetrinitramine)  crystals  have  indicated  the  presence  of  shock-generated 
nanoscale  shear  cracks  with  displacement  vectors  smaller  than  the  length  of  individual 
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molecules  [7].  More  recent  AFM  observations  have  included  preliminary  examination 
of  nanoindentation  deformations. 
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ABSTRACT 

Enhanced  burn  rale  lesults  are  presented  for  ammonium  perchlorate/AI  nanoparticle  strand  burners  at  atmospheric  (and  higher)  pressure  and 
lor  the  comparative  combustion  in  a  high  pressure  closed  vessel  of  a  solid  propellant  containing  15%  of  either  conventional  micrometer-scale 
Al  or  nanomelric  Al.  The  burn  rate  at  the  smallest  nanometric  Al  particle  size  appears  to  be  asymptotically  approaching  an  inverse  particle- 
diameler-squared  dependence. 


Surface-Controlled  Decomposition.  The  enhancements 
at  nanomelric  particle,  grain,  crystal,  or  other  unit  sizes  of 
either  structural  material  strength  properties  or  electronic 
dcvicc/circuil  performances  appear  to  he  relatively  mature 
research  topics  compared  to  the  achievement  of  greater 
energy  release  rates  al  similar  ulirafine  particle  dimensions. 
Such  potential  achievement  of  substantial  inercases  in  energy 
release  rates  for  the  chemical  decomposition  of  ullraftnc 
particles  might  be  foreseen,  however,  because  such  decom¬ 
positions  arc  often  surface  area  dependent,  as  well- 
demonstrated.  for  example,  in  the  use  of  ullmnne  particles 
for  catalytic  assistance  of  chemical  reactions.  Here,  we  arc 
concerned  with  the  achievement  of  enhanced  propellant  (and 
explosive)  energy  release  rales  al  nanomelric  particle  sizes, 
particularly,  involving  Al  panicles  whose  surface-controlled 
oxidation  to  alumina  is  strongly  exothermic.* 

Bum  Rate  Dependence  on  Al  Particle  Diameter.  Recent 
early  indications  of  increasing  the  bum  rate  of  solid  rocket 
propellant  materials  with  addition  of  Russian-produced 
“Alex"  nanoscalc  Al  powder  ingredient  has  been  credited 
by  Baschung  et  al.^  to  results  reported  by  Kuo  and  col¬ 
leagues.^  Baschung  ct  al.  made  their  own  comparison  al  ISL, 
France,  of  the  combustion  behaviors  exhibited  by  a  new 
double-base  gun  propellant,  al  different  pressures  in  a  high- 
pressure  closed  vessel  of  50  cm’  volume.  The  propellants 
employed  15%  of  Alex  material  suhsiiiuicd  for  the  same 
amount  of  conventional  micrometer-sized  Al  in  otherwise 
the  same  Pechincy  YX  76  propellant  formulation  of  55% 
nitr(x.'ellulnsc  and  .T0%  dioxyelhyinitramine  dinitnite  (DINA). 
Also  recently,  Technanogy,  LLC  has  prixluced  another  Al 
nanopartiele  material  and  obtained  experimental  suand-type 

*  Cum»p(Muiing  aulhor. 

»  AFRITMNMK 

<  Pronth-Gemuin  RcMruvh  InMitme  of  Saim  Louis 

*  Technanogy.  LLC 

10.1021/n)02S905k  CCC:  $25.00  C  2003  Amencan  Chemical  Society 
Published  on  Web  01/09/2003 


Figure  1.  Bum  rate  as  a  function  of  aluminum  particle  size,  also 
showing  indicated  pressure  ranges  for  the  vertical  lines,  and 
providing  indication  at  the  smallest  particle  sizes  of  an  inverse 
square  dependence  of  bum  rale  on  particle  diameter. 

bum  rale  measurements,  both  as  a  function  of  particle 
diameicr  at  ambicni  pressure  and.  separately,  for  200  nm 
particle  diameter  material  burned  at  different  higher  pres- 
surcs. 

First,  Figure  1  shows  a  comparison  of  the  combined  Alex 
and  Technanogy  results  in  terms  of  the  bum  rate  dependence 
on  particle  diameter.  The  aluminum  material  employed  in 
the  ISL  experimenis  had  a  range  of  panicle  sizes:  for  the 
conventional  Al  between  10  and  100  micrometers;  and.  for 
the  Alex  Al  between  50  and  200  nm.  As  noted  by  the  vertical 
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line  segments  at  each  panicle  si/e.  pressures  between  4()  and 
200  MPa  were  realized  in  the  closed  chamber  for  these  bum 
rates  that  were  monitored  in  terms  of  the  dynamic  vivacity 
property  ol  the  propellants.  The  dynamic  vivacity  L  is  denned 
as  the  derivative  of  the  pressure  dp/dr  divided  by  the  product 
of  the  pressure  p  times  the  maximum  pressure  achieved 
in  the  experiment,  as  the  expression  (.STANAG  41 15): 

L  =  dp/dr  X  1/p  X  1/P^, 

Monitoring  the  vivacity  during  the  pressure  histories 
provided  an  account  of  the  change  in  burning  behavior  during 
consumption  of  the  material  mass.  For  the  ISL  comparison 
of  the  two  different  sized  A1  ingredients,  emphasis  was  given 
to  the  Alex-containing  propellant  exhibiting  a  pressure- 
dependent  steep  increase  in  vivacity  at  the  beginning  of 
combustion.  ITte  initial  burn  rate  behavior  was  found  to 
depend  only  on  pressure  and  not  on  loading  density, 
apparently  eliminating  such  mechanical  effect  as  porosity 
from  playing  a  major  role  in  determining  the  higher  A1 
nanoparticle  burn  rate  property. - 

Also  shown  in  Figure  I  arc  the  Technanogy  strand-type 
burn  rate  measurements  obtained  for  different  Al  particle 
diameter  ingredients  within  a  standard  ammonium  perchlo¬ 
rate  formulation.  The  strands  were  3.175  mm  in  diameter 
with  88%  solids  and  of  6.35—8.47  mm  length  pres.sed  to 
80-85%  density.  Here,  as  indicated,  the  particle  diameters, 
found  by  high  resolution  scanning  electron  microscopy  to 
be  essentially  spherical,  were  ciuefully  sieved  to  the  indicated 
particle  diameters  and  the  bum  rate  measuremenLs  were  made 
at  ambient  pressure.  As  shown,  substantial  enhancement  of 
the  bum  rate  occurred  in  the  particle  diameter  range  between 
200  and  40  nm.  for  example,  giving  an  indicated  exponential 
slope  dependence  of  {—2}  marked  near  to  the  smallest 
particle  burn  rate  measurement.  For  this  Al  material  a 
nominal  oxide  thickness  of  2  nm  applies  for  the  different 
particle  sizes.  Similar  bum  rate  measurements  have  been 
made  with  strands  containing  greater  solids  content,  even 
approaching  100%.  and  with  other  oxidizer  ingredients.''  For 
the  Technanogy  Al  product,  separate  pressure-dependent 
measurements  were  made  also  with  200  nm  Al  particles.  In 
Figure  I,  the  vertical  line  segment  is  shown  to  cover  an 
incrca.se  in  bum  rate  from  ~  12-580  mm/s  for  a  pressure 
incrca.se  from  0.14  to  15  MPa.  respectively,  and  shown  at 
the  higher  pressure,  to  give  a  greater  bum  rate  than  achieved 
with  the  ISL  gun  propellant. 

Burn  Ratc/Pressure  Dependencies.  The  pre.ssure.  P. 
dependence  of  the  propellant  bum  rate.  r.  is  normally  of 
interest  for  de.scription  via  a  power  law  dependence,  that  is. 
in  accordance  with  the  so-called  Vieille’s  law,  expressed  a.s¬ 


r  =  /jP“ 

in  which  /?  and  a  arc  experimental  constants.  The  general 
form  of  the  power  law  dependence  has  a  theoretical  inter¬ 
pretation.*  In  particular,  small  values  of  the  pressure  expo¬ 
nent.  a,  are  aimed  for.  in  order  to  avoid  transition  to  unstable 
burning  at  high  pressures.  Figure  2  shows  a  comparison  of 
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Pressure  (MPa) 

Figure  2.  Pressure  dependence  of  the  bum  rate  for  200  nm  panicle 
diameter  aluminum  and  for  conventional  and  Alex  containing 
aluminum  gun  propellant. 


the  pressure  dependent  results  obtained  by  Bu.schung  et  al. 
for  both  gun  propellant  materials  in  comparison  with  results 
obtained  in  the  Technanogy  experiments. 

A  smaller  (exponential)  slope,  a,  for  the  Alex-containing 
gun  propellant  of  0.66  was  deduced  by  Baschung  ct  al..  in 
comparison  with  a  equal  to  0.8  for  the  conventional  Al  sized 
propellant,  and  an  enhanced  p  value  of  9.85  mtn/s  was 
obtained  for  the  new  formulation  with  Alex  material.  For 
the  Technanogy  product,  the  power  law  dependence  is 
obeyed  strictly  only  at  the  intermediate  pressures  between 
1.1  and  11.3  MPa.  giving  an  a  value  near  to  that  of  the 
conventional  gun  propellant  along  with  an  appreciably 
increased  P  value  of  25  mm/s.  Otherwise,  an  average  long- 
dashed  line  of  larger  p  and  smaller  a  is  shown  in  Figure  2 
in  order  to  adjust  to  the  lowest  pressure  measurements.  The 
lowest  pressure  measuremenLs.  which  are  high  relative  to 
the  Vicille  law  extrapolation,  may  be  compared  to  the 
separately  determined  particle  diameter  dependence  measure¬ 
ments  shown  in  Figure  1.  More  interestingly,  the  highest 
measured  bum  rate  583  mm/s  at  15  MPa  in  Figure  2  is  shown 
to  be  significantly  above  the  Vicille  dependence  for  the  lower 
pressure  measurements,  thus  being  indicative  of  unstable 
burning. 

.Model  Consideration.  A  bum  rate  dependent  on  the 
inverse  square  of  the  Al  particle  diameter  should  he  expected 
at  very  small  particle  sizes  because  of  the  energy  release 
rate  being  determined  by  the  available  amount  of  particle 
surface  area,  say.  specified  by  the  number  of  Al  particles 
per  unit  surface  area.  The  surface  area  change  with  particle 
diameter  is  innocuous  at  large  particle  diameters  where 
continuation  of  individual  particle  burning  is  controlled  more 
so  by  further  diffusion  of  Al  atoms  through  the  coherent 
alumina  “.skin".  Such  inverse  square  of  the  nanoparticle  size 
dependence  compares,  for  example,  to  a  lesser  particle- 
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iliamctcr-dclcrmincd  bum  rale  controlled  by  the  amount  of 
particle  surface  area  per  unit  volume  (within  the  bum  /.one), 
then  corresponding  to  a  reciprocal  dependence  of  the  burn 
rate  on  particle  diameter.  A  higher  bum  rate  dependence  even 
than  the  inverse  square  dependence  would  be  determined 
for  control  by  the  number  of  particles  per  unit  volume, 
corresponding  to  a  reciprocal  cube  dependence  on  the  particle 
diameter.  For  any  of  these  proposed  bases,  the  current  results 
give  pnrmising  indication  of  orderfs)  of  magnitude  enhance¬ 
ment  of  energy  release  rates  in  the  new  regime  of  nanometric 
particle  ingredients,  and  current  experimcnts/analyses  are  in 
progress  in  pursuit  of  firmer  grounding  of  the  potential 
enhancements.' 

Acknowledgment.  Appreciation  is  expressed  to  Byron 
Atimon.  at  AFRL/MNME,  for  construction  of  figures  and 
assistance  with  the  manuscript. 

Note  .Added  in  Proof:  Recent  contact  with  Professor 
Dana  DIott.  University  of  Illinois,  has  led  to  receipt  of  a 
proof  copy  of  the  coming  Chemical  Physics  Letters  article 
“Fast  spectroscopy  of  energy  release  in  nanomcu-ic  explo¬ 


sives”.  co-authored  by  .Shufeng  Wang.  Yanqiang  Yang. 
Zhaoyong  Sun.  and  Dana  D.  DIott,  in  which  important  laser- 
induced  measurements  arc  described  of  nanosecond-scale 
energy  release  rate  characteristics  for  a  dilute  1%  Alex/ 
nitrocellulose  system. 
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GRAIN  SIZE  DEPENDENCE  OF  SHOCK  INDUCED  TWINNING 
STRESSES  FOR  ARMCO  IRON 
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Measurements  fitted  to  a  linear  inverse  square  root  of  grain  size  dependence  for  the  deformation-twinning- 
deteimined  Hugoniot  Elastic  Limit  (HEL)  stresses  of  shocked  Armco  iron  specimens  of  differing  average  grain  sizes 
are  shown  to  compare  favorably  with  a  Hall-Petch  (H-P)  dependence  established  for  the  combination  of  split 
Hopkinson  pressure  bar  (SHPB)  and  conventional  deformation  test  results  known  to  be  controlled  by  deformation 
twinning.  The  comparison  of  shock  test  results,  in  the  pressure  range  of  ~0.6  to  ~-2.0  GPa,  are  suggested  to  give  a 
counterpart  match  of  flEL  and  conventional  twinning  stress  values  analogous  to  the  calibration  of  reference  static 
and  shock  pressure  levels  required  for  the  a  -  e  phase  transition,  in  the  latter  case,  corresponding  to  a  volumetrically- 
better-specified  pressure  of  1 3.5  GPa.  A  further  comparison  of  the  essentially  athermal  twinning  stresses  with  those 
computed  for  plastic  flow  on  a  dislocation  mechanics  model  basis  provides  for  a  grain  size  dependent  transition 
between  twinning-  or  slip-determined  plastic  flow  stresses. 


INFRODIJCTTON 

At  higher  strain  rates  and  larger  polycrystal  grain 
diameters,  defonnation  twinning  is  an  increasingly 
likely  defonnation  mechanism  in  structural  metals, 
particularly,  body  centered  cubic  (bcc)  metals  and 
alloys  (1).  Current  interest  centers  on  the  role  of 
twinning  under  shock  loading  conditions  (2).  In  the 
present  report,  pioneering  laboratory  test  results  (3-5) 
and  more  recent  shock  test  results  (6,7)  obtained  on 
deformation  twinning  in  low  carbon  and  Armco  iron  are 
brought  together  on  a  Hall-Petch  (H-P)  stress  versus 
inverse  square  root  of  average  grain  diameter  basis  to 
give  a  rational  explanation  of  shock-induced 
deformation  twinning  and  its  relationship  to  alternative 
deformation  by  dislocation  mechanics  based 
viscoplastic  slip.  The  compared  deformation  behaviors 
relate  also  to  particular  solid  cylinder  impaet  (Taylor) 
test  results  where  essentially  athermal  twinning  occurs 
on  first  impact  and  is  followed,  on  a  microsecond  time 
scale,  by  later  initiation  of  dislocation  slip  within  the 
twin-hardened  material  (8). 


TWINNING  VS.  SLIP 

Figure  1  shows,  first,  addition  of  shock-related 
strength  measurements  for  specimens  of  different 
average  grain  diameters  (7)  to  a  previous  flall-Petch  (I  I- 
P)  based  compilation  (8)  of  experimental  and 
theoretical  slip  and  twinning  stresses  collected  for 


annealed  Armco  iron  materials  and,  in  one  case, 
annealed  SAE  1010  steel  material,  as  obtained  too  on 
test  specimens  with  different  polycrystal  grain  sizes.  In 
the  Figure,  the  Armco  iron  twinning  stress  points  listed 
for  the  five  different  temperatures  in  the  range  between 
88  and  248  K  [3],  were  reported  separately  at  each 
temperature  in  the  original  work  as  proportional  limit 
(p.l.)  stresses,  as  measured  in  split-Hopkinson  pressure 
bar  tests  at  a  strain  rate  of  I .  I  x  1 0^  s  ‘. 

On  a  combined  data  basis,  the  different  temperature 
dependent  twinning  results  show  that  preference  for 
twinning  over  slip  was  restricted  to  lower  temperatures 
for  smaller  polycrystal  grain  sizes.  That  is,  the 
temperature-dependent  transition  from  an  athermal 
twinning-controlled  yield  stress  to  athennally-activated 
slip-controlled  one  moves  to  lower  temperatures  for 
smaller  grain  size  materials.  In  Figure  1 ,  the  relatively 
high  stress,  cross-within-circle,  datum  at  large  grain 
size,  was  obtained  by  Rohde  (4)  in  a  shock-loading  test 
where  extensive  twinning  occurred.  The  1010  steel 
measurement  in  the  Figure,  designated  MWA  (5)  and 
spanning  a  vertieal  range  in  stress  between  proportional 
limit  (p.l.)  and  load  drop  (l.d.)  stresses,  was  obtained  by 
compression  testing  at  a  conventional  laboratory  strain 
rate  of2.8xl0'*s''  and  at  liquid  helium  temperature,  4.2 
K..  The  higher  (l.d.)  stress  value  was  determined  for 
pronounced,  twin-produced,  load  drops  that 
accompanied  yielding  after  further  stress  loading  past 
the  lower  p.l.  stress.  This  higher  raised  stress  for 
general  yielding,  after  the  p.l.,  was  employed  in  relation 
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to  all  of  the  p.l.  data  then  for  designation  of  the  H-P 
linear  stress  dependence  shown  for  the  twinning  yield 
stress  and  expressed  in  the  H-P  relation: 

ot  =  Oto  +  kr  E  '  ^  ( I ) 

where  Oto  =  330  MPa  and  kt  =  90  MPa.mm'^. 
Equation  ( I ),  with  these  constants,  was  employed  in 
model  calculations  that  predicted  satisfactorily  the 
involvement  of  twinning  in  cylinder  impact  (Taylor) 
test  results  (8). 

A1.SO  shown  in  Figure  1  are  several  computed  11-P 
dependencies  for  dislocation  mechanics  based  model 


computations  of  the  so-called  Z-A  slip-controlled 
plastic  yield  stress  at  different  strain  rates,  expressed  as 

Oj  =  0(1  +  kC''^  (2) 

where  the  slope,  k,  though  being  smaller  than  kj,  is  also 
athennal,  and  Oo  is  given  by 

0,1  =  0,1  +  B  e  '*^  (3) 

where 

P  =  P,i  -  Pi  In  (s  dot)  (4) 


Figure  I .  Comparison  of  experimental  shock  test  results  with  previous  high  strain  rate  or  low  temperature  laboratoiy  test 
measurements,  including  computed  viscoplastic  slip-controlled  yield  stress  results  at  different  strain  rates  and,  also, 
converted  reference  hardness  flow  stress  values. 


with  material  constants  previously  reported  for  Armco 
iron  (9)  as  Oc  ~  0,  B  =  1033  MPa,  =  0.00698  K  ', 
and  Pi  =  0.000415  K  '.  The  cross-within-circle  yield 
stress  point  (4),  mentioned  earlier  in  Figure  I,  for 
example,  applies  tor  twinning  detected  in  the  shock 
wave  test  at  an  estimated  strain  rate  of  I O'  s  ',  that  is  in 
agreement  otherwise  with  the  somewhat  higher  stress 


predicted  for  slip-controlled  yielding  at  that  same  strain 
rate. 

The  slip-controlled  yield  stress  M-P  lines  at  different 
strain  rates  in  Figure  1,  with  k  =  22  MPa.mm' give 
increasing  yield  stresses  that  are  rai.sed  non-linearly 
from  the  lowest  conventional  laboratory  strain  rate  of 
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0.014  s  '  to  the  indicated  higher  strain  rates,  all  in 
accordance  with  the  thermal  dependence  in  Cd  in  eqs.  3 
and  4. 

The  second  lowest  H-P  line  in  Figure  1,  however, 
was  computed  for  the  conversion  of  Vickers  HVIO 
hardness  measurements  that  were  made  initially  on  the 
annealed  Amico  iron  material  prior  to  shocking  (7).  In 
this  case,  with  multiplication  of  the  hardness  pressure 
by  a  factor  of  (1/3),  the  measurements,  made  at  1 0  kgf 
load,  are  converted  to  estimated  compressive  flow  stress 
values  corresponding  to  a  true  strain  of  0.075. 

The  predicted  increase  in  flow  stress  shown  for  the 
second-lowest  hardness-determined  H-P  line  in  Figure 
1  was  obtained  by  accounting  for  the  strain  hardening  at 
0.075  strain  past  yielding  with  the  stress  increment 

Ao=Ke"  (5) 

where,  for  Armco  iron  (9),  the  experimental  constants 
are  obtained  as  K  =  266  MPa  and  n  =  0.289.  As 
indicated  in  the  Figure,  a  computed  increase  in  stress 
value  of  126  MPa  is  predicted  for  the  raised  H-P  flow 
stress  line,  that  is  somewhat  above  the  converted 
hardness  points  but  otherwise  provides  a  reasonable 
basis  for  relating  the  conventionally  assessed  plastic 
flow  measurements  for  the  pre-shocked  material  and 
those  results  to  be  described  now  tor  shock  testing. 

Relatively  higher,  twin-controlled,  both  cross- 
within-circle-,  -within-triangle-,  and  four  -within- 
diamond-identified  yield  stress  values  are  shown  in 
Figure  1.  as  determined  from  Hugoniot  elastic  limit 
(HHL)  measurements  obtained  in  shock-type,  plate 
impact,  tests  (6,7,10).  The  HEL  measurements  have 
been  transformed  to  unidirectional  yield  stress 
measurements  according  to  the  relationship 

a  =  [(l-2u)/(l-u)]0HEL  (6) 

with  u  =  0.28,  being  Poisson’s  ratio. 

Despite  the  observation  being  made  that  there  is 
significant  scatter  from  the  linear  dashed  ll-P 
dependence  shown  for  the  six  converted  HEL-type 
measurements  in  Figure  I ,  the  measurements  are  taken 
anyway  to  give  the  suggested  greater  ll-P  line 
dependence  that  is  indicated.  There  are  reasons  for 
expected  test  variations.  The  HEL  is  not  easily 
determined  experimentally,  say,  relative  to  the  better 
volumetrically-determined  a  -  £  phase  transformation 
in  iron,  also,  occurring  at  a  greater  pressure  of  13.5 
GPa.  And,  in  Figure  1,  for  example,  the  largest  grain 
size  converted  HEL  measurement,  that  is  shown  to  fall 


below  the  dashed  H-P  twinning  line  for  this  set  of  data, 
is  very  probably  accounted  for  by  stress  lowering 
associated  with  a  reduced  number  of  grains  in  the 
specimen  cross-section  (7)  and,  therefore,  a  reduction  in 
true  polycrystalline  constraint  associated  with  restricted 
accommodation  of  the  twinning  strain  (1 ). 

THE  TWINNING  kj 

The  H-P  dependencies  shown  in  Figure  I  for  the  slip 
and  twinning  stresses  are  fully  consistent  with  a  model 
dislocation  pile-up  explanation  for  the  relation, 
including  explanation  of  the  slope,  kx,  being  a 
microstructural  stress  intensity  that  is  greater  for 
nucleation  oftwinning  at  grain  boundaries  as  compared 
with  the  smaller  k  for  the  propagation  of  slip  ( I ).  A 
dislocation  pile-up  based  model  de.scription  was  given 
earlier  for  the  nucleation  of  twinning  (II),  with 
parametric  evaluation  of  kj  being  logarithmically 
dependent  on  the  strain  rate  raised  to  a  small  power. 
And  numerical  estimates  were  given  previously  also  for 
kx  depending  only  weakly  on  temperature  and  even 
more  weakly  on  the  strain  rate. 

Arnold  reported  elastic  strain  rate  values  in  the  range 
of  I  O'-to- 1  o'’  s  '  for  his  HEL  test  results,  except  for  the 
one  lower  measurement  for  the  largest  grain  size  test 
(7).  According  to  the  weak  twinning  dependence 
derived  from  (II),  mentioned  above,  an  increase  in  kx 
between  9  and  1 4%  of  the  conventionally-determined 
value  is  estimated.  The  predicted  increase  at  the  shock- 
type  strain  rates  for  kx,  although  large  enough  to  be 
experimentally  measured,  is  smaller  than  that  indicated 
by  the  da,shed  line  for  kx  =  130  MPa.mm"^  from  the 
shock  wave  measurements. 

Such  consideration  of  the  kx  value  leads  to  an 
available  alternative  upper  limit  for  such  H-P  k  values, 
say,  by  comparison  with  theoretical  and  experimental 
upper-limiting  k  values  for  cleavage.  The  tensile 
cleavage  value,  kf,  is  estimated  in  one  way  by  the 
circular  pilc-up  stress  intensity  for  dislocation 
coalescence  (12)  as 

kc  =  (mn/4a)  Gb''^  (7) 

where  m  is  an  orientation  factor,  a  =  2(  I  -v)/(2-v),  G  is 
the  shear  modulus  and  b  is  the  dislocation  Burgers 
vector.  With  m  =  2,  a  value  of  kc  ~  80  MPa.mm'^  is 
obtained  (12)  that  compares  with  the  kx  value  of  90 
MPa.mm'^  from  Figure  I  and  with  experimental 
cleavage  values  (13,14)  between  90  and  120 
MPa.mm''^.  No  cracking  was  found  in  the  shock- 
loaded  specimens  of  Arnold  at  the  HEL  and  so  the 
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drawn  dashed  estimation  ofkT  in  Figure  I  is  very 
probably  an  overestimate  of  the  true  kt  value. 

CYLINDER  IMPACT  RESULTS 

A  major  application  of  the  separate  constitutive 
equations  ( I )  and  (2),  for  tw  inning  versus  slip  in  Amico 
iron,  has  been  to  prediction  of  solid  cylinder  (Taylor) 
impact  test  results,  for  example,  reported  by  Johnson 
and  Cook  ( 1 5),  and  modeled  with  their  Lagrangian 
(Elastic  Plastic  Impact  Computation)  EPIC  material 
dynamics  code  (8,9).  With  the  deformation  twinning 


stress  in  eqn  (I )  taken  as  essentially  athermal,  then 
twinning  should  definitely  occur  on  first  impact,  if  the 
stress  is  high  enough,  and  should  be  followed  then  by 
viscoplastic  slip  in  the  “new”  twin  hardened  structure. 
McKirgan  (16)  has  verified  the  sequential  nature  ofthe 
overall  deformation  process  by  step-wise  model 
assessment  of  the  Taylor-type  deformation  state. 

More  recently,  T aylor  test  results  have  been  reported 
for  the  same  shock-tested  Annco  iron  material  as 


Figure  2.  Cylinder  impact  (Taylor)  test  results  obtained  at  300  or  600  K,  in  comparison  with  twin-  or  slip-controlled 
yield  stress  relations. 


employed  for  the  multiple  HEL  measurements  shown 
in  Figure  I .  The  Annco  iron  material  has  now  been 
tested  at  the  two  temperatures  of  300  and  600  K  as 
shown  by  the  measured  points  labeled  NHA  (1 7)  and 


plotted  in  Figure  2.  In  this  case,  comparison  is  also 
made  at  300  K.  with  the  twinning  and  high  strain  rate 
slip  determined  FI-P  dependencies  taken  from  Figure  I . 
An  additional  computed  H-P  line  has  been  put  onto  the 
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Figure  for  600  K  and  a  strain  rate  of  10^  s  ',  that  is 
appropriate  also  to  such  Taylor  test  results  (9). 

An  apparent  curvature  in  an  H-P  dependence  is 
indicated  in  Figure  2  for  the  total  VISAR-determined 
yield  stress  points.  However,  the  superposition  of  the 
separate  H-P  dependencies  for  twinning  and  slip 
indicates  that  the  curvature  is  more  likely  a 
consequence  of  twinning  being  the  controlling 
mechanism  foryielding  at  the  largest  grain  size  and  slip, 
at  the  smallest  grain  sizes.  Such  interpretation  is  in 
line  with  metallographic  observations  made  on  the 
impacted  specimens  (1 7).  Again,  the  large  grain  points, 
although  showing  the  yield  stress  to  be  essentially 
athermal,  are  relatively  low  because  of  the  lack  of  full 
polycrystal  constraint  associated  with  few  grains  being 
present  in  the  specimen  cross-section.  Further  support 
for  this  interpretation  comes  from  the  reported  remark 
that  the  deformation  of  individual  grains  showed  up  in 
the  velocity  curve  records  of  the  large  grained 
specimens,  thus  connecting  with  the  earlier  explanation 
of  the  large  grain  HF.L  result  in  Figure  I .  An  important 
metallographic  observation  on  the  deformed  specimens 
led  to  the  conclusion  that  deformation  twinning  had 
preceded  follow-on  slipdefonnation  on  the  basis  of  the 
distorted  shapes  of  the  otherwise  expected  straight 
deformation  twin  boundaries. 

Finally,  such  consideration  as  described  here  of 
added  twinning  complication  to  deciphering  the  shock- 
induced  deformation  properties  of  structural  materials  is 
but  one  of  several  challenges  necessarily  involved  with 
the  special  regime  of  shock-induced  deformations.  For 
example,  the  very  .substantial  displacements  occurring 
at  nanometric  dimensions  at  all  points  along  the  shock 
front  are  not,  in  general,  capable  of  being  relaxed  by 
remote  displacement  of  the  more  widely-spread, 
initially  present,  resident  dislocation  population. 
Instead,  even  without  twinning,  a  nano-scale  dislocation 
structure  is  produced  sequentially  at  the  propagating 
shock  front,  and  left  behind  for  interaction  with  the 
lagging  movement  of  the  originally  resident  dislocation 
structure  (18,19).  Thus,  another  additional 
consideration  of  a  “changed  material”  to  be  dealt  with 
in  describing  subsequent  viscopla.stic  flow  obtains.  An 
interesting  consequence  of  such  consideration  is  the 
description  of  simultaneous  increases  in  both  strength 
and  ductility  under  particular  high  strain  rate  conditions 
(20). 


SUMMARY 


A  Hall-Petch  dependence  is  indicated  for  the  shock 
wave  determined  HEL  yield  stress  of  Armco  iron 
material  whose  shock  response  is  controlled  by 
defomiation  twinning.  Under  such  high  strain  rates,  in 
the  range  of  I  O'  to  10'’  s  ',  as  occur  in  shock  loading 
compared  to  conventional  and  split  Hopkinson  pressure 
bar  (SHPB)  tests,  the  further  indication  is  that  a 
somewhat  greater  I  l-P  microstructural  twinning  stress 
intensity,  kt,  value  is  obtained. 

The  relatively  different  H-P  dependence  fortwinning 
as  compared  with  slip,  when  applied  in  constitutive 
equation  descriptions  of  the  two  behaviors, 
demonstrates  that  twinning  is  favored  at  larger  grain 
sizes  in  both  shock  loading  and  cylinder  impact 
(Taylor)  tests.  During  impact  in  a  Taylor  test,  twinning 
occurs  before  slip,  if  the  stress  level  is  sufficiently  high, 
and  then  further  deformation  is  achieved  by  viscoplastic 
slip  in  the  newly  twinned  (and  hardened)  material. 
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Deformation  shear  influences,  rather  than  hydrostatic  pressure  considerations,  are  described  for 
impact/shock  initiation  behaviors  of  energetic  crystals  in  terms  of  differentiating  between  discrete 
dislocation  and  continuum  mechanics  descriptions  of  shear  banding  properties.  The  dislocation  mechanics 
view  extends  downward  to  inter-molecular  interactions  at  the  crystal  lattice  scale  and  includes  the  proposed 
creation  of  nanometer-scale  dislocation  stmctures  at  sufficiently  strong  shock  fronts.  An  initial  comparison 
is  given  between  dislocation  pile-up  avalanche  and  pore  collapse  models  for  impact/shock  initiations. 


INTRODUCTION 

Modem  research  study  has  been  directed  to 
assessing  the  relative  importance  of  shear 
stress/strain  behaviors  of  solid  energetic  material 
formulations  as  compared  with  hydrostatic  pressure 
influences  on  initiation  behaviors.  Not  so  long  ago, 
Eyring  (1)  pointed  out  in  discussion,  presumably  on 
an  hydrostatic  pressure  basis,  that  initiation  in  the 
solid  state  must  involve  an  effective  space 
associated  with  holes,  dislocations,  or  the 


separations  of  grains  where  the  breaking  up  of 
a(n  energetic)  molecule  can  occur  with  respect 
to  the  same  strong  bond  in  the  gas  phase.  In 
this  same  vein,  Gilman  and  Armstrong  (2)  gave 
emphasis  in  developing  model  results  for  the 
shear-induced  polymerization  of  benzene  to  the 
pioneering  “mechanochemistry”  experiments  of 
Bridgman  demonstrating  that  shear  stresses 
combined  with  compression  enhanced  the  onset 
of  chemical  reaction. 


Figure  1.  Particle  size  dependence  of  model  piston-driven  impacted  beds  of  HMX  particles,  after  computed 
results  presented  by  Gonthier  (5). 
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In  the  current  report,  the  mechanism  of  hot  spot 
development  accomplished  via  internal  stress 
concentrations  leading  even  to  crystal  fracture 
associated  release  of  a  dislocation  pile-up  avalanche 
(3)  is  considered  in  further  detail,  including  certain 
comparisons  with  a  continuum  mechanics  model 
view  of  shear  banding.  Connection  of  the 
dislocation  mechanics  model  with  crystal-lattice- 
scale  considerations  of  molecular  stability  and  bond 
breaking  is  described.  And,  lastly,  preliminary 
comparison  is  made  between  the  shear-based 
dislocation  avalanche  mechanism  and  a  pressure- 
based  pore  collapse  mechanism  for  hot  spot 
generation. 


MODEL  RESULTS 
Particle  Size  Influences 

The  dislocation  pile-up  avalanche  model  has 
been  characterized  by  a  parametric  inverse 
square  root  of  crystal  or  grain  size  dependence 
for  the  drop-weight  impact  sensitivity  property. 
Supporting  drop-weight  impact  sensitivity 
results  were  presented  (4)  for  RDX  (cyclotri- 
methylenetrinitramine).  Additional  crystal  size 
dependent  results  for  RDX  and  a  related 
energetic  crystal  structure,  CL- 12,  were 
similarly  presented  later  (5). 


Log  (Ax,).  Ax,  in  nm 

Figure  2.  Theoretical  pile-up  avalanche  temperatures,  at  two  crystal  sizes  for  RDX  and  PETN,  in 
comparison  with  their  previously  computed  thermal  explosion  temperatures. 
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In  Figure  1,  a  similar  model  interpretation  to  that 
mentioned  for  (4,5),  is  given  for  somewhat  distantly 
related  test  results  analyzed  more  recently  by 
Gonthier  (6).  In  the  Figure  are  plotted  computed 
results  for  impacting  piston  velocity  versus 
reciprocal  square  root  of  particle  diameter  for 
different  compacts  of  HMX  (cyclotetra- 
methylenetetranitramine)  particle  beds  and  with 
demarcation  indicating  that  the  beds  with  larger 
particle  diameters  were  being  impacted  at  lesser 
velocity  than  needed  to  sustain  combustion  within 
the  bed.  The  computations  were  based  on  earlier 
reported  experimental  results  obtained  by  Jacobs, 
Sandusky  and  Elban  (7). 

Figure  1  has  been  constructed  on  the  following 
basis:  (1)  the  Hertzian-type  contact  stresses 
between  spherical  particles  in  the  bed  correspond  to 
those  required  for  full  yielding  of  the  particles  and 
thus  follow  an  inverse  square  root  of  particle  size 
dependence,  as  for  the  impact  sensitivity  prediction 
mentioned  above;  and  (2)  at  high  impacting 
velocity  (relating  to  the  material  strain  rate),  the 
material  flow  stress  depends  simply  on  the 
logarithm  of  the  strain  rate.  Consequently,  the 
graph  is  proposed  to  represent,  at  small  particle 
sizes,  an  effective  yield  stress  dependence  on  the 
reciprocal  square  root  of  particle  diameter.  At  the 
largest  particle  diameters  in  Figure  1,  the  condition 
of  general  yield  has  not  been  obtained  and 
combustion  is  not  sustained. 

Lattice/Molecule  Stabilities 

Comparison  is  made  in  Figure  2  between  the 
predicted  pile-up  avalanche  initiated  hot  spot 
properties  of  RDX  and  PETN  (pentaerythritol 
tetranitrate)  crystals,  relating  also  to  previous  model 
calculations  given  by  Boddington  (8)  for  thermal 
explosion. 

The  pile-up  avalanche  temperatures  are  computed 
on  an  upper  limiting  temperature  basis  (3) 
stemming  from  previous  difficulty  in  accounting 
quantitatively  for  any  significant  temperature  rise 
on  a  continuum  mechanics  model  description  of 
shear  banding  behavior.  The  discrete  dislocation 
pile-up  versus  continuum  mechanics  comparison  of 


such  behaviors  is  accounted  for  by  neglect, 
normally,  of  any  crack-like  stress-concentrating 
behavior  in  the  continuum  mechanics  modeled 
shear  band.  The  consideration  enters  into  the 
comparison  of  such  continuum  and  discrete  (or 
smeared)  dislocation  model  computations  in 
two  ways:  (1)  the  continuum-type  average  strain 
within  a  shear  band  is  assumed  generally  to  be 
uniform;  and,  (2)  the  non-uniformly-spaced 
dislocations  in  a  pile-up  produce  a  locally 
magnified  stress  that  is  enhanced  at  the  pile-up 
tip  in  direct  proportion  to  the  number  of 
dislocations. 

In  Figure  2,  the  difference  in  elastic  constant 
parameters  involved  in  determining  the 
dislocation  pile-up  avalanche  temperatures  for 
RDX  and  PETN  is  not  so  great,  even  producing 
predicted  hot  spot  temperatures  lower  for  PETN 
than  RDX  at  the  same  hot  spot  size.  This  is 
opposite  to  PETN  being  known  generally  to  be 
more  impact  sensitive  than  RDX.  However, 
intersection  of  the  computed  avalanche 
temperature  curves  with  their  respective 
explosion  temperatures  shows  that  PETN  is 
predicted,  anyway,  to  be  more  sensitive  than 
RDX,  that  is,  PETN  undergoes  initiation  at  a 
lower  temperature  than  RDX  because  of  the 
lesser  stability  of  the  PETN  molecule  on  a 
thermal  explosion  basis! 

Consideration  in  Figure  2  of  the  fundamental 
stabilities  of  the  different  energetic  molecules  is 
not  unrelated  to  the  discrete  dislocation  model 
application  to  evaluation  of  shear  displacements 
associated  with  dislocation  movement  at  the 
inter-molecular-scale  in  the  crystal  lattice 
structures  (9).  The  case  for  the  orthorhombic 
Pbca  crystal  structure  of  RDX  has  been 
considered  in  detail.  For  example,  identified 
[100](021)  slip  in  RDX  has  been  connected  on  a 
model  crystal  lattice  basis  with  reported  nitroso- 
compound  formations  identified  in  both  drop- 
weight  impact  and  combustion  test  specimens. 
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Figure  3.  Comparison  of  model  hot  spot  temperatures  for  dislocation  pile-up  avalanches  and  pore 
collapses,  spanning  estimated  hot  spot  dimensions  appropriate  to  impact  and  shock  regimes. 


The  point  here,  relative  to  the  earlier  description  of 
Eyring’s  comment  (1),  is  that  dislocation  (shearing) 
movement  in  the  relatively  dense  crystal  lattices  of 
these  energetic  crystals  is  quite  capable  of  bringing 
out-cropping  appendages  of  adjacent  molecules  into 
critical  reaction  coordinate  distances  (5). 

A  further  crystal-lattice-scale  consideration 
involves  application  of  the  discrete  dislocation 
model  to  analysis  of  proposed  nano-scale  events 
transpiring  in  shock  initiation  (10).  Here,  the  issue 
raised  is  complication  for  shocks  of  sufficient 
amplitude  that  are  considered  to  generate  a  nano¬ 
scale  dislocation  structure  at  the  shock  front.  Shock 
hardening  of  metallic  materials  is  well-known,  even 
to  the  degree  of  recent  description  being  reported  of 
three-dimensional  dislocation  dynamics  modeling 
of  unstable  dislocation  “channeling”  that  occurs  in 
neutron-irradiated  metals  (11).  The  same  type  of 


model  had  been  proposed  earlier  for  immediate 
occurrence  in  shocked  energetic  materials  in 
that  nano-scale  dislocation  loops  were  proposed 
to  form  at  the  propagating  shock  front,  then  to 
be  swept  through  by  subsequent  lagging 
movement  of  the  original  dislocation  network 
within  the  material  (12).  The  key  aspect  of  such 
mechanistic  consideration  is  that  the  very 
localized  subsequent  channeling  behavior  of  the 
originally  resident  dislocation  population 
through  the  nano-scale  network  promotes  shear 
banding  development,  and,  hence,  is  an 
additional  factor  to  be  considered  in  evaluating 
localization  of  the  dislocation  work  leading  to 
enhancement  of  any  temperature  rise. 
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Dislocation  shear  comparison  with  pressnred 
pore  collapse 

Recent  advances  in  dislocation  model 
descriptions,  as  have  been  described  in  the  modem 
literature,  and,  likewise,  recent  quantitative 
assessment  of  hot  spot  temperature  rises  to  be 
expected  from  the  pressure-induced  collapse  of 
pores  has  led  to  the  suggestion  of  making 
comparison  between  the  two  model  considerations 
(13).  Such  comparison  is  made  possible,  at  least 
on  a  preliminary  basis,  by  recent  quantitative 
assessment  made  by  Tarver  and  colleagues  for  the 
case  of  pore  collapse  (14). 

Figure  3  shows  a  first  attempt  at  comparison 
between  dislocation  pile-up  avalanche  and  pore 
collapse  models.  The  dislocation  basis  of  applying 
such  computation  to  FIMX  is  the  same  as  made  for 
RDX  and  PETN  in  Figure  2.  The  pile-up  avalanche 
dependence  has  been  superposed  onto  the  graphical 
results  developed  by  Tarver  and  colleagues  for 
initiation  of  HMX  and  TATB  (triaminotrinitro- 
benzene),  taken  to  be  caused  in  these  cases  by 
pressure-induced  pore  collapses.  The  critical 
conditions  for  the  upper  pore-size-equivalent  hot 
spot  temperatures  shown  in  the  Figure  for  HMX 
were  reportedly  determined  on  the  basis  of 
“multistep  Arrhenius  kinetic  chemical 
decomposition  models  derived  from  thermal 
explosion  and  kinetic  experiments”. 

The  indicated  attempt  at  comparing  the  two-type 
hot  spot  model  mechanisms  in  Figure  3  must  be 
assessed  with  the  thought  in  mind  that  “apples  and 
oranges”  are  possibly  being  compared.  Whereas 
the  hot  spot  diameter  for  the  pore  collapse 
mechanism  in  Figure  3  is  taken  as  the  pore  diameter 
for  the  collapse  mechanism,  the  same  abscissa 
diameter  in  the  Figure  is  taken  as  the  much  smaller 
estimated  hot  spot  diameter  produced  for  a  pile-up 
avalanche  initiated  within  a  1.0  mm  diameter 
crystal. 

Consider,  as  Tarver  has  mentioned,  that  the  pore 
sizes  may  be  comparable  to  the  crystal  sizes.  Then 
collapse  of  a  1.0  mm  pore  in  an  HMX  compact 
would  initiate  at  a  temperature  of  ~800K.  From  the 
same  Figure  results,  a  much  smaller  hot  spot 
diameter  would  lead  to  explosion  at  a  somewhat 
lower  temperature  marked  by  the  lower  filled  circle 
at  the  intersection  of  the  mechanical  work  and 
thermal  explosion  curves.  Such  pile-up  release  at 
the  greater  stress  intensity  associated  with  cracking 
would  be  capable  of  generating  the  highest 


temperature,  extending  to  the  upper  fdled  circle 
point  on  the  vertically-extended  pile-up 
avalanche  model  curve.  Thus,  the  preliminary 
indication  is  that  dislocation  pile-up  avalanches 
produced  by  a  shear  component  of  the  applied 
stress  tensor  may  be  capable  of  generating  hot 
spot  initiations  quite  comparable  to,  and  perhaps 
even  higher  than  ,  those  associated  with  similar 
sized  pore  collapses. 

Finally,  the  reminding  comment  should  be 
stated  that  such  comparison  as  has  been 
described  suffers  additionally  from  the  pile-up 
collapse  mechanism  being  evaluated  until  now 
by  overestimating  the  degree  of  hot  spot 
heating.  Recently,  more  realistic  model 
calculations  are  being  pursued  for  the  pile-up 
release  mechanism  (15),  as  monitored  by  the 
same  type  of  computationally  intensive  3- 
dimensional  parametric  dislocation  mechanics 
approach  mentioned  for  reference  1 1 . 
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Abstract 

Dislocations  in  molecularly  bonded  energetic  crystal  lattices  are  associated  with  an  interesting  contrast  in 
mechanical  properties:  the  crystals  are  elastically  soft,  plastically  hard,  and  brittle.  The  combination  of 
properties  is  shown  here  to  relate  to  the  important  consideration  of  mechanically  enhanced  chemical 
decomposition  possibly  leading  to  explosive  detonation.  A  variety  of  dislocation  Burgers  vectors  and  line 
vectors  occur  in  individual  crystals  because  of  the  characteristically  weak  intermolecular  bonding  among 
oddly-shaped  molecules  that  are  packed  into  otherwise  low-symmetry  crystal  structures.  Large  Burgers 
vector  dimensions  occur  as  a  result  of  the  weak  molecular  bonding  and  outweigh  the  low  magnitude  of 
elastic  stiffnesses  to  give  relatively  large  dislocation  self-energies,  consonant  with  low  dislocation  densities 
being  observed  in  such  crystals. 

Despite  the  weak  bonding,  dislocation  movement  is  reasoned  to  be  intrinsically  hindered  in  the  lattice 
structures  because  of  cross-blockages  among  the  intertwined,  oddly-shaped  molecules.  An  important 
consequence  is  that  the  crystals  exhibit  limited  operation  of  restricted  slip  or  twinning  systems.  Models  of 
the  molecular  displacements  occurring  for  identified  slip  systems  show  molecular  interactions  in  agreement 
with  detected  chemical  decompositions.  Alternatively,  energetic  crystals  have  relatively  low  (anisotropic) 
values  of  surface  energies  that  lead  to  low  values  of  theoretical  cleavage  stresses.  Thus,  energetic  crystals 
are  elastically  compliant,  relatively  hard,  and  brittle. 

The  relatively  high  crystal  yield  stresses  and  low  cleavage  stresses  provide  only  a  narrow  range  in  stress 
level  for  dynamic  dislocation  behavior.  In  this  range,  localized  hot  spot  heating  from  the  dislocation  work, 
though  involving  only  small  numbers  of  dislocations  in  pile-up  avalanches,  is  significantly  enhanced 
because  of  characteristically  low  material  thermal  conductivities.  In  shock  wave  studies,  hindrance  to 
dislocation  movement  at,  and  behind,  the  shock  wave  front  is  an  important  consideration,  as  is  the 
generation  of  nanoscale  dislocation  clusters  at  large  shock  pressures  and  the  role  of  anharmonic  lattice 
strains  at  the  dislocation  cores  in  allowing  either  "in  situ"  hot  spot  heating  or  direct  electronic  excitations. 

1 .  Crystal  dislocation  characteristics 

The  combination  of  large  Burgers  vectors  and  low  elastic  stiffnesses  for  energetic  crystal  lattices  leads  to 
large  dislocation  self-energies  and  consequently  low  dislocation  densities  in  crystals  that  are  mostly  grown 
from  solution.  The  low  symmetry  of  the  crystal  lattice  structures  leads  also  to  a  variety  of  dislocation 
Burgers  vectors  and  line  vectors;  however,  the  intertwined  packing  of  oddly-shaped  molecules  produces 
strong  intrinsic  resistance  to  dislocation  movement  that  allows  only  few  slip  or  twinning  systems  to 
operate.  The  crystal  surface  energies  are  low  and,  with  restricted  plastic  flow,  cleavage  is  relatively  easy. 

1.1.  Large  Burgers  vectors 

The  large  Burgers  vector  dimensions  in  molecularly  bonded  lattice  structures,  and  in  related  crystal 
structures  characterized  by  large  interatomic  or  intermolecular  lattice  parameters,  provided  an  opportunity 
in  the  early  days  of  dislocation  studies  of  optically  detecting  dislocation  end  points  at  intermolecular-scale 
ledge  heights  on  otherwise  flat  crystal  surfaces  exposed  during  growth  from  the  vapor  [1].  Very  fine 
striated  ledge  structures  had  been  observed  with  the  optical  microscope  on  the  surfaces  of  vapor-grown 
crystals  subjected  to  supersaturations  less  than  the  theoretical  value  required  for  growth  on  a  perfectly  flat 
surface.  Cabrera  and  Burton  showed  that  such  ledges  had  to  be  a  non-equilibrium  phenomenon  [2].  Frank 
[3]  gave  a  model  explanation  for  the  generation  of  such  ledges  spiraling  outward  from  emergence  sites  of 
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screw  dislocations  and  thus  explained  the  growth  of  crystals  at  low  vapor  supersaturations  through  the 
continuous  availability  of  such  ledges  for  attachment  of  atoms  or  molecules. 

Griffin  [4]  first  reported  from  the  Royal  Holloway  College  (RHC)  a  screw  dislocation-type  spiral  marking 
on  the  surface  of  a  beryl  crystal.  Verma  [1]  followed  shortly  thereafter,  also  while  at  the  RHC,  with  optical 
phase  contrast  microscopy  observations  of  spiral  ledges  terminating  at  presumed  dislocation  positions 
threading  otherwise  molecularly  flat  surfaces  of  silicon  carbide  crystals  brought  from  India.  With 
supervision  from  S.  Tolansky,  the  optical  microscope  technique  was  pushed  to  provide  quantitative 
measurement  using  multiple  beam  interferometry  of  a  1.5  nm  step  height  along  the  [0001]  growth  surface 
normal,  equivalent  to  a  c-axis  unit  cell  length  for  a  dislocation  Burgers  vector  in  a  6H  (hexagonal)  crystal 
[5].  Figure  1  shows  an  example. 

It  is  interesting  to  note  the  important  consequence  of  there  being  relatively  few  dislocations  with  such 
[0001]  Burgers  vectors  in  these  silicon  carbide  or  related  hexagonal  crystals  because  the  optical  microscope 
observations  and  measurements  were  facilitated  by  the  low  areal  density  of  dislocations.  Such  growth 
spirals  were  observed  on  silicon  carbide  crystals  also  by  Amelinckx  [6],  and  on  the  (0001)  surfaces  of  n- 
paraffin  crystals  by  Dawson  and  Vand  [7],  and  on  the  surface  of  an  n-nonatriacontane  crystal  by  Anderson 
and  Dawson  [8].  Not  long  after  report  of  these  optical  microscope  observations,  Menter  [9]  obtained 
electron  microscope  pictures  of  dislocations  in  copper  and  platinum  phthalocyanine  crystals  by  imaging  the 
Burgers  vector  displacements  along  1.2  nm  separations  of  (201)  lattice  planes. 

The  preceding  observations  were  made  before  dislocations  were  directly  observed  in  the  transmission 
electron  microscope  (TEM)  via  strain  field  imaging  of  the  defects  in  thin  metal  foils  [10,11].  The  TEM 
technique  allowed  high  magnification  observations  of  dislocations  to  be  made  through  strain-induced 
changes  in  the  local  diffracted  intensity,  say,  over  regions  extending  radially  tens  of  nms  from  the 
dislocation  cores,  even  with  the  much  smaller  Burgers  vectors  of  dislocations  in  metals  and  alloys.  Whereas 
a  low  dislocation  density  is  generally  required  to  make  successful  dislocation  step  height  measurements  on 
crystal  growth  surfaces  with  an  optical  microscope,  much  higher  dislocation  densities  are  conveniently 
studied  with  the  higher  magnifications  employed  in  the  TEM.  Unfortunately,  electron  beam  damage  often 
interferes  with  molecular  crystal  defect  observations  in  the  TEM. 

The  same  low  dislocation  density  requirement  of  optical  step  height  measurements  applies  for  the 
observation  of  dislocations  by  the  two  techniques  of  chemical  etch  pitting  or  x-ray  diffraction  topography. 
X-ray  topography  involves  optical  microscope  enlargement  of  strain-induced  contrast  in  the  diffracted  x- 
ray  intensity,  thus  occurring  at  dislocations  for  the  same  reason  as  in  the  TEM  technique,  but  here  in  Lane 
or  Bragg  diffraction  images  recorded  on  fine  grained  nuclear  emulsions.  Reference  [12]  gives  a  description 
of  dislocation  observations  by  optical,  x-ray,  electron  and  field  ion  microscopy.  The  chemical  etch  pitting 
and  x-ray  topography  techniques  have  been  applied  very  successfully  to  observing  dislocations  in  energetic 
crystals.  Connick  and  May  [13]  first  reported  etch  pitting  observations  of  dislocations  in 
cyclotrimethylenetrinitramine  (RDX)  crystals  and  Lang  [14]  first  reported  that  dislocations  had  been 
observed  in  RDX  crystals  by  transmission  x-ray  topography.  RDX  has  an  orthorhombic  crystal  structure 
with  space  group  Pbca  [15]. 

1.2.  Chemical  effects 

Nabarro  [16]  points  out  that  Faraday  [17]  reported  the  observation  of  efflorescence  occurring  at  mechanical 
scratches  put  into,  for  example,  surfaces  of  otherwise  perfectly  stable  crystal  surfaces  of  sodium  sulfate 
decahydrate.  Reference  to  the  original  article  shows  that  the  same  observation  was  made  on  other  related 
decahydrate  crystals.  Faraday  described  the  effects  "as  a  curious  illustration  of  the  influence  of  mechanical 
forces  over  chemical  affinity".  Nabarro  suggested,  because  efflorescence  involves  only  a  change  in  part  of 
the  molecular  unit  cell  stmcture,  that  imperfect  dislocations  might  be  involved  in  the  scratching 
deformation  as  compared  with  perfect  dislocations  participating  in  the  crystal  growth  process,  and  thereby 
being  relatively  inert.  He  also  pointed  out  that  etch  pits  and  dehydration  points  were  reported  to  be 
uncorrelated  in  gypsum. 
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The  segregation  of  solute  atoms  or  moleeules  at  the  eores  of  disloeations  in  erystals  of  any  type  is  a  natural 
eonsideration  beeause  of  the  severe  state  of  strain  at  the  disloeation  eenters,  and  oeeurs  frequently  enough 
to  have  led  to  eoining  of  the  term  "disloeation  deeoration".  This  influenee  is  well  known  in  atomie  diffusion 
studies  and,  also,  in  deformation  studies  of  diseontinuous  yielding  and  follow-on  plastie  flow.  For  the  ease 
of  large  Burgers  veetor  disloeations,  however,  the  strain  state  is  relatively  so  severe  as  to  allow  disloeation 
eores  in  vapor  grown  erystals  to  be  hollow.  The  eondition  for  this  oeeurring  was  theoretieally  predieted  by 
Frank  [18].  Also,  for  erystals  grown  from  solution,  as  is  the  main  method  of  produetion  of  energetie 
erystals,  the  disloeation  eores  may  eontain  a  residual  amount  of  solvent.  Thomas  [19],  while  eommenting 
also  on  Faraday's  observations  with  deeahydrate  erystals,  pointed  out  that  suffieient  water  oeeurred  in  the 
ehannels  of  disloeation  eores  in  suerose  erystals  to  appreeiably  affeet  their  eleetrieal  eonduetivities. 

Figure  2  shows  a  modem  x-ray  seetion  topograph  that  reveals  the  strain  field  of  a  serew  disloeation  in  a  6Ft 
silieon  earbide  erystal  plate  [20].  The  disloeation  eore  is  demareated  by  a  hollow  "mieropipe".  Seanning 
eleetron  mieroseope  (SEM)  images  were  obtained  of  the  emergent  holes  at  the  mieropipe  interseetions  with 
the  erystal  surfaee.  Different  mieropipe  diameter  measurements  were  shown  to  be  proportional  to  the 
square  of  the  x-ray  derived  disloeation  Burgers  veetors,  as  predieted  by  Frank  [18].  In  faet,  the  mieropipes 
were  only  observed  at  double  Burgers  veetor  or  larger  disloeation  eenters;  and,  the  largest  mieropipes  took 
on  an  hexagonal  eross-seetion.  Flirth  [21]  has  reeently  extended  the  Frank  model  eonsideration  to  inelude 
the  ease  of  a  mieropipe  being  represented  by  an  elliptieal  hole. 

1.3.  Self-energies 

Beyond  the  disloeation  eore  straeture  and  its  assoeiated  energy,  there  is  the  important  far  reaehing  linear 
elastie  strain  energy  eomponent  of  the  total  self-energy.  The  well-known  formula  for  the  elastie  strain 
energy  per  unit  length  of  the  disloeation  line  is  expressed  in  terms  of  the  produet  of  a  strain  energy 
eoeffieient,  involving  the  material  shear  modulus  and  square  of  the  Burgers  veetor,  and  the  logarithmie 
ratio  of  outer  and  inner  eut-off  radii.  The  inner  eut-off  radius  is  determined  by  the  disloeation  eore 
property;  however,  the  logarithmie  term  is  generally  thought  to  be  a  minor  variable  in  determining  the 
disloeation  energy. 

Table  1,  from  Armstrong  and  Elban  [22],  gives  in  eolumn  2  eomputed  values  of  the  strain  energy 
eoeffieient,  multiplied  by  unit  intermoleeular  length  along  the  disloeation  line,  for  a  [100](040)  edge 
disloeation  in  RDX  eompared  to  [1 10](1 10)  edge  disloeations  in  LiF  and  MgO.  The  elastie  strain  energy  is 
seen  to  be  greater  for  RDX.  In  eolumn  3,  the  elastie  strain  energy  eoeffieient  is  divided  by  the  heat  of 
formation  per  moleeule  for  eaeh  of  the  materials  to  show  that  the  disloeation  strain  energy  is  greater  for 
RDX.  Thus,  the  lower  shear  modulus  of  RDX  is  more  than  eompensated  by  the  larger  Burgers  veetor 
influenee  and,  also,  the  larger  referenee  length  along  the  disloeation  line.  Also  in  Table  1,  eomputed  values 
are  given  for  the  eonventional  Peierls-Nabarro  (P-N)  stress  that  depends  linearly  on  the  shear  modulus  and 
exponentially  on  the  ratio  of  interplanar  spaeing  and  Burgers  veetor  [23].  Although  RDX  has  a  lower 
eomputed  P-N  stress,  it  is  shown  to  be  a  larger  fraetion  of  the  shear  modulus.  Of  eourse,  sueh  estimate  is 
vitiated  by  the  eross-bloekage  between  adjaeent  moleeules  eonsideration  and,  so,  an  even  higher  disloeation 
movement  stress  should  be  expeeted. 

2.  Crystal  morphologies  and  perfeetions 

Energetie  erystals  form  a  subgroup  of  the  broader  elass  of  moleeular  organie  erystals  that  are  mostly  grown 
in  solution  from  supersaturated  liquids.  Gem-like  erystals  with  flat,  nearly  speeular,  external  habit  planes 
are  produeed  in  the  growth  proeess.  Disloeation  bundles,  as  oeeur  even  in  the  natural  produetion  of 
diamonds,  pieree  the  individual  planar  faeets  and  faeilitate  growth  by  produeing  moleeular  seale  step 
heights  in  aeeordanee  with  the  Frank  [3]  model  of  erystal  growth.  Internally,  the  erystals  are  built  up  of 
striated  growth  bands,  parallel  to  the  planar  faeets,  that  meet  at  impurity-assoeiated  growth  seetor 
boundaries. 

Conniek  and  May  [13],  in  a  study  of  RDX  erystal  growth,  provided  an  example  of  the  variety  of  external 
erystal  morphologies  that  ean  be  produeed  from  different  supersaturated  solvents.  Thus,  tabular  erystals 
with  dominant  (001)  habit  planes  were  produeed  by  growth  from  eyelohexanone  solution  while  prismatie 
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crystals  with  prominent  <001>  axes  were  grown  from  solution  in  acetone  or  dimethyl  formamide. 
Armstrong  and  Elban  [24]  reported  crystals  having  both  the  morphology  described  by  Connick  and  May 
and  a  new  crystal  morphology  with  prominently  exposed  {210}  surfaces  when  grown  by  slow  evaporation 
of  reagent-grade  acetone  with  Holston  RDX  Class  D  (specifying  1-2  mm  size  crystals)  put  into  solution. 
The  relatively  large  crystals  were  several  mm  in  size.  Klapper  [25]  has  produced  beautiful  color  pictures  of 
the  natural  appearances  of  such  inorganic  and  organic  crystals  as  cesium-alum,  benzil,  salol,  and 
benzophenone. 

Figure  3  shows,  from  Klapper  [26],  a  schematic  vapor-  or  solution-growth  representation  of  the  relation 
between  external  growth  planes  and  an  internal  microstructure,  first,  of  growth  sectors  with  growth  banding 
parallel  to  the  external  crystal  facets,  secondly,  of  growth  sector  boundaries  defining  the  intersection  of 
growth  sectors,  and,  lastly,  separate  “sprays”  of  straight  dislocation  bundles  spreading  through  individual 
crystal  growth  sectors.  As  indicated,  the  dislocation  bundles  are  oriented  roughly  perpendicular  to  the 
external  crystal  facets.  Figure  4  shows  a  sectioned  (001)  surface  of  an  RDX  crystal  that  has  been  etched  to 
reveal  its  internal  growth  sector  pattern  [27].  The  planar  (Miller-index  type)  surfaces  of  the  growth  bands 
are  identified  in  the  Figure  from  the  trace  intersections  based  on  a  stereographic  projection  method  of 
analysis  [28].  Such  internal  growth  structures  are  particularly  well-documented  by  x-ray  topography 
studies  of  both  solution-grown  and  vapor-grown  crystals  [26,29],  including  silicon  carbide  and  sapphire 
crystals. 

The  straight  line  dislocation  bundles,  that  are  very  often  impurity  decorated,  are  those  generally  associated 
with  the  type  of  ledge  structure  described  by  Frank  [3]  to  participate  fundamentally  in  the  crystal  growth 
process.  These  dislocations  have  been  shown  to  either  propagate  from  a  crystal  seed,  to  nucleate  at  the  seed 
interface,  particularly  at  inclusion  particles,  or  to  arise  from  perturbations  in  the  crystal  growth  conditions. 
In  the  latter  case,  the  observation  of  incorporated  liquid  inclusions  has  been  attributed  to  local  variations  of 
supersaturation  and  solvent  flow  at  the  liquid-crystal  interface  [26].  Gross  [30]  reported  the  observation  of 
solvent  inclusions  in  RDX  crystals.  Klapper  accounts  for  the  straight  line  character  of  the  dislocation 
bundles  on  a  minimum  energy  basis,  thus  favoring  the  identification  of  screw  dislocation  growth,  as 
described  by  Frank  [3].  Flowever,  bundled  dislocations  with  edge  and  mixed  screw  plus  edge  characters 
have  been  identified  also  to  participate  in  the  crystal  growth  process. 

Both  the  dislocation  core  energy  and  the  anisotropic  character  of  elastic  constants  that  determine  the 
coefficient  for  the  elastic  self-energy  [31],  as  mentioned  for  Table  1,  enter  into  the  determination  of  the 
dislocation  growth  line  orientations.  The  striated  growth  bands  themselves  are  caused  by  local  variations  in 
the  impurity  content  grown  into  the  crystals  through  the  ledge  structures.  The  growth  sector  boundaries 
appear  to  mark  small  localized  impurity-dependent  mismatches  in  lattice  parameter  or  slight  planar 
misorientations.  Van  Enckevort  and  Klapper  [32]  report  a  modem  example  of  matching  x-ray  topography 
measurements  with  highly  sensitive  differential  interference  contrast  optical  microscopy  showing  0.9  and 
1.8  nm  (unit  cell  and  half-unit  cell,  respectively)  steps  on  the  growth  faces  of  nickel  sulfate  hexahydrate 
crystals,  relating  to  the  observations  mentioned  above  from  Faraday  [17]  and  Verma  [1].  In  this  case,  pure 
screw  dislocation  out-crops  were  responsible  for  the  growth  features. 

McDermott  and  Phakey  [33]  showed  by  x-ray  topography  and  etch  pitting  a  one-to-one  match  of  a  small 
dislocation  bundle  with  etch  pits  of  the  individual  dislocations  on  the  crystal  growth  surface  of  an  RDX 
crystal,  as  shown  in  Figure  5.  An  opposite  case  of  out-cropping  pill  boxes  at  dislocation  line  end  points  was 
identified  also  on  the  surfaces  of  the  same  batch  of  RDX  crystals  studied  in  cross-section  by  van  der  Steen 
and  Duvalois,  as  shown  with  identified  plane  indices  in  Figure  6  [27,28].  Rather  detailed  x-ray  topography 
results  were  reported  for  dislocation  Burgers  vectors  and  line  vectors  in  laboratory-grown  RDX  crystals  by 
Flalfpenny  et  al.  [34].  A  large  variation  in  both  dislocation  Burgers  vectors  and  line  vectors  was 
determined.  The  variety  of  dislocation  line  shapes  in  carefully  grown  crystals  of  cm  scale  dimensions, 
including  the  observation  of  dislocation  helices,  was  interpreted  to  indicate  that  significant  post-growth 
dislocation  movement  had  occurred.  A  combination  of  x-ray  topography  and  etch  pitting  results  were 
reported  also  to  lead  to  expectation  of  a  variety  of  dislocation  types  in  energetic  pentaerythritol  tetanifrate 
(PETN)  ciystals  [35]. 

3.  Mechanical  properties 
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Despite  the  large  variety  of  dislocation  Burgers  vectors  and  line  vectors  in  energetic  crystals,  the  crystals 
are  hard  and  easily  cracked.  The  behavior  has  been  traced  to  the  difficulty  of  mechanically-forced 
dislocation  movement.  Such  movement,  whose  dynamic  characteristics  are  narrowly  limited  between 
relatively  high  yield  and  low  cleavage  stresses,  is  capable  of  conveying  outcropping  appendages  of  adjacent 
molecules  into  critical  reaction  coordinate  distances.  Chemical  reaction  is  enhanced  by  the  localized  "hot 
spot"  heating  produced  by  the  dislocation  work. 

3.1.  Early  mechanical  testing  of  anthracene/naphthalene. 

Kochendorfer  [36]  did  early  work  on  the  deformation  of  naphthalene  as  an  example  of  determining  slip 
systems  for  a  molecular  crystal,  in  this  case,  having  a  monoclinic  crystal  structure  and  nm  scale  lattice 
parameters  (a  =  0.86  nm,  b  =  0.60  nm  and  c  =  1.12  nm).  Robinson  and  Scott  [37]  followed  later  with  tensile 
stress-strain  measurements  of  the  critical  shear  stress  for  structurally  similar  anthracene  crystals  as  a 
function  of  temperature  and  strain  rate.  Although  the  ends  of  the  melt-grown  crystal  rods  were  cleaved 
along  (001)  surfaces  to  provide  shorter  tensile  specimens,  the  thus  tested  crystals  were  quite  ductile  with 
single  slip  occurring  on  the  (001)  [010]  system  or,  with  greater  strain  hardening,  the  (001)[110]  system. 
Robinson  and  Scott  pointed  out  that  the  same  slip  systems  were  known  then  to  occur  in  pure  naphthalene 
crystals. 

Despite  lack  of  knowledge  of  the  thirteen  elastic  constants  for  the  anthracene  crystal  structure,  the  P-N 
stress  was  reasoned  by  Robinson  and  Scott  to  be  an  important  factor  in  determining  the  crystal  strength 
properties.  Of  greater  importance  in  this  study  was  the  schematic  demonstration  on  a  lattice  model  basis,  as 
indicated  by  identified  plane  projections  and  directions  in  Figure  7a,b,  that  slip  should  occur  relatively 
easily  across  (001)  along  [010]  but  should  occur  with  greater  difficulty  along  [110].  Although  arranged  in 
a  layered  herring-bone  lined-up  array  of  uniaxial  molecules,  the  crystal  model  nevertheless  shows 
reasonably  unobstructed  pathways  for  the  proposed  [010]  shear  displacements  between  the  rows  of 
molecules. 

Perhaps  more  closely  related  to  results  to  be  described  for  an  energetic  crystal,  di  Persio  and  Escaig  [38] 
reported  "in-situ"  x-ray  topographic  observations  made  during  mechanical  testing  of  a  reasonably  simple 
body-centered  cubic  structure  obtained  for  hexamethylene  tetramine  crystals.  Deformation  in  constant 
strain  rate  compression  tests  and  in  tensile  creep  tests  was  monitored.  The  crystals  were  free  of  growth 
bands  and  contained  very  low  dislocation  densities.  Cracking  occurred  often  in  the  compression  tests,  even 
at  the  lowest  machine  testing  speed.  From  the  compression  tests  and  tensile  creep  observations,  it  was 
concluded  that  the  dislocations  must  have  a  large  P-N  stress,  also  reasoned  to  occur  because  of  molecular 
misfit  difficulties  associated  with  displacements  in  the  dislocation  cores.  Difficulties  proposed  for  adjacent 
intermolecular  nitrogen  -  hydrogen  interactions  during  dislocation  movement  in  a  <1 1 1>  slip  direction  were 
briefly  discussed. 

3.2.  Micro-indentation  hardness  testing. 

The  brittleness  of  most  energetic  crystals  and,  also,  interest  in  testing  small  material  samples  for  personal 
safety  led  to  micro-indentation  hardness  testing  providing  an  excellent  means  of  obtaining  information  on 
crystal  plasticity  and  cracking  mechanisms  in  these  materials.  Elban  and  Armstrong  [24]  first  reported  etch 
pitting  evidence  of  highly  localized  plastic  flow  surrounding  Knoop  micro-indentation  hardness 
impressions  in  RDX  crystals.  Figure  8  provides  an  example  of  etch  pits  surrounding  various  Knoop  and 
diamond  pyramid  indentations  put  into  the  (210)  surface  of  an  RDX  crystal  [24,39,40].  The  areas  covered 
by  the  hardness-produced  dislocation  pits  in  RDX  were  approximately  eleven  times  greater  than  the 
residual  impression  sizes,  as  compared  with  an  area  about  ninety  times  larger  for  the  x-ray  determined 
plastic  strain  field  of  an  indentation  in  LiF.  Previous  measurements  on  MgO  crystals  showed  etch  pit  and  x- 
ray  contrast  zones  were  roughly  comparable  [41].  Flalfpenny  et  al.  [35]  also  reported  etch  pitting  evidence 
of  restricted  plastic  flow  at  indentations  in  RDX  crystals  and,  in  PETN  crystals,  somewhat  larger  areas  of 
very  directionally  distributed  dislocation  pits  were  observed  at  indentations.  In  a  subsequent  study, 
Gallagher  et  al.  [42]  investigated  the  orientation  dependence  of  Knoop  microhardness  measurements  in 
both  RDX  and  PETN  crystals  and  correlated  the  results  with  a  model  of  the  hardness  anisotropy  being 
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determined  by  the  orientation  dependenee  assoeiated  with  aehieving  the  eritieal  resolved  shear  stress  for 
disloeation  slip. 

Well-defined  erystallographie  eraeking  was  observed  at  the  miero-indentation  sites  in  eaeh  of  the  preeeding 
studies  [24,35,40].  Earlier,  Hagan  and  Chaudhri  [43]  measured  the  eraek  size  dependenee  of  diametral 
eraeking  at  diamond  pyramid  miero-indentations  in  RDX  and  PETN  erystals  and,  in  line  with  a  model  of 
the  predieted  eraek  size  dependenee  on  load  [44],  determined  fraeture  surfaee  energies  for  sueh  eraeking. 
These  experiments  led  Elban  [45]  to  make  a  eomparison  of  the  meehanieally  determined  surfaee  energies 
and  those  determined  separately  from  eontaet  angle  measurements  in  wettability  studies.  The  fraeture 
surfaee  energies  were  found  to  be  larger  by  less  than  a  faetor  of  two  for  RDX  and  less  than  a  faetor  of  three 
for  PETN.  Possible  explanations  for  the  differenees  in  the  energies  were  given  of  plastie  work  being 
assoeiated  with  the  eraeking  proeess  or  signifieant  anisotropy  in  energy  being  assoeiated  with  different 
surfaees.  Hagan  and  Chaudhri  eommented  that  substantial  plastie  flow  had  oeeurred  in  assoeiation  with  the 
eraeking  in  PETN.  Also,  they  found  a  higher  fraeture  surfaee  energy  in  RDX  for  (100)  eraeking  than  (001) 
eraeking,  in  agreement  with  the  qualitative  observation  by  Conniek  and  May  [13]  that  (001)  eraeking  was 
easiest. 

There  are  two  disloeation  flow  eonsiderations  that  bear  on  affeeting  meehanieally-indueed  fraeture  surfaee 
energy  measurements:  eraeks  are  formed  or  small  eraeks  made  larger  beeause  of  disloeation  flow;  and, 
disloeation  generation  and  movement  ean  oeeur  at  eraek  tips  to  blunt  or  otherwise  restriet  eraek  growth. 
The  role  of  disloeation  flow  in  forming  a  eraek  is  well-illustrated  by  the  ease  of  {110}<110>  slip  on 
adjaeent  slip  planes  in  MgO  erystals  eombining  to  form  eleavage  eraeks  on  the  relatively  unfavorable 
{110}  eleavage  planes  [46].  Figure  9  shows  a  eomparison  of  residual  diamond  pyramid  diagonal  lengths 
and  tip-to-tip  diagonal  eraek  extensions  measured  for  indentations  on  (001)  MgO  and  (210)  RDX  erystal 
surfaees  [22,40,47,48].  The  em-size  RDX  erystal  on  whieh  the  major  (triangular  point)  measurements  had 
been  made  was  supplied  by  H.H.  Cady,  Los  Alamos  National  Laboratory.  Note  the  reasonable  agreement 
with  Hagen  and  Chaudhri’s  reported  eraek  size  measurements  [43]  on  RDX,  although  the  lowest  load 
datum  measured  by  them  is  shown  in  the  Figure  to  deviate  from  the  theoretieal  indentation  exponential 
(3/2)  dependenee  predieted  on  an  indentation  fraeture  meehanies  basis  [44].  This  low  indentation  load 
datum  plus  the  low  load  measurements  shown  for  MgO  are  in  the  regime  where  visible  eraeking  first 
oeeurs. 

In  Figure  9,  the  foree  dependeneies  of  the  indentation  diagonal  lengths  and  diagonal  tip-to-tip  eraek  sizes  of 
{110}  eraeks  in  MgO  and  (001)  eraeks  in  RDX  present  an  interesting  eomparison.  The  (left-side)  line  of 
slope  2  in  the  Figure  applies  for  a  eonstant  hardness.  For  MgO,  the  dependenee  of  the  indentation  foree  on 
indentation  size  follows  the  same  indentation-fraeture-meehanies-based  (3/2)  dependenee  as  shown  for  the 
eraek  size  measurements.  The  eraeking-eontrolled  deerease  in  hardness  at  larger  indentation  sizes  was 
attributed  to  the  positive  influenee  of  eraeking  on  dissipating  the  pent-up  disloeation  pile-up  strain  energies. 
X-ray  topography  evidenee  was  obtained  on  MgO  showing  redueed  strain  fields  assoeiated  with  sueh 
eraeking  in  larger  ball  indentation  tests  [47].  By  eontrast,  RDX  shows  in  Figure  9  an  apparently  negligible 
influenee  of  eraeking  on  the  material  hardness,  that  is  seen  to  be  essentially  eonstant  at  all  eraek  sizes.  The 
eombined  indentation  and  eraek  size  measurements  for  RDX  are  eontrolled  by  the  disloeation  flow, 
eonsistent  with  the  proposed  intermoleeular-meehanism-based  diffieulty  of  disloeation  movement. 

Another  hardness-based  method  for  eomparison  of  plastie  flow  or  eraeking  behavior  that  is  informative  of 
the  position  of  RDX  and  related  energetie  materials  among  a  wider  range  of  materials  is  shown  in  Figure 
10,  as  developed  by  Hammond  and  Armstrong  [49].  Here,  a  hardness  stress,  defined  as  the  mean  pressure 
on  a  ball  indentor  of  diameter,  D,  is  plotted  versus  a  measure  of  the  hardness  strain,  defined  by  the  ratio  of 
the  eireular  eontaet  diameter,  d,  divided  by  D.  For  elastie  loading  behavior,  a  linear  hardness  stress 
dependenee  on  (d/D)  ean  be  eomputed  from  the  elastie  eonstants  of  the  ball  and  test  material  in  aeeordanee 
with  the  theory  of  Hertz  [50].  In  Figure  10,  all  of  the  linear  dependeneies  are  eomputed  for  a  steel  ball 
indenting  the  various  materials  that  are  identified:  AP  is  ammonium  perehlorate;  and,  anthraeene,  as 
diseussed  above,  is  shown  on  the  left  side  of  the  Figure  to  have  the  lowest  elastie  loading  eurve,  that  is  to 
be  eompared,  on  the  right-side  Figure  loeation,  with  indieation  of  a  low  value  for  the  eonventional  diamond 
pyramid,  or  Viekers  hardness  number,  VHN.  The  elastie  loading  lines  are  raised  in  order  of  inereasing 
material  elastie  stiffnesses.  Also,  at  the  end  points  of  the  elastie  loading  lines  for  a  number  of  materials  are 
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the  theoretical  hardness  stresses  for  cracking,  ac,  as  computed  on  an  indentation  fracture  mechanics  basis 
for  different  ball  sizes  [51,52].  The  Hertzian  cracking  stress  is  seen  to  be  lower  at  larger  ball  sizes. 

Experimental  ball  indentation  test  results  are  shown  in  Figure  10  (as  open  square  points)  for  an  RDX  (210) 
crystal  surface  indented  with  a  steel  ball  of  1.59  mm  D  value  and  d  values  determined  from  optical 
measurements  of  the  residual  indentation  sizes.  As  noted,  several  smaller  indents  were  free  of  cracking. 
The  hardness  values  with  associated  cracks  are  lower  than  the  theoretical  Hertzian  stress  for  cracking, 
shown  at  the  terminus  of  the  elastic  loading  line,  because  of  the  role  of  dislocations  on  promoting  cracking. 
However,  the  plastically-cracked  material  measurements  are  not  very  far  below  the  theoretical  elastic 
cracking  stress,  and  so  Figure  10  gives  indication  of  only  a  relatively  limited  range  in  stress  level  for 
dynamic  dislocation  behavior  preceding  cracking  in  RDX,  say,  compared  to  the  much  larger  range  shown 
for  NaCl. 

Diamond  pyramid,  or  Vickers,  hardness  test  results  are  supposed  to  give  a  plastic  hardness  stress  equivalent 
to  that  obtained  in  a  ball  indentation  test  with  (d/D)  =  0.375  [53],  and  so,  on  the  right  side  of  Figure  10, 
representative  diamond  pyramid  hardness,  or  VHN,  values  are  shown  for  a  variety  of  materials.  Consider 
the  comparative  hardness  curves  shown  for  MgO,  RDX,  and  NaCl  crystals,  in  the  latter  case,  where  a 
continuous  hardness  stress-strain  determination  produced  excellent  agreement  with  the  measured  value  of 
the  conventional  hardness  [50,49].  MgO  is  at  the  top  of  the  vertically  distributed  conventional  hardness 
values  with  an  indicated  hardness  of  approximately  7  GPa  (1  GPa  =  98.1  kgf/mm^),  falling  below  a 
theoretical  Hertzian  stress  value,  for  the  same  hardness-determined  ball  size,  in  excess  of  20  GPa.  The 
lower  hardness  is  obtained  because  of  the  dislocation  promoted  cracking  that  has  been  described.  NaCl  has 
a  lower  elastic  stiffness  than  MgO,  of  course,  but  greater  elastic  stiffness  than  RDX.  However,  RDX  has  a 
greater  plastic  hardness  than  NaCl  and  lower  theoretical  cracking  stress.  The  comparison  establishes  a 
quantitative  basis  for  concluding  that  RDX  is  elastically  compliant,  plastically  hard,  and  relatively  brittle. 

3.3.  Drop-weight  impact  testing  for  crystal  sensitivity. 

Early  work  on  the  initiation  of  explosives  by  mechanical  loading,  particularly,  by  drop-weight  impact 
testing,  led  to  an  explanation  of  localized  hot  spot  heating  being  responsible  for  the  effect  [54,55].  Field  et 
al.  [56]  set  out  the  main  mechanisms  for  the  generation  of  hot  spots  as:  (1)  adiabatic  compression  of  gas 
pores;  (2)  viscous  heating  of  extruded  material  between  grains  or  at  external  surfaces;  (3)  frictional  heating 
between  contacting  surfaces;  and,  (4)  localized  adiabatic  shear  during  mechanical  failure.  Modem 
demonstration  of  such  hot  spots  was  produced  by  Coffey  and  Jacobs  in  drop-weight  impact  tests  of 
explosive  particles  placed  onto  heat  sensitive  film  [57]. 

Important  evidence  for  discontinuous  load  drop  behavior,  analogous  to  yield  point  behavior,  being 
associated  with  the  hot  spot  heating  and  chemical  decomposition  was  provided  by  Heavens  and  Field  [58] 
who  monitored  the  transmitted  pressure  wave  in  impact  tests  on  thin  layers  of  powder  particles.  The  results 
are  summarized  in  Table  2  in  columns  under  schematically-drawn  impact  stress-time  responses  shown  for 
an  assortment  of  energetic  and  inert  materials  that  were  tested.  Explosive  initiations  were  obtained  for  the 
“sharp  drop”  type  curves  characteristic  of  the  energetic  materials  shown  in  the  right  hand  column  of  Table 
2,  and  for  ammonium  perchlorate  (NH4CIO4,  AP)  in  the  second-from-right  hand  column.  Other  evidence 
for  the  need  of  the  deformation  to  be  concentrated  in  narrow  shear  bands  associated  with  the  sharp  load 
drop  behavior  and,  also,  with  adiabatic  heating,  was  obtained  in  a  companion  study  by  Winter  and  Field 
[59],  who  investigated  the  reverse  situation  of  small  particles  impacting  onto  larger  explosive  crystal 
surfaces.  Such  adiabatic  shear  bands  are  associated  with  discontinuous  load  drops. 

Armstrong  et  al.  [60]  pointed  out  that  achievement  of  appreciable  hot  spot  temperatures  at  the  slip  band 
level  required  the  non-uniformly  distributed  dislocations  to  be  blocked  by  strong  obstacles  and,  also,  that 
the  pile-up  energy  had  to  be  dissipated  by  an  avalanche -type  of  collapse  of  the  obstacles.  Theoretical 
description  of  the  dynamic  release  of  a  pile-up  had  been  described  on  a  dislocation  group  basis  by  Head  for 
the  case  where  there  is  a  linear  dislocation  velocity  dependence  on  the  shear  stress  [61].  Figure  11  gives 
an  example  of  sequential  steps  of  spreading  (to  the  left)  of  an  initial  pile-up,  viewed  on  a  continuous 
distribution  of  infinitesimal  dislocations  basis,  after  the  pile-up  block  was  removed  suddenly.  The  local 
dislocation  density,  compared  to  the  average  value  in  the  normalized  initial  slip  length,  is  shown  to  spread 
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rapidly  on  initial  release  of  the  pile-up.  At  one  hundredth  of  the  time  it  would  take  for  a  single  dislocation 
to  travel  the  length  of  the  pile-up  when  driven  by  the  external  stress,  and  with  a  linear  stress-velocity  law, 
about  one-fifth  of  the  released  dislocations  have  raced  ahead  to  a  distance  of  about  one  tenth  of  the  unit 
pile-up  length.  Thus,  the  consideration  obtains  that  an  initial  build-up  of  strain  energy  might  proceed 
isothermally  until,  at  cristical  microstmctural  stress  intensity,  breakthrough  occurs  and  the  stored  energy  is 
released  adiabatically  [60].  Other  related  model  descriptions  of  dynamic  dislocation  pile-up  properties 
have  been  given  by  Gerstle  and  Dvorak  [62]  and  Ockendon  and  Ockendon  [63].  Such  pile-ups  are  now 
known  to  be  an  integral  part  of  slip  band  propagation  against  planar  boundary  obstacles  in  all  materials 
including,  for  example,  ice  polycrystals  [64]. 

An  important  feature  of  significant  dislocation  pile-up  activity  coming  into  play  in  determining  the 
mechanical  properties  of  a  material  is  that  crystal  size  or,  within  a  polycrystal,  the  grain  size  influence 
becomes  very  recognizable  [65].  For  polycrystalline  metals,  the  onset  of  yield  point  behavior  for  plastic 
yielding  [66,67]  and  subsequent  plastic  flow  stress  [68]  both  follow  a  Flall-Petch  dependence  on  the 
reciprocal  square  root  of  the  polycrystal  grain  diameter.  An  increased  slope  for  the  grain  size  dependence, 
now  termed  the  microstmctural  stress  intensity  after  fracture  mechanics  connection,  is  associated  with 
greater  pile-up  activity  [65].  Somewhat  analogous  pile-up  formation,  and  sudden  release,  at  slip  band 
intersections  is  directly  observable  via  birefringent  transmitted  light  experiments  during  tensile  deformation 
of  MgO  crystals.  Also,  the  subject  of  crystal  size  influence  on  initiation  and  explosion  of  energetic  crystals 
has  been  of  scientific  interest  since  early  impact  deformation  results  were  reported  by  Bowden  and  Singh 
[69]. 

Application  of  the  pile-up  model  to  drop-weight  impact  test  results  obtained  on  RDX  crystals  led  to 
prediction  of  greater  drop-weight  heights  being  required  to  initiate  smaller  size  crystals  on  a  logarithmic 
basis  of  drop-weight  height  versus  reciprocal  square  root  of  crystal  size  [70].  Figure  12  shows 
experimental  measurements  demonstrating  the  effect  for  RDX  crystals.  A  strongest  obstacle  of  cleavage 
cracking  was  assumed  to  release  the  pile-up  avalanche  for  hot  spot  generation.  Bowden  et  al.  [71]  had 
proposed  that  the  speed  of  cracking  itself,  with  an  associated  plastic  zone  at  its  tip,  could  be  reponsible  for 
chemical  decomposition.  Fox  and  Soria-Ruiz  [72]  produced  evidence  for  the  consideration  from 
experiments  on  several  carbonate  and  azide  crystals.  In  the  latter  study,  an  interesting  association  of  crystal 
hardness  and  estimated  crack  tip  temperature  was  described. 

Connection  of  the  pile-up  avalanche  model  for  deformation-induced  crystal  initiation  was  made,  also,  with 
theoretical  work  of  Boddington  [73]  who  derived  a  relationship  between  the  hot  spot  size  and  temperature 
required  for  direct  initiation  of  thermal  explosion.  Reasonable  estimations  were  determined  from  the  pile- 
up  model  of  deformation-induced  temperature  rises  [74]  and  of  hot  spot  sizes  and  lifetimes  [70]. 
Comparison  of  calculated  RDX  and  PETN  crystal  deformation-induced  hot  spot  temperatures,  relative  to 
their  thermal  explosion  temperatures,  led  to  the  conclusion  of  PETN  being  more  sensitive  because  of  the 
lessor  stability  of  its  molecular  structure.  For  more  ductile  metallic  materials,  the  dislocation  pile-up 
avalanche  model  has  been  suggested  to  provide  the  fundamental  mechanism  for  adiabatic  shear  banding  at 
small  plastic  strains  [75].  Based  on  the  model.  Figure  13  shows  a  derived  graphical  basis  for  comparison  of 
the  shear  banding  susceptibilities  of  a  number  of  metallic,  ionic  and  energetic  materials  [76].  In  the  Figure, 
a  greater  slope  to  the  material  point  is  indicative  of  greater  shear  band  susceptibility. 

The  two  dependencies  shown  in  Figure  13  for  the  theoretical  avalanche-induced  temperature  rise  apply  for 
different  material  property  combinations,  among  which  the  microstructural  stress  intensity,  ks,  and  thermal 
conductivity,  K,  show  the  greatest  material  variations  [60].  For  metals  and  ionic  solids,  then,  lines  of 
greater  slope  to  the  different  (open  circle)  material  points  in  the  Figure  establish  a  greater  shear  band 
susceptibility,  that  reasonably  agree  with  the  (closed  circle)  upper  limiting  temperature  rise  predictions 
shown  using  the  right-hand  ordinate  scale  of  the  Figure.  In  the  Figure,  RDX,  PETN  and  AP  follow  the 
second  inequality  and  are  seen,  at  the  necessarily  reduced  scale  dimensions,  to  be  relatively  shear  band 
prone  even  compared,  for  example,  to  a-iron  whose  low  temperature  shear  banding  behavior  is  very  well- 
known. 

3.4.  Shock  loading. 


Ill 


Shock  wave  phenomena,  modeled  hydrodynamically  to  involve  discontinuities  in  pressure,  temperature, 
and  density,  across  an  idealized  shock  “front”,  are  associated  with  explosive  decompositions  of  energetic 
materials  and  their  formulations  [76].  A  one-dimensional  state  of  strain  obtains  in  the  shock  propagation 
direction.  For  the  real  case,  the  shock  front  is  taken  to  be  a  narrow  zone  of  nanometric-scale  dimensions 
experiencing  shear  stresses  comparable  to  the  theoretical  strength. 

Dislocation  model  aspects  of  shock  wave  loading  began  with  a  pioneering  description  by  Smith  of  the 
creation  of  dislocations  at  the  shock  wave  front  in  metals  [77].  Meyers  [78]  has  extended  the  model  in  a 
modem  account.  Figure  14  shows  another  model  description  of  distributed  nanometer-scale  dislocation 
loops  proposed  to  produce  a  residual  state  of  one-dimensional  strain  in  the  shock  propagation  direction 
across  a  presumed  non-planar  shock  front  [79].  An  analysis  has  been  given  of  the  dislocation  stracture 
[80].  The  shear  stress  state  across  the  compressional  front  is  shown  in  the  top  left  comer  of  the  Figure  near 
to  the  indicated  sectioned  crystal  edge.  Assumed  nanometer-scale  dislocation  shear  loops  are  distributed  so 
as  to  give,  at  the  curved  arrow  positions,  reacted  dislocations  with  Burgers  vectors  parallel  to  the  shock 
propagation  direction,  hence,  the  residual,  essentially  one-dimensional,  strain-state. 

The  presence  of  such  a  nanometer-scale  stmcture  should  promote  shear  banding,  also,  for  subsequent 
plastic  flow  from  the  originally  resident  dislocation  stmcture  whose  thermally-activated  movement  lags 
behind  the  shock  front  that  travels  at  a  velocity  elastically-determined  for  one-dimensional  strain  [81].  The 
situation  is  proposed  to  be  analogous  to  dislocation  channeling  that  occurs  in  post-deformation  of  neutron- 
irradiated  metals,  where  follow-on  dislocation  movement  is  made  easier  through  slip  band  channels  cleared 
by  first  dislocation  movements  sweeping-out  the  neutron-induced  dislocation  loops.  Otherwise, 
considerable  research  has  been  performed  by  Gupta  and  colleagues  [82,  83]  on  the  shock  wave  deformation 
of  LiF  crystals,  for  example,  as  described  in  a  later  paper  on  investigating  details  of  stress  relaxation  at 
shocked  crystal  surface  layers  previously  hardened  by  solute  diffusion  [84]. 

Special  challenges  occur  for  dislocation  mechanics  in  the  regime  of  shock  wave  loading  both  because  of 
the  nanometric  size  of  events  at  the  shock  front  and  its  speed  of  passage  that  requires  dislocation-connected 
reaction  times  from  picoseconds  to  nanoseconds  [85].  Results  in  the  low  shock  regime  are  presumed  to 
connect  with  normal  high  strain  rate  deformation  results,  for  example,  as  achieved  in  split  Flopkinson 
pressure  bar  tests  and  Taylor  cylinder  impact  tests  [81].  At  these  strain  rates,  deformation  twinning  is 
favored  over  dislocation  slip  for  a  number  of  metals  and  alloy  systems  [86].  Follow-on  thermally-activated 
slip  trails  behind  the  shock  front.  For  larger  amplitude  shocks,  creation  of  a  dislocation  nanostructure 
occurs  at  the  shock  front  and,  then  interaction  of  the  trailing  viscoplastic  flow  with  the  newly  created 
nanoscale  dislocation  clusters  is  an  important  consideration.  At  the  highest  shock  amplitudes,  "in  situ" 
initiation  is  envisioned  to  occur  within  the  strain  fields  of  dislocation  cores  [85,81].  At  the  most 
fundamental  level,  the  mechanism  of  vibrational  energy  (up  pumping)  enhancement  should  be  favored  at 
dislocation  cores  because  of  the  importance  of  anharmonic  coupling  to  the  process  [87-89]. 

An  early  experiment  on  the  shock  initiation  of  individual  cm-size  RDX  crystals,  measured  on  a 
microsecond  time  scale,  was  reported  by  Adams  and  colleagues  [90].  The  crystals  were  found  to  be  less 
sensitive  to  shock  initiation  than  powdered  material;  however,  once  shocked  the  crystals  were  significantly 
more  sensitive  to  a  second  shock.  The  result  was  explained  on  the  basis  that  the  first  shock  produced 
internal  damage  that  provided  sources  of  hot  spots.  Sandusky  et  al.  [91]  obtained  sub-microsecond  optical 
observations  of  luminous  slip  and  cracking  systems  trailing  behind  the  shock  front  of  AP  crystals  supplied 
by  T.L.  Boggs,  Naval  Air  Warfare  Center,  China  Lake,  CA.  Significant  light  generation  occurred  from  the 
shock  interaction  with  microhardness  impressions  put  into  the  crystal  surface  before  shocking.  X-ray 
photoelectron  spectroscopy  (XPS)  measurements  made  on  the  shocked  AP  crystals  demonstrated  that 
perchlorate  decomposition  occurred  at  the  shocked  indentation  sites  [92].  Fligher  shock  pressures  applied  to 
RDX  crystals  were  also  shown  to  produce  deformation-associated  luminosity  in  RDX  and 
cyclotetramethylenetetranitramine  (FIMX)  crystals  [93]. 

The  importance  of  crystal  orientation  on  the  shock  initiation  sensitivity  of  cm-size  body-centered  tetragonal 
PETN  crystals  was  established  in  a  series  of  papers  by  Dick  and  colleagues  [94-96].  The  experiments  built 
onto  very  complete  information  being  known  about  the  linear  and  volume  compression  of  shocked  PETN 
crystals  [97],  their  elastic -plastic  properties  [98]  and  elastic  constants  [99].  The  times  and  shock 
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propagation  distances  from  the  impacted  crystal  surface  for  “mn-to-detonation”  were  very  different  for 
different  crystal  orientations.  Match  with  the  dislocation  Burgers  vectors  and  slip  systems  led  to  correlation 
of  this  measure  of  crystal  sensitivity  for  <110>  directed  shocks  with  difficult  shearing  across  {100}  planes 
of  interleaved  molecules,  where  it  was  reasoned  significant  "steric  hindrance"  should  oppose  molecular 
displacements. 

For  <100>  directed  shocks  that  should  produce  shearing  across  the  primary  {110}  slip  planes,  having 
essentially  no  interleaved  molecules,  relatively  easy  slip  occurred  and  the  crystals  did  not  detonate. 
Intermolecular  distances  between  outcropping  arms  of  adjacent  molecules  were  computed  for  different 
shear  displacements  across  the  crystal  planes  to  demonstrate  the  occurrence  of  steric  hindrance.  In  [96],  the 
PETN  results  were  connected  with  a  similar  dislocation  model  analysis  for  RDX  of  intermolecular 
interactions  proposed  to  explain  the  formation  of  nitroso  compounds  [100].  In  a  follow-on  analysis,  Dick 
suggests  that,  because  of  the  low  stiffness  coefficients  of  PETN  compared  to  metallic,  ionic  and  covalent 
crystals,  the  theoretical  perfect  crystal  strength  might  have  been  reached  in  his  crystals  so  as  to  produce 
homogeneous  slip  and  thus  explain  equally  well  the  experimental  results  [101].  The  model  of  steric 
hindrance  was  also  applied  by  Dick  [102]  to  explain  the  explosive  behavior  reported  for  nitromethane 
crystals  when  rapidly  loaded  in  a  diamond  anvil,  high  pressure  cell  [103]. 

3.5.  Laser  damage  experiments. 

The  use  of  a  laser  beam  to  probe  damage  processes  within  individual  energetic  crystals  or  powder  compacts 
is  also  of  interest  to  the  energetic  materials  community  [104].  For  dislocation-based  interest  in  the  present 
report,  brief  descriptions  are  given  here  of  very  low  power  ignition  of  melt  spots  on  an  RDX  crystal 
surface,  in  order  to  investigate  thermal  stress  influences  on  the  combustion  process  [105],  and,  at  the 
opposite  extreme,  high  energy  input  for  initiation  of  detonation  [106,107].  Figure  15  shows,  from  the 
thermomechanically-based  combustion  study,  a  re-solidified  melt  spot  produced  with  an  impinging  laser 
beam  directed  onto  an  RDX  crystal  surface  at  a  grazing  angle  of  about  ten  degrees.  A  micrometer-scale 
network  of  cracking  is  observed  underneath  the  circumferentially-enclosed  melt  spot-crystal  interface. 
Cracking  was  reasoned  to  occur  more  easily  inside  the  melt  spot  because  only  the  liquid-solid  interfacial 
energy  had  to  be  supplied  [108]. 

A  very  different  appreciably  laser-damaged  situation  is  shown  in  Figure  16a,b,  where  heavy  damage  is 
shown  at  a  hole  (in  16b)  drilled  into  an  AP  crystal  surface  from  a  straight-on  focused,  nanosecond  pulsed, 
laser  beam  [109].  In  this  case,  the  crystal  was  laser-heated  through  the  low-temperature  orthorhombic-to 
rocksalt  cubic  phase  transformation  and  cracking  occurred  sequentially  in  both  crystal  phases.  The  initial 
macro-scale  cracking  on  the  ambient  temperature  (010)  cleavage  plane  of  the  orthorhombic  structure  is 
observed  in  Figure  16a,  for  a  magnified  region  of  16b,  while  residual  microscopic  cracking  is  observed  in 
Figure  16a  on  the  previous  {010}  planes  of  the  laser-heating-induced  transformed  cubic  (rocksalt  ) 
structure.  Figure  17  gives  a  dislocation  reaction  description  of  the  compressive  stress-state  leading  to  the 
{010}  cubic  cracking,  based  on  the  complementary  crystal  axes  shown  in  both  Figures.  The  hydrostatic 
tensile  stress  state  of  the  reacted  <100>  dislocations  was  reasoned  to  contribute  to  sub-surface  hollow 
decomposition  sites  observed  with  the  scanning  electron  microscope.  X-ray  photoelectron  spectroscopy 
(XPS)  measurements  gave  evidence  of  decomposition  to  ammonium  chlorate  and  atomic  force  microscopy 
images  gave  indication  of  transformation  product  [109]. 

4.  Chemical  reactivity 

4. 1 .  Cleavage-induced  chemical  reaction 

As  mentioned  earlier,  the  concept  of  chemical  decomposition  of  thermally  unstable  solids  resulting  from 
elastic  strain  energy  released  during  crack  propagation  breaking  intramolecular  bonds  was  introduced  by 
Bowden  et  al.  [71].  The  amount  of  decomposition  gas  generated  during  fracture  within  a  closed  chamber 
was  measured  for  P  -lead  azide,  sodium  azide,  PETN,  and  several  inert  solids,  as  given  in  Table  3.  The 
results  and  their  interpretation  were  presented  subsequently  in  more  detail  [72].  Particularly  noteworthy 
was  the  obtainment  of  estimates  of  the  effective  temperature  at  the  various  crack  tips,  calculated  assuming 
that  the  mechanism  occurring  is  the  same  as  for  conventional  thermal  decomposition. 
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In  a  separate  study,  Fox  and  Soria-Ruiz  [1 10,1 1 1]  assoeiated  differenees  in  deeomposition  kineties  of 
ealeite  as  a  funetion  of  temperature  with  the  movement  (rather  than  mere  presenee)  of  disloeations  on 
several  different  slip  planes.  Easy  deeomposition  oeeurred  on  those  planes  in  whieh  the  earbonate  (C03^') 
ion  was  favorably  oriented,  being  essentially  perpendieular  to  the  plane.  As  sueh,  this  provided  an  early 
indieation  that  solid-state  ehemieal  reaetivity  involved  partieular  erystallographie  arrangements  of  ions  (or 
moleeular  speeies,  more  generally),  rather  than  being  just  a  bulk  phenomena. 

More  reeent  fraeture-indueed  deeomposition  results  for  PETN  were  obtained  by  Hauser  et  al.  [112].  Using 
mass  speetroseopy  distinetly  different  deeomposition  produets  were  obtained  depending  on  the  rate  of 
eraek  propagation.  Slow,  low  energy  eraeks,  aehieved  by  a  maehine-driven  ehisel,  eaused  the  eleavage  of 
the  outer  moleeule  CH20-to-N02  bonds,  sueh  as  oeeurs  in  eonventional  thermal  deeomposition.  Rapid, 
high  energy  eraeks,  resulting  from  an  explosively-driven  ehisel,  eaused  the  breaking  of  the  inner  and 
stronger  C-C  bonds  at  the  eenter  of  the  moleeule.  Thus,  it  was  eoneluded  that  the  reaetion  pathway  for 
deeomposition  was  influeneed  by  the  amount  of  energy  transmitted  to  the  eraek  tip. 

4.2.  Polymerization  related  to  suerose  deeomposition. 

Enhaneed  reaetivity  at  disloeations  in  inert  (non-energetie)  solids  has  been  extensively  reviewed  by 
Thomas  [113,114].  The  emphasis  has  been  primarily  on  inorganie  materials,  but  photomierographs  were 
presented  [19]  that  indieate  eommon  aspeets  [115]  exist  between  the  initial  stages  of  solid-state 
polymerization  and  the  beginning  stages  of  deeomposition  of  moleeular  erystal  suerose,  also  assoeiated 
with  the  presenee  of  disloeations. 

4.3.  AP  deeomposition 

The  most  extensive  work  on  an  energetie  (oxidizer)  solid  to  eonneet  ehemieal  deeomposition  with  the 
presenee  of  disloeations  has  been  reported  for  single  erystal  ammonium  perehlorate,  AP.  Pioneering  results 
were  obtained  by  Galwey  and  Jaeobs  [116],  showing  that  deeomposition  of  as-grown  orthorhombie  (stable 
room  temperature  phase)  AP  erystals  is  initiated  by  defeets  and  oeeurs  primarily  at  mosaie  erystal 
boundaries.  Deeomposition  then  proeeeds  in  a  three-dimensional  manner,  eonsuming  the  boundary  regions. 
Subsequently,  Herley  and  Jaeobs  [117]  investigated  the  effeet  of  exposing  similar  erystals  to  y  -ray 
irradiation  before  thermal  deeomposition  oeeurred.  Assuming  the  same  deeomposition  meehanism,  the 
authors  eoneluded  that  irradiation  shortened  the  induetion  period  for  initial  deeomposition  while  the 
reaetion  rate  eonstant  was  found  to  be  enhaneed  for  the  aeeeleratory  period  of  deeomposition.  The  ehange 
in  rate  eonstant  was  attributed  to  an  inerease  in  the  number  of  deeomposition  nuelei  (presumably  assoeiated 
with  disloeations)  that  were  produeed  by  the  irradiation. 

Mueh  more  detailed  thermal  deeomposition  results  were  reported  by  Herley  and  Jaeobs  [1 18,1 19]  for 
millimeter  size  and  powdered  AP  erystals  unirradiated  and  exposed  to  x-  or  y  -ray  irradiation.  Previous 
observations  eoneeming  the  effeet  that  irradiation  has  on  induetion  period  and  aeeeleratory  period 
deeomposition  rate  were  eonfirmed.  Separate  mieroseopie  examination  of  partially  deeomposed  erystals 
revealed  that  the  number  of  deeomposition  nuelei  inereased  in  the  irradiated  erystals.  Irradiation  also 
inereased  the  rate  that  potential  deeomposition  sites  beeame  aetive.  The  eombined  kinetie  and 
mierostmetural  results  show  that  the  same  deeomposition  proeesses  take  plaee  in  both  unirradiated  and 
irradiated  AP,  but  that  irradiation  inereases  the  rate  eonstants. 

The  most  direet  evidenee  relating  thermal  deeomposition  sites  in  laboratory-grown  single  erystals  of  AP  to 
disloeations  was  obtained  by  Herley  et  al.  [120].  Initial  nueleation  proeesses  involved  in  establishing 
thermal  deeomposition  sites  were  investigated  using  both  light  and  (transmission  and  seanning)  eleetron 
mieroseopies.  Both  unstrained  and  plastieally  deformed  (i.e.,  from  pre-eompression)  erystals  were  heated  at 
relatively  mild  elevated  temperatures  (170-210°  C)  in  order  to  study  the  initial  stages  of  deeomposition  by 
surfaee  nueleation.  On  the  (001)  surfaee  of  unstrained  erystals,  isolated,  eireular  nueleation  sites  appear  to 
oeeur  at  random.  This  was  also  observed  for  some  sites  on  the  (210)  surfaee,  but  other  sites  were  aligned 
erystallographieally. 
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For  the  strained  crystals,  the  nucleation  sites  were  well  aligned  in  both  surfaces  and  were  mostly  in  the 
same  crystallographic  directions  as  pits  that  formed  in  chemical  etch  pitting  studies  on  separate  crystals. 
The  close  correspondence  in  the  two  sets  of  observations  led  to  the  conclusion  that  the  decomposition 
nuclei  form  preferentially  at  dislocations  that  intersect  the  surfaces.  The  initial  decomposition  nuclei  stop 
growing  when  a  critical  size  is  reached,  and  new  nuclei  (originating  a  little  below  the  surface)  form  and 
grow  until  the  top  surface  layer  has  completely  reacted.  Nucleation  of  deeper  interior  decomposition  sites  is 
also  reasoned  to  be  associated  with  dislocations  based  on  observations  that  relate  to  the  activity  of  known 
slip  systems. 

The  effects  that  x-  and  y-ray  irradiation  have  directly  on  the  chemical  decomposition  of  AP  was  also 
investigated  by  Flerley  and  Levy  [121].  As  with  thermal  decomposition,  circular  decomposition  nuclei 
(resembling  etch  pits)  formed  having  approximately  the  same  size  and  surface  density  (2-40  x  lOVcm^). 
Many  of  the  nuclei  formed  by  irradiation  were  aligned  in  the  same  crystallographic  directions  observed  in 
the  thermal  decomposition  studies  of  pre-irradiated  and  pre-stressed  crystals.  Thus,  it  was  concluded  that 
dislocations  are  also  involved  in  this  type  of  decomposition.  Irradiation-induced  decomposition  sites  form 
initially  throughout  the  exposed  volume,  in  a  roughly  uniform  distribution,  while  thermal  decomposition 
sites  appear  on  the  surface  initially  and  grow  inward  subsequently. 

Flaving  established  the  direct  involvement  that  dislocations  have  in  the  radiation-induced  decomposition  of 
AP,  Flerley  and  Levy  [122,123]  reported  measurements  of  dislocation  density  as  a  function  of  radiation 
level.  Figure  18  shows  etch  pit  observations  on  matching  crystal  surfaces  indicating  dislocation  movement 
caused  by  the  irradiation.  The  threshold  for  dislocation  generation  was  found  to  be  about  10^  rad.  Between 
10^  and  10^  rad,  the  dislocation  density  increased  linearly  with  irradiation,  and  above  10^  rad,  the  density 
was  too  high  to  measure  reliably.  It  was  concluded  that  the  observed  accelerated  decomposition  rates 
occurred  because  additional  decomposition  nuclei  formed  at  dislocations  generated  during  irradiation. 
Further,  it  was  proposed  that  interior  gaseous  decomposition  products  created  stress  levels  in  the  lattice 
high  enough  for  Frank-Read  sources  to  operate,  thus  providing  an  explanation  for  radiation-induced 
dislocation  generation.  The  interesting  results  are,  no  doubt,  very  particular  to  the  stability  of  the  AP 
molecule,  its  lattice  energy,  and  the  gaseous  decomposition  products. 

4.4.  Spectroscopy 

Tang  et  al.  [124]  have  reported  mass  spectrographic  results  for  the  decomposition  of  RDX,  PETN,  silver 
azide  and  lead  azide  crystal  surfaces  impacted  with  a  focused  laser  beam.  The  work  followed  on  from 
interest  at  the  University  of  Cambridge  Cavendish  Laboratory  [125]  in  detecting  molecular  fragments  via 
time-of- flight  mass  spectroscopy  after  fracturing  crystals  of  PETN  and  lead  azotetrazole.  The  work  relates 
as  well  to  the  results  reported  by  Fox  and  Soria-Ruiz  [1 10].  Although  significant  damage  sites  were  shown 
by  scanning  electron  microscope  observations  of  cratering  and  cracking  on  the  various  crystal  surfaces,  no 
dislocation  connection  via  slip  traces  or  other  crystallographic  markings  were  reported.  Rather,  the  results 
were  interpreted  in  terms  of  the  reaction  pathways  taken  by  the  different  molecules  in  the  decomposition 
process.  In  this  regard,  such  measurements  relate  as  well  to  the  important  consideration  of  the 
decomposition  of  energetic  crystals  in  combustion  studies  of  propellant  formulations.  Brill  [126]  has  given 
a  description  of  surface  chemistry  considerations  employing  Fourier  infrared  (FTIR)  spectrometry  to 
identify  the  chemical  reaction  steps  for  decomposition  of  AP,  RDX,  FIMX  and  other  energetic  compounds 
for  different  bum  rate  conditions.  Behrens  and  Bulusu  [127]  have  detected  early  nitroso  compound 
formation  via  mass  spectroscopy  observations  made  during  thermal  decomposition  and  pointed  out  that  the 
measurments  were  consistent  with  decomposition  in  the  solid  phase  as  reported  also  from  mechanical  test 
results  [128]. 

4.5.  X-radiation  induced  decomposition 

Synchrotron  x-ray  topography  was  used  by  Bhat  et  al.  [129]  to  investigate  beam-induced  low-level 
decomposition  (about  1%)  in  large  AP  single  crystals.  Significant  stractural  damage  was  observed  in  the 
outer  layers  of  the  crystal  causing  bending  to  occur  by  slip  on  previously  reported  systems.  An 
enhancement  in  the  decomposition  rate  was  reported  with  dislocations  serving  as  the  decomposition  nuclei. 
The  precipitation  of  decomposition  products  at  the  nuclei  causes  crystals  (initially  transparent)  to  become 
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opaque  in  a  process  that  resembles  the  formation  of  a  colloidal  suspension.  Importantly,  the  synchrotron 
white  beam  is  able  to  image  the  resultant  defect  structure  at  higher  levels  of  decomposition  than  was 
possible  in  conventional  laboratory  experiments  using  monochromatic  radiation. 

The  use  of  x-ray  topography  to  examine  sub-structural  changes  occurring  during  chemical  reactions  in 
several  crystalline  solids,  including  single  crystal  AP,  was  reported  by  Sheen  and  Sherwood  [130].  The 
suitability  of  using  laboratory  and  synchrotron  sources  to  study  solid-state  decomposition  was  characterized 
in  more  detail.  Using  white  beam  synchrotron,  transmission  Lane  diffraction  topographs  of  AP  were 
obtained  as  a  function  of  absorbed  doses.  Graded  absorption  of  the  available  wavelengths  across  the  crystal 
thickness  caused  the  crystal  to  bend  initially  with  the  development  of  a  large  amount  of  dislocation 
structure.  Eventually,  the  strain  equalized  and  flexure  of  the  crystal  was  eliminated,  followed  by  detection 
of  chemical  decomposition  based  on  image  broadening  and  subsequent  formation  of  product  powder  rings. 

4.6.  Shock  reactivity 

The  roles  that  deformation,  fracture,  and  material  microstructure  have  on  the  shock  reactivity  of  single 
crystal  AP  were  investigated  by  Sandusky  et  al.  [91].  Experiments  were  conducted  in  a  mineral-oil  bath 
with  the  shock  being  provided  by  a  small  explosive  donor.  This  arrangement  allowed  samples  to  be 
recovered  for  subsequent  chemical  analysis  and  microstructural  characterization.  Particularly  noteworthy 
experiments  involved  shocking  crystals  possessing  a  large  surface  strain  center  created  by  a  Vickers 
indenter.  High-speed  photographic  coverage  revealed  distinct  diagonal  lines  immediately  behind  the  shock 
front  corresponding  to  (010)[001]  slip,  a  moving  luminous  band  believed  associated  with  a  propagating 
crack,  and  light  near  to  the  hardness  impression  that  was  conjectured  to  be  due  to  chemical  reaction. 

A  subsequent  experiment  was  conducted  [131]  on  an  AP  crystal  containing  several  Vickers  impressions 
and  shocked  at  the  reaction  threshold.  High-speed  photographs  again  showed  luminosity  near  impressions. 
The  recovered  crystal  was  cleaved  through  two  of  the  impressions,  allowing  spatial  chemical  analysis  to  be 
performed  in  the  crystal  interior  by  using  the  XPS  technique  to  scan  1 .0  x  1.0  mm  areas.  The  line  widths  of 
XPS  spectra  were  correlated  [92]  with  dislocation  densities  that  were  measured  directly  using  a  chemical 
etch  pit  technique  and  related  to  somewhat  analogous  irradiation  results  obtained  by  Herley  and  Levy  [121- 
123]  that  were  discussed  earlier.  The  linewidth  of  the  Cl(2/73/2)  spectra  ranged  from  1.70  eV  for  the  region 
of  greatest  visible  damage  to  1 .22  eV  for  the  region  with  no  visible  damage  or  for  an  as-grown  (control) 
crystal. 

The  broadening  of  the  photoelecfron  lines  was  attributed  to  the  lattice  distortion  present  in  the  immediate 
vicinity  of  the  dislocations.  On  the  same  freshly  cleaved  surfaces,  the  largest  amounts  of  chemical 
decomposition  were  associated  with  the  hardness  impressions.  However,  the  highest  levels  were  not 
detected  in  the  areas  immediately  adjacent  to  the  impressions  where  dislocation  densities  were  found 
(based  on  line-width  measurements)  to  be  highest.  Rather,  the  highest  decomposition  levels  were  found  in 
areas  where  the  dislocation  densities  had  intermediate  values  and  appeared  to  be  associated  with  a 
prominent  (001)  crack  or  irregular  cleavage  steps  emanating  from  an  impression.  It  was  concluded  that  the 
observed  luminosity  resulted  from  chemical  reaction  occurring  either  at  the  tip  of  the  (001)  crack  or  in  its 
plastic  zone. 

5.  Molecular  modeling 

The  various  methods  of  molecular  modeling  have  been  applied  mainly  to  prediction  of  crystal  structures  and 
densities,  determination  of  the  P-N  stress,  dislocation  displacements  relating  to  intermolecular  interactions,  and 
shock  wave  properties,  including  electronic  or  vibrational  excitations  and  shear  wave  characteristics. 

5.1.  Crystal  densities 

Prediction  of  energetic  crystal  densities  is  important  because,  in  addition  to  interest  in  greatest  energy  density 
for  explosiveness,  both  the  detonation  pressure  and  velocity  increase  with  density  [132].  The  topic  relates  to  x- 
ray  and  neutron  diffraction  methods  employed  in  crystal  structure  determinations.  Computer  assisted  modeling 
and  graphics  programs  are  employed  with  energy  minimization  techniques.  Stine  [133]  has  summarized 
empirical  methods  for  predicting  densities  and  heats  of  formation  of  proposed  high  explosives.  Gilardi  [134] 
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has  given  a  crystal  structure -based  description  of  interest  in  polycyclic  energetic  compounds  for  enhancing 
performance  and/or  safety.  An  interesting  connection  with  such  modeling  efforts  is  provided  by  the  description 
given  by  Armstrong  et  al.  [135]  for  deformation  twinning  in  the  monoclinic  beta-HMX  crystal  structure  where 
the  crystal  structure  derived  potential  energy  was  followed  in  identifying  a  type  II  {101}<101>  twinning 
mechanism;  see  Figure  19.  Twinning  was  concluded  to  be  favorable  over  slip  in  the  beta-FIMX  structure 
because  of  the  flexibility  of  the  FIMX  molecule.  Density  considerations  have  been  investigated  on  a  molecular 
fragmentation  basis  also  so  as  to  understand  proposed  explosive  versus  non-explosive  behaviors  of  related 
molecular  structures  [136]. 

5.2.  The  P-N  stress  for  anthracene 

The  computer  simulation  of  an  [010](001)  edge  dislocation  in  anthracene  has  been  reported  [137].  The 
structure  relates  to  that  described  by  Robinson  and  Scott  [37]  for  Figure  7.  Figure  20  shows  equipotential 
contours  in  the  (001).  A  dislocation  width  of  6.2  times  the  Burgers  vector  was  obtained,  indicative  of  a 
relatively  low  P-N  stress  for  the  material.  The  hardness  stress-strain  result  indicated  for  anthracene  in  Figure  10 
is  consistent  with  a  low  shear  strength  being  observed  for  the  material,  presumably  because  of  only  minor 
intermolecular  blockages  occurring  in  the  slip  process.  An  interesting  consideration  in  regard  to  the 
intermolecular  blockage  of  adjacent  molecules  is  the  suggestion  of  Nabarro  that  taking  account  of  the 
displacement  normal  to  the  slip  plane  in  evaluating  the  P-N  stress  should  give  a  more  accurate  value.  Lee  and 
Dundurs  [138]  found  in  evaluating  the  suggestion  on  a  conventional  basis  that  the  stress  was  lowered. 

5.3.  Model  dislocation  displacements 

A  schematic  picture  is  shown  in  Figure  21  for  a  [100]  Burgers  vector  dislocation  with  line  direction  along  [001] 
one  half-cell  deep  in  an  RDX  crystal.  The  dislocation  has  been  reasoned  to  slip  on  (040)  because  of  the 
minimized  intermolecular  interferences  for  shear  displacements  during  movement  [100].  The  indicated 
molecular  positions  at  the  dislocation  core  are  clearly  far  out  of  order  relative  to  normal  nearest  neighbor 
associations.  Based  on  the  consideration  that  dislocation  cores  in  such  molecular  crystals  might  be  hollow,  as 
predicted  by  Frank  [18],  see  section  1.2,  the  hollow  core  radius  computed  for  RDX  is  greater  than  7  nm,  that  is 
about  6  unit  cell  distances,  thus  relating  as  well  to  the  relatively  large  values  of  the  dislocation  self-energies. 

The  issue  of  molecular  interferences  associated  with  dislocation  movement  is  demonstrated  as  well  for  the 
RDX  lattice  stracture  by  examination  of  slip  across  the  (021)  slip  plane  that  was  identified  at  hardness  micro¬ 
indentation  sites  by  Elban  et  al.  [139].  Figure  22  gives  a  (100)  view  of  the  molecular  blockages  occurring  for 
shearing  across  the  (021)  slip  plane,  in  and  out  of  the  projected  plane  of  observation,  along  [100].  The  nitrogen- 
oxygen  interactions  indicated  by  the  cross-hatched  areas  of  overlapping  molcecules  when  displaced  along  [100] 
are  just  those  favoring  the  creation  ofnitroso  compounds  detected  in  "go/no-go"  drop-weight  impact  tests  [128] 
and  in  thermal  decomposition  experiments  [127]. 

Dick  and  Ritchie  [140]  followed  up  also  on  the  earlier  work  by  Dick  and  colleagues  [94-96],  and  in  extension 
of  individual  reports  of  their  own  [141,142],  to  evaluate  the  importance  of  molecular  interferences,  or  steric 
hindrances,  in  causing  shorter  run-to-detonation  distances  in  PETN  crystals  shocked  along  different  directions, 
that  is,  the  distance  along  a  crystal  rod  traversed  by  a  shock  wave  before  detonation  is  initiated.  Also,  elastic 
precursor  stress  levels  specifying  different  shock  strengths  were  measured  along  the  different  directions  using  a 
light  gas  gun  facility.  Unhindered  [100]  and  [101]  directions  were  tested  along  with  hindered  [110]  and  [001] 
tests.  The  elastic  precursor  strength  increased  in  that  same  order.  The  results  are  in  agreement  with  the 
reasoning  for  dislocation  passage  being  blocked  by  interferingly-close  intermolecular  distances.  Intermolecular 
interactions  were  evaluated  with  a  molecular  mechanics  program  AMBER  develped  for  biochemical  systems. 
Snapshots  of  the  differently  sheared  lattice  were  presented.  The  computations  were  stated  to  support  the  elastic 
precursor  measurements;  and,  intramolecular  bond  rotations  and  angular  strains  were  found  to  be  relatively 
important  as  compared  with  little  changes  occurring  in  bond  lengths. 

5.4.  Flot  Spots  in  MD  simulations 

Molecular  dynamics  (MD)  methods  lend  themselves  to  shock  wave  modeling  considerations  because  of  the 
nm-scale  sharpness  of  the  shock  front  and  presumed  femtosecond-to-picosecond  time  scale  of  fundamental 
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events.  Karo  and  eolleagues  performed  pioneering  eomputations  relating  to  energy  partitioning  and  transfer 
proeesses,  thermal  versus  athermal  eharaeteristies,  bond  mpturing,  and  defeet  influenees  [143].  Tsai  [144,145] 
has  deseribed  the  influenee  of  a  eluster  of  vaeaneies  on  enhaneing  hot  spot  development  through  struetural 
relaxation  at  the  elustered  defeet  site.  During  sueh  relaxation,  apparent  slip  bands  developed  at  the  eluster  site 
and  spread  aeross  the  speeimen.  Bandak  et  al.  [146]  elaborated  on  the  model.  Figure  23  shows  an  (001)  body- 
eentered  eubie  lattiee  eell  of  Morse-potential-type  atoms,  also  eontaining  a  ten-vaeaney  eluster  that  relaxed 
under  the  MD  proeess  of  rapid  eompression  to  form  a  pair  of  disloeation  dipoles.  The  stress-strain  behavior  of 
the  system  was  followed  to  show  apparent  yielding  during  relaxation  of  the  vaeaney  eluster.  In  a  follow-up 
study  of  a  simple  moleeular  lattiee  eontaining  vaeant  moleeules,  Tsai  and  Armstrong  emphasized  that  sueh 
permanent  relaxation  was  assoeiated  with  eonversion  of  the  higher  potential  energy  of  moleeules  adjaeent  to  the 
defeet  strueture  to  thermal  hot  spot  energy  [147].  The  important  point  is  that  signifieant  hot  spot  heating  is 
always  assoeiated  with  meehanieal  relaxation  and  so  a  hot  spot  influenee  is  always  present  at  sueh  sites.  The 
same  eonelusion  was  obtained  for  adaptation  of  the  MD  method  to  a  lattiee  of  eight-membered  moleeular  units 
[148]. 

The  issue  of  an  MD  deseription  being  appropriate  for  the  eomprehensive  deseription  of  the  generation  and 
propagation  of  shoek  waves,  eovering  the  gamut  from  mateh  with  the  eontinuum  theory  of  planar  detonation  to 
eonneetion  with  loeal  defeet  properties,  as  deseribed  above,  was  taken  up  by  C.T.  White  and  eolleagues  [149]. 
The  assoeiated  hot  spot  meehanism  of  void  eollapse  was  eonfirmed  for  a  model  diatomie  lattiee  [150].  Other 
MD  aspeets  have  been  reported  for  shoek  eompressions  and  hot  spot  properties,  as  follows:  model  assessment 
of  the  multi-phonon  up-pumping  time  seale  [151],  vibrational  heating  of  naphthalene  moleeules  [152],  and 
pieoseeond  MD  resolution  of  events  behind  a  modeled  shoek  front  in  polymethylmethaerylate  to  eompare  with 
dye-monitored  experimental  measurements  [153]. 

5.5.  Shoek  meehanoehemistry  and  atomie  foree  mieroseopy 

Gilman  [154]  has  pointed  to  historieal  evidenee,  stemming  from  Bridgman's  experiments,  for  the  importanee  of 
shear  deformations  in  promoting  ehemieal  reaetions.  Based  on  an  analysis  of  detonation  front  veloeity,  Gilman 
argues  that  the  front  propagation  is  too  fast  for  thermally-aetivated  ehemieal  reaetions  to  oeeur  [155]  and 
homogeneous  shear-indueed  eleefronie  exeitations  akin  to  metallization  then  provide  the  meehanism  for 
ehemieal  deeomposition.  The  issue  relates  to  earlier  diseussion  [85,76]  of  the  various  stages:  disloeation 
movement  behind  a  shoek  front  eontributing  to  initiation  of  deeomposition  for  weak  shoeks;  that  disloeation 
strueture  moving  behind  the  shoek  front,  for  stronger  shoeks,  and  interaeting  with  nanoseale  disloeation  elusters 
formed  at  the  front;  and,  to  "in-situ"  disloeation  eore  enhaneement  of  deeomposition  for  even  stronger  shoeks. 
In  prineiple,  the  importanee  of  anharmonie  eoupling  to  the  meehanism  of  vibrational  up  pumping  so  as  to 
provide  ehemistry  in  the  shoek  front  is  enhaneed  within  the  strain  field  of  the  disloeation  eores  [89].  The  Frank 
hollow  eore  eonsideration  [18]  that  appears  to  apply  at  the  nm-seale  for  energetie  erystals  should  be  involved. 

Reeent  results  on  energetie  erystal  surfaees  obtained  with  atomie  foree  mieroseopy  (AFM)  by  Sharma,  and 
analyzed  with  eolleagues,  may  help  to  further  quantify  deformation  and  thermal  influenees  on  energetie 
material  deeompositions,  for  example,  relating  to  reeent  deseription  ofpost-shoek  examination  of  RDX  erystals 
having  undergone  an  allotropie  phase  transformation  [156].  Sub-moleeular-sized  shear  eraek  displaeements 
were  deteeted  within  nanometrie-sized  erystals  formed  by  multiple  fraeturing  of  the  initial  bulk  erystal 
speeimen  [157].  Reeent  measurements  have  ineluded  examination  of  nanoindentations  made  within  the  AFM 
[158].  Figure  24  shows  an  AFM  pieture  of  the  RDX  erystal  lattiee  at  nm-seale  resolution  of  individual 
moleeule  "spots"  [159].  Other  AFM  examinations  of  laboratory-grown  RDX  erystals,  in  an  as-eleaved 
eondition  and  after  various  deformations  and/or  heat  treatments,  have  shown  differently  oriented  surfaee  step 
struetures  and  elustered  zones  of  disordered  moleeules. 

6.  Summary 

The  enhaneed  visibility  of  steps  assoeiated  with  the  large  Burgers  veetors  of  disloeations  in 
energetie/moleeular  erystals  is  aided  indireetly  by  the  oeeurrenee  of  low  disloeation  densities  also 
assoeiated  with  the  large  Burgers  veetors  that  lead  to  relatively  high  self-energies.  A  signifieant  variation  in 
disloeation  Burgers  and  line  veetors  oeeurs  during  erystal  growth  beeause  of  the  low  symmetry  of  the 
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molecular  crystal  structures,  conformational  requirements  of  positioning  oddly-shaped,  relatively  large, 
molecules,  and  their  relatively  weak  intermolecular  bonding. 

Dislocation  movement,  whether  in  static  indentation  tests,  under  drop-weight  impact,  or  in  shock  wave 
loading,  is  relatively  hindered  by  intermolecular  cross-blockages  of  adjacent  juxtaposed  molecules,  thus 
giving  plastic  flow  stress  values  near  to  the  theoretical  limit  for  cleavage  cracking.  The  high  self-energies 
of  the  dislocations  make  dislocation  nucleation  prohibitive  while,  once  present,  their  strain  fields  can  assist 
chemical  decomposition  under  suitable  thermal  and  environmental  circumstances. 

Local  concentration  of  strain  energy  within  the  relatively  more  stable  energetic  crystals  is  required  to 
induce  intramolecular  decomposition,  relating  to  (intermolecular-type)  hydrogen  and  (intramolecular) 
nitrogen/oxygen  bond  breakages,  for  example,  as  indicated  in  Figures  19,  21,  and  22.  Dislocation  pile-ups 
provide  an  explanation  of  crystal  size  effects  in  energetic  crystal  decompositions.  Such  molecular  crystals 
give  relative  enhancement  of  local  heating  from  plastic  work  because  of  their  insulating  character. 

An  interesting  recent  connection  is  with  the  thermomechanical  aspect  of  energetic  crystal  combustion  when 
serving  as  an  ingredient  of  a  propellant  motor.  Such  combustion  occurs  over  a  melt  layer  formed  on  the 
individual  crystal  surfaces.  Micro-cracking,  as  indicated  in  Figure  15,  is  enhanced  at  the  liquid-solid 
interface  because  of  need  to  supply  only  the  lower  liquid-solid  interfacial  energy  of  cracks,  as  compared 
with  the  far  larger  solid-vapor  value  [108].  Recent  model  fracture  mechanics  estimates  of  the  pressures 
required  to  propagate  such  micro-crack  networks,  relating  to  the  indentation  fracture  mechanics 
measurements  in  Figure  9,  give  indication  of  possible  crystal  shattering  instabilities  for  combustion  under 
high  pressure  conditions  [160]. 

The  role  of  dislocations  in  shock  wave  initiation  of  detonations  is  not  resolved.  Possible  mechanisms  begin 
with  dislocation  activity  following  behind  the  shock  front  for  low  shocks,  shock  front  generation  of 
dislocations  for  greater  shocks  also  connecting  with  follow-on  movement  of  the  previously-resident 
dislocation  structure,  and,  finally,  “in-situ”  initiation  at  dislocation  cores  for  high  shocks. 
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Figure  9.  Logarithmic  dependence,  according  to  constant  hardness  or  constant  indentation  fracture  stress 
intensity,  of  applied  indenter  force  on  diagonal  indentation  and  crack  lengths  for  RDX  and  MgO  crystals; 
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Figure  10.  Flardness  stress-strain  basis  for  comparing  the  elastic,  plastic  and  cracking  behaviors  of  crystal 
materials,  including  calculation  of  the  terminal  elastic  cracking  stresses  for  different  ball  sizes;  after 
Flammond  and  Armstrong  [49]. 

Figure  1 1 .  Release  of  a  dislocation  pile-up  from  its  blocking  obstacle,  showing  the  shape  of  the  dislocation 
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including,  for  0.25  mm  size  crystals,  an  indication  of  influence  of  variation  of  crystal  quality  in  different 
batches,  after  Armstrong,  Coffey,  DeVost  and  Elban  [70]. 

Figure  13.  Graphical  basis  for  illustrating  hot  spot  susceptibility  of  materials  on  a  dislocation  pile-up 
avalanche  model  description,  including  different  hot  spot  temperature  expressions  for  metallic  and  ionic 
materials  as  compared  with  energetic  RDX  and  PETN,  after  Armstrong,  Coffey  and  Elban  [60]. 

Figure  14.  Shock  model  for  nanometer-scale  dislocation  loops  that  produce  a  residual  deformation  state  of 
one-dimensional  compression,  after  Armstrong,  Miller,  and  Sandusky  [79],  as  elaborated  by  Bandak, 
Armstrong,  and  Douglas  [80]. 

Figure  15.  Combustion  residue  and  microcracking  at  an  RDX  (210)  crystal  surface  bum  spot  after  laser 
irradiation;  the  vertical  separation  of  parallel  microcracks  is  ~  5  micrometers.  After  Armstrong, 
Ramaswamy,  and  Field  [105],  as  further  analyzed  by  Elban,  Armstrong,  and  Russell  [108]. 

Figure  16.  Laser-impacted  (OOl)AP  crystal  surface  (a)  at  high  magnification  to  show  the  sequential 
macroscopic  (210)  cracking  in  the  low  temperature  orthorhombic  phase  and  higher  temperature,  finer  scale, 
{OlOjcracking  retained  from  the  rocksalt-type  cubic  phase;  and  (b)  lower  magnification  view  of  the  central 
drilled  hole  exhibiting  a  whitish  halo  of  sub-surface  (001)  cracking.  After  Ramaswamy,  Shin,  Armstrong, 
Lee,  and  Sharma  [109]. 

Figure  17.  Schematic  dislocation  model  for  explanation  of  the  {010}  cracking  shown  in  Figure  16a  for  the 
higher  temperature  rocksalt-type  cubic  stmcture  of  AP;  after  Ramaswamy  et  al.  [109]. 

Figure  18.  Photomicrographs  of  matching  AP  (210)  cleavage  surfaces  showing  radiation-induced 
dislocation  movement.  Side  a  was  not  irradiated  but  side  b  received  a  dose  of  10^4  rad  gamma-ray 
exposure.  A  set  of  dislocations  which  have  moved  during  irradiation  is  labeled  M  while  set  S  has  remained 
stationary.  After  Flerley  and  Levy  [122,123]. 

Figure  19.  Crystallography  of  Type  II  {I01}<101>  twinning  in  HMX;  on  a  stereographic  projection  basis 
at  top  left,  and  (010)  plane  projection  of  pre-and  post-twin  positions  at  top  right,  with  tabulated  sequential 
potential  energy  evaluations  for  the  indicated  lattice  point  displacements,  also  involving  required  N-N  bond 
rotation.  After  Armstrong,  Ammon,  Du,  Elban  and  Zhang  [135]. 

Figure  20.  A  projection  in  the  (100)  plane  of  molecular  centers  at  an  edge  dislocation  core  in  anthracene 
after  relaxation  of  the  lattice  positionsand  with  equipotential  curves  (in  units  of  eV).  The  crosses  denote 
comer  molecules  and  circles,  center  molecules.  After  Ide,  Okada,  and  Kojima  [137]. 

Figure  21.  Schematic  view  of  RDX  molecule  positions  for  an  [001]  edge  dislocation  on  the  (040)  plane, 
after  Armstrong  [28]. 

Figure  22.  For  RDX,  dinifroso-related  interactions  from  [100](021)  slip,  in  and  out  of  the  surface-projected 
(100);  after  Armstrong  [28]  and  relating  to  experimental  results  both  from  Hoffsommer,  Glover,  and  Elban 
[128]  and  Behrens  and  Buluso  [127]. 
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Figure  23.  Molecular  dynamics  modeling  of  the  collapse  of  a  vacancy  cluster  into  a  pair  of  edge 
dislocation  dipoles,  after  Bandak,  Tsai,  Armstrong,  and  Douglas  [146]. 

Figure  24.  Atomic  force  microscope  image  of  molecule  positions  on  an  identified  RDX  (001)  cleavage 
surface,  after  Sharma,  unpublished;  see  J.  Sharma  et  al.  [159]. 
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Table  1.  Comparison  of  dislocation  strain  energies  and  Peierls-Nabarro  Stresses  for  RDX,  LiF,  and 
MgO,  after  Armstrong  and  Elban  [22]. 


Material 

Dislocation 
Strain  Energy 
Coefficient 

(El/p). 

J 

i  'l 

Peierls- 
Nabarro 
Stress,  (Tp), 
N/mm2 

VAHf/nry 

Ig/ 

RDX 

6.3x10-19 

12 

580 

0.077 

LiF 

0.9x10-19 

2.1 

3100 

0.050 

MgO 

2.4x10-19 

1.9 

8000 

0.054 

E,,=(Gb2Ae/4;rflr) 

rp(=  Tp.n)  =  [2G/1-v)]exp[-2;rd/(1-  v)b] 


Table  2.  Tbe  types  of  pressure  -  time  curves  sbowu  by  various  materials  subjected  to  impact,  after 
Heaveus  aud  Field  [58]. 


NO  PRESSURE  SLIGHT  PRESSURE  SERIES  OF  PRESSURE 
OISTURBANCE  OISTURBANCE  FLUCTUATIONS 


SHARP  OROP  IN 
PRESSURE 


SHARP  OROP. 
PRESSURE  FALLING 
TO  ZERO 


CHARCOAL 

NaCI 

TETRYLt 

GRAPHITE 

KBr 

PICRIC  ACID 

GLASS 

BORAX 

COMPOSITION  B 

CaCOj 

Hgl 

OXALIC  ACID 

NajCOj 

UF 

POLYVINYL  ALCOHOL 

AMMONIUM 

OXALATE 

CELLULOSE 

AMMONIUM 

(SAWDUST) 

ACETATE 

NGt 

NCt 

BLASTINGt 

GELATINE 

NUMEROUS 
ORGANIC  SOLIDS 

NH4CIO4* 

PETN* 

NaNj 

ROX* 

BaNOj 

HMX* 

SULPHUR 

TNT 

PHENOLPHTHALEIN 

NH4NO3 

KNO3 

KN3 

NH4CI 

TARTARIC  ACID 

SUCROSE 

GUNPOWDER 

POLYSTRYENE 

*  THESE  MATERIALS  IGNITED  AT  THE  INSTANT  OF  THE  PRESSURE  DROP 
t  THESE  MATERIALS  IGNinO.  BUT  NOT  AT  THE  INSTANT  OF  A  PRESSURE  DROP 
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Table  3.  Fracture-induced  cleavage  surface  chemical  decompositions,  after  Bowden,  Fox,  and  Soria- 
Ruiz  [71]  and  Fox  and  Soria-Ruiz  [72]. 


Material 

Magnesium  oxide 
Calcium  carbonate  (calcite) 
Magnesium  carbonate  (magnesite) 
Lead  carbonate  (cerussite) 

Sodium  azide 
/}-Lead  azide 
PETN 


No.  of 

Amount  of 

Average  No. 

decomposition 

of  monolayers 

expts. 

(molecules  cm-*) 

decomposed 

9 

<12« 

<  0001 

14 

2  X  10>» 

0-09 

14 

1  X 10“ 

0-37 

8 

~  7  X  10“ 

~2 

5 

1  X 10“ 

0-02 

10 

~3xl0“ 

~10 

7 

2  X  10“ 

(for  mass  28) 

~2 
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Figure  1.  A  phase  contrast  micrograph  of  the  (0001)  face  on  a  6H  SiC  crystal  showing  the  growth 
pattern  due  to  three  dislocations  of  like  sign  close  together  and  giving  rise  to  three 
cooperating  growth  spirals.  The  central  or  the  last  part  of  the  growth  is  circular  while  the 
rest  is  hexagonal.  (X500).  After  Verma  [1]. 


Figure  2.  An  (0006)  section  x-ray  transmission  topograph  showing  an  image  of  a  screw  dislocation 
micropipe,  M;  after  Dudley,  Si,  Wang,  Carter,  Glass,  and  Tsvetkov  [20]. 
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Figure  3.  Schematic  illustration  of  (a)  division  of  a  crystal  into  growth  sectors  according  to  different 
growth  directions  including  indication  of  (dashed)  contours  at  different  stages  of  growth; 
and  (h)  typical  geometry  of  grown-in  dislocation  lines  initiated  at  inclusion  particles.  After 
Klapper  [26]. 


Figure  4.  Growth  sector  structure  within  a  sectioned  RDX  crystal,  after  van  der  Steen  and  Duvalois 
[27],  with  growth  hand  (planar)  trace  identifications,  after  Armstrong  [28]. 
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Figure  5.  A  {210}  x-ray  transmission  topograph  of  an  RDX  crystal  showing  dislocation  bundles 
spreading  from  a  central  seed  position  to  the  external  crystal  growth  surfaces;  after 
McDermott  and  Phakey  [33]. 


Figure  6.  A  scanning  electron  micrograph  of  an  RDX  crystal,  after  van  der  Steen  and  Duvalois  [27], 
with  crystallographically-identified  planar  growth  surfaces,  also  with  A  to  D  locations 
marked  for  growth  protuberances  at  positions  of  dislocation  bundle  emergence,  after 
Armstrong  [28]. 
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Figure  7.  Schematic  plane  sections  of  the  crystal  structure  of  anthracene:  (a)  (010)  projection  of  the 
unit  cell,  with  a  and  c  axes  identified,  and  showing  carhon  atoms  and  carhon-carhon 
bonds;  (h)  (001)  projection,  with  a  and  h  axes  marked,  and  showing  carhon-carhon  bonds, 
with  the  light  lines  applying  for  molecules  at  height  c  above  the  dark  Une  molecules;  (c)  a 
faulted  (001)  structure  relating  to  (b).  After  Robinson  and  Scott  [37]. 
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Figure  8.  Dislocation  etch  pits  centered  on  microindentations  in  the  (210)  plane  growth  surface  of 
RDX:  (top)  Knoop  indentation  with  long  axis  parallel  to  [001];  and  (bottom)  identified 
Vickers  diamond  pyramid  indentation  with  one  equiaxed  diagonal  length  parallel  to  [001  ]. 
The  long  hlack  lines  in  the  figures  are  cracks.  After  Elhan  and  Armstrong  [24]. 


139 


FORCE, 


Figure  9.  Logarithmic  dependence,  according  to  constant  hardness  or  constant  indentation  fracture 
stress  intensity,  of  applied  indenter  force  on  diagonal  indentation  and  crack  lengths  for 
RDX  and  MgO  crystals;  after  Armstrong  and  Elhan  [22]. 


140 


d/D 

Strain,  (indentation  dia./ball  dia.) 


Figure  10.  Hardness  stress-strain  basis  for  comparing  the  elastic,  plastic  and  cracking  behaviors  of 
crystal  materials,  including  calculation  of  the  terminal  elastic  cracking  stresses  for 
different  ball  sizes;  after  Hammond  and  Armstrong  [49]. 
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Figure  11.  Release  of  a  dislocation  pile-up  from  its  blocking  obstacle,  showing  tbe  shape  of  the 

dislocation  distribution  as  it  moves  from  right  to  left,  at  sequential  (dimensionless)  time 
steps,  t  =  0,  0.01,  0.03,  0.10  and  0.30,  after  Head  [  61]. 


Figure  12.  Drop-weight  impact  height  for  50%  probability  of  initiation  as  it  depends  on  crystal  size, 
including,  for  0.25  mm  size  crystals,  an  indication  of  influence  of  variation  of  crystal 
quality  in  different  batches,  after  Armstrong,  Coffey,  DeVost  and  Elban  [70]. 
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Figure  13.  Graphical  basis  for  illustrating  hot  spot  susceptibiUty  of  materials  on  a  dislocation  pile-up 
avalanche  model  description,  including  different  hot  spot  temperature  expressions  for 
metallic  and  ionic  materials  as  compared  with  energetic  RDX  and  PETN,  after 
Armstrong,  Coffey  and  Elban  [60]. 


Figure  14.  Shock  model  for  nanometer-scale  dislocation  loops  that  produce  a  residual  deformation 
state  of  one-dimensional  compression,  after  Armstrong,  Miller,  and  Sandusky  [79],  as 
elaborated  by  Bandak,  Armstrong,  and  Douglas  [80]. 
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t,  -0 


Figure  15.  Combustion  residue  and  microcracking  at  an  RDX  (210)  crystal  surface  burn  spot  after 
laser  irradiation;  tbe  vertical  separation  of  parallel  microcracks  is  ~  5  micrometers. 
After  Armstrong,  Ramaswamy,  and  Field  [105],  as  further  analyzed  by  Elban, 
Armstrong,  and  Russell  [108]. 


Figure  16.  Laser-impacted  (OOl)AP  crystal  surface  (a)  at  high  magnification  to  show  the  sequential 
macroscopic  (210)  cracking  in  the  low  temperature  orthorhombic  phase  and  higher 
temperature,  finer  scale,  {010}cracking  retained  from  the  rocksalt-type  cubic  phase;  and 
(b)  lower  magnification  view  of  the  central  drilled  hole  exhibiting  a  whitish  halo  of  sub¬ 
surface  (001)  cracking.  After  Ramaswamy,  Shin,  Armstrong,  Lee,  and  Sharma  [109]. 
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DISLOCATION  MECHANISM  FOR  CRACKING.  R=ROCKSALT 
0=0RTH0RH0MBIC 


Figure  17.  Schematic  dislocation  model  for  explanation  of  the  {010}  cracking  shown  in  Figure  16a 
for  the  higher  temperature  rocksalt-type  cubic  structure  of  AP;  after  Ramaswamy  et  al. 
[109]. 
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Figure  18.  Photomicrographs  of  matching  AP  (210)  cleavage  surfaces  showing  radiation-induced 
dislocation  movement.  Side  a  was  not  irradiated  hut  side  h  received  a  dose  of  lO'^d  rad 
gamma-ray  exposure.  A  set  of  dislocations  which  have  moved  during  irradiation  is 
labeled  M  while  set  S  has  remained  stationary.  After  Herley  and  Levy  [122,123]. 
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Position 

A 

(TIPE*)^^ 

Kcal/mol 

Position 

A 

(T2PE**)^^ 

Kcal/mol 

0,000 

-41.97 

0.000 

-29.95 

0.265 

-40.91 

0,265 

-32.43 

0.530 

-36.79 

0.530 

-35.91 

0.790 

-36.40 

0.790 

-37.58 

1.060 

-34 , 19 

1.060 

-38.06* 

*  Type  1  Potential  Energy 
**  Type  2  Potential  Energy 

4.  Type  2  Potential  Energy  +  21°  rotation  of  N-N  bond 

Figure  19.  Crystallography  of  Type  II  {101}<101>  twinning  in  HMX;  on  a  stereographic  projection 
basis  at  top  left,  and  (010)  plane  projection  of  pre-and  post-twin  positions  at  top  right, 
with  tabulated  sequential  potential  energy  evaluations  for  the  indicated  lattice  point 
displacements,  also  involving  required  N-N  bond  rotation.  After  Armstrong,  Ammon, 
Du,  Elban  and  Zhang  [135]. 


146 


Figure  20.  A  projection  in  the  (100)  plane  of  molecular  centers  at  an  edge  dislocation  core  in 

anthracene  after  relaxation  of  the  lattice  positionsand  with  equipotential  curves  (in  units 
of  eV).  The  crosses  denote  corner  molecules  and  circles,  center  molecules.  After  Ide, 
Okada,  and  Kojima  [137]. 


Figure  21.  Schematic  view  of  RDX  molecule  positions  for  an  [001]  edge  dislocation  on  the  (040) 
plane,  after  Armstrong  [28]. 
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Figure  22.  For  RDX,  dinitroso-related  interactions  from  [100](021)  slip,  in  and  out  of  the  surface- 
projected  (100);  after  Armstrong  [28]  and  relating  to  experimental  results  both  from 
Hoffsommer,  Glover,  and  Elban  [128]  and  Behrens  and  Buluso  [127]. 
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Figure  23.  Molecular  dynamics  modeling  of  the  collapse  of  a  vacancy  cluster  into  a  pair  of  edge 
dislocation  dipoles,  after  Bandak,  Tsai,  Armstrong,  and  Douglas  [146]. 
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Figure  24.  Atomic  force  microscope  image  of  molecule  positious  ou  au  ideutified  RDX  (001)  cleavage 
surface,  after  Sharma,  uupublished;  see  J.  Sharma  et  al.  [159]. 
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CRYSTAL  IMPERFECTIONS  SEEN  BY  X-RAY 
Dll' FRACTION  TOPOGRAPHY 


R.W.  Armstrong 

AFRL/MNMli,  2306  Perimeter  Road,  Eglin  Air  Force  Base,  FL  32542-5910 


Abstract 

Internal  lattice  strains  imaged  via  x-ray  diffraction  topography  (XRDT)  methods  are  presented  for 
relatively  low  dislocation  density,  or  dislocation-free,  crystals,  much  as  investigated  for  higher 
dislocation  density  materials  by  Professor  Gareth  Thomas  and  colleagues  and  students  via  transmission 
electron  microscopy  (TEM)  methods.  Relating  to  the  broad  range  of  materials  covered  today  by  The 
Minerals,  Metals  and  Materials  Society  (TMS),  various  characterization-type  results  are  shown  to  be 
obtainable  for  millimeter-to-centimeter  size  crystals  exhibiting  cither  metallic,  ionic,  covalent,  or 
molecular  (energetic)  type  bonding  between  their  respective  constituent  atoms  or  molecules. 

Introduction 

In  retrospect,  an  especially  important  technical  conference  was  held  on  the  topic:  “Direct  Observation  of 
Imperfections  in  Crystals”  [1],  that  was  sponsored  by  the  Chemistry  and  Physics  of  Metals  Committee 
of,  then.  The  Metallurgical  Society  (TMS),  now.  The  Minerals,  Metals  and  Materials  Society,  American 
Institute  of  Mining,  Metallurgical  and  Petroleum  Engineers,  in  St.  Louis,  MO,  during  March  1-2,  1961. 
At  the  time,  rather  early  dislocation  observations  were  being  reported,  in  connection  with  property 
measurements,  via  the  several  techniques  of  transmission  electron  microscopy  (TEM),  field  ion 
microscopy  (FIM),  and  x-ray  diffraction  topography  (XRDT),  as  were  recently  pioneered  along  with 
dislocation  etch  pitting  investigations.  Professor  Gareth  Thomas  and  colleagues  reported  results  on 
dislocation  structures  observed  in  deformed  and  recovered  molybdenum  polycrystals,  presaging  many 
other  observations  to  come  on  a  wide  variety  of  engineering  materials  involved  in  multiple  technical 
uses,  as  evidenced  by  the  papers  delivered  in  the  present  meeting  by  colleagues  and  students.  Also, 
Professor  George  Smith,  Oxford  University,  U.K.,  is  reporting  in  the  present  TMS  symposium  on 
advances  made  via  FIM. 

I'he  present  paper  provides,  on  the  one  hand,  updated  information  on  characterization-type  results 
obtained  via  several  XRDT  techniques,  following  a  number  of  post-founding  reports  made  in  reference 
1 1],  for  example,  for  transmission  XRD  T  of  sliced  silicon  wafers  by  A.E.  Jenkinson  and  A.R.  Lang.  The 
technical  conference  proceedings  was  co-edited  by  J.H.  Wernick  and  J.B.  Newkirk,  who,  in  the  latter 
case,  provided  encouragement,  because  of  his  research  effort  with  Berg-Barrett  (B-B)  reflection 
topography,  to  the  present  author  and  colleague,  J.M.  Schultz,  for  inclusion  of  first  results  reported  for 
dislocation  observations  in  cleaved  zinc  crystals.  At  the  conference,  G.  Borrmann  and  K.  Lehmann 
presented  XRDT  results  obtained  on  a  very  nearly  perfect  gennanium  crystal  via  .so-called  anomalous 
transmission  of  x-rays;  and,  U.  Bonse  reported  on  crystal-monochromated  XRDT  results  obtained  on  a 
similar  low  dislocation  content  silicon  crystal. 

_ nano  AKO 

submitted  for  “IVII(:U()STRUCTUIU\L  DESIGN  OF  ADVANCED  MATERIALS;  A 
Commemorative  Volume  on  Professor  C.  Thomas’  Seventieth  Birthday,  Elsevier  Science,  2003. 
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Here,  reflection,  transmission,  and  monochromatic  (beam  conditioned)  topographic  observations  are 
presented  for  characterization  of  dislocation  influences  in  zinc,  sapphire,  cyclotrimethylenetrinitramine 
(RDX),  and  ammonium  perchlorate  (AP)  crystals,  plus  layer  deposition  and  diffusion  implantation 
generated  strains  in  otherwise  perfect  silicon  crystals;  in  the  latter  case,  achieved  via  line  modified 
asymmetric  crystal  topography  (LM-ACT).  In  addition,  mention  is  made  that  the  very  appreciable 
advances  now  achievable  via  XRDT  and  related  x-ray  techniques  were  reviewed  in  1995  in  a  Special 
Issue  of  Physics  Today  [2],  including  reference  to  synchrotron  radiation  source  capabilities  such  as  have 
been  available,  for  example;  CCLRC  at  Daresbury,  U.K.;  SSRL  at  Stanford  University;  CHESS  at 
Cornell  University;  NSLS  at  Drookhaven  National  Laboratory;  LUR£  at  Orsay,  FR;  HASYLAB  at 
DESY,  Hamburg,  DE;  ESRF  at  Grenoble,  FR;  and,  the  Photon  Factory,  Japan. 

Results 

Berg-Barretl  (B-B)  Observations  in  Zinc 

A  stereographic  projection  description  of  the  B-B  reflection  XRDT  technique,  as  applied  to  dislocation 
observations  made  through  the  (0001)  cleavage  surfaces  of  zinc  crystals,  is  shown  in  Figure  1.  In  the 
Figure,  the  inset  drawing  depicts  the  physical  set-up  corresponding  to  the  left-side,  zero  layer  case,  of 
cobalt  Ka  radiation  being  incident  to  the  crystal  (0001)  at  6.1°  and  reflected  from  the  inclined  (T013)  at  a 
Bragg  angle  of  41.6°  so  as  to  give  a  diffracted  beam  at  77.1°  from  the  (0001)  specimen  surface,  thus 
traveling  a  distance  as  small  as  possible  to  the  recording  nuclear  emulsion  film.  Point-by-point 
connection  between  the  local  diffracted  intensity  and  the  real  crystal  surface  region  is  achieved  by 
enlargement  of  the  developed  image  in  an  optical  transmission  microscope  [3]. 


Figure  1:  Stereographic  standard  (0001)  projection,  with  inset  physical  set-up,  for  Berg-Barrett 
reflection  x-ray  topography  of  zinc  employing  characteristic  cobalt  Ka  radiation. 
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liutiviJual  dislocations  and  small  angle  dislocation  suhgrain  boundaries  fonned  during  solidilication 
along  the  dilficult  (0001  j  crystal  growth  direction  for  an  eight  millimeter  diameter  zinc  crystal  are 
shown  in  i-igure  2a,  to  be  compared  with  the  characterized  dislocation  boundary  structure  in  b'igure  2b. 
I  he  individual  ilislocation  lines  within  the  subgrain  volumes  arc  revealed  by  the  (blackened)  localized 
enhancement  ot  the  rellected  intensity  caused  by  reduction  in  (the  perfect  crystal)  primary  extinction. 
I  he  occurrence  ot  dislocation  loops  [4|  in  such  cr)'stals  when  chemically  poli.shcd,  and  occurring  also  in 
separate  (0001)  surface  oxidation  experiments  |5].  are  observed  at,  say,  smallest  ten  micrometer 
iliamcters.  to  be  comparable  to  dislocation  loop  observations  made  via  TF.M.  The  dislocation  subgrain 
boundaries  are  generally  observed  to  have  reduced  zones  of  extinction  contrast  but  are  cither  white  (via 
reduced  intensity)  or  black  according  to  whether  the  adjacent  subgrain  reflections  are  separated  or 
overlapped  in  the  I'tlm  record.  Stereographic  (6|  and  diffracted  vector  |7|  analyses  have  been  given  for 
s|>ecillcation  of  subgrain  inisorientations  in  zinc  and  nickel  crystals,  respectively.  Different  levels  of 
subgrain  structures  were  ob.served  in  the  nickel  crystals,  as  first  reported  via  Debye-Scherrer  x-ray 
results  for  nickel  polycrystals  by  Weissmann  jS]. 


Figure  2:  Individual  dislocation  and  subgrain  boundary  structure  within  a  ( 10T3)  H-B  image  through  the 
(0001 )  cleavage  surface  of  a  zinc  crystal  solidified  along  the  (0001 )  crystal  growth  direction. 


Lang  Transmission  XRDT  of  (Magnesia  ami)  Sapphire 

Hmployment  of  the  Lang  transmission  XRDT  method  [9],  as  utilized  with  a  micro-focus  x-ray  generator, 
is  shown  via  the  stereographic  projection  method  lor  {200}  diffraction  of  molybdenum  Ku  radiation 
directly  incident  at  a  small  angle  to  the  (001)  crystal  surface  nonnal  in  Figure  3,  also  whereby 
synchronized  traversing  ot  the  crystal  specimen  and  recording  film  occurs  on  either  side  of  collimating 
lead  screen  slits.  With  this  technique,  depth  perception  of  the  dislocation  line  structures  is  achieved  via 
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stcrco-pair  topographic  images  indicated  to  be  obtainable  in  Figure  3  by  the  inset  physical  set-up 
actually  applying,  relative  to  the  identified  diffracting  plane  normal  in  the  projeetion,  for  an  opposite 
(negative)  x-ray  source  direction  being  positioned  on  the  right  equatorial  side  of  the  [001]  center  of  the 
standard  projection. 


Figure  3:  Stereographic  standard  (001)  projection,  with  conjugate  inset  physical  set-up  for  Lang 
transmission  x-ray  topography  of  magnesia  employing  characteristic  molybdenum  Ka  radiation. 

Figure  4  shows  individual  dislocation  bundles  spreading  outward  from  a  central  seed  as  accomplished 
for  a  sapphire  crystal  grown  by  a  chemical  vapor  deposition  technique  [10].  The  crystal  growth  process 
was  interrupted  at  certain  stages  to  show  the  nucleation  of  addition^  dislocation  bundles  then  spreading 
outward  also  from  new  points  established  during  the  stopped  operation,  as  modeled  by  Klapper  [11], 
mainly,  for  organic  crystal  growth.  Comparative  B-B  and  Lang  topographs  were  obtained  also  on  the 
same  crystal  slice  [12].  By  reduction  of  the  Lang-type  lead  exit  slits  to  smaller  width  than  a  reasonably 
thick  sapphire  crystal  slice,  anomalous  transmission  imaging  of  a  lesser  dislocation-induced  diffracted 
intensity  was  demonstrated  in  the  same  image  for  dislocation  lines  close  to  the  beam  entrance  side  of  the 
crystal  superposed  on  the  enhanced  intensity  of  dislocation  lines  close  to  the  crystal  exit  side  [13]. 

Hardness  Strains  in  RDX  and  A  P  Crystals 

A  stereographic  projection  for  (fS  00)  B-B  reflection  with  copper  Ka  radiation  of  crystal  regions 
surrounding  diamond  pyramid  hardness  indentations  put  at  various  applied  load  values  into  an  RDX 
(210)  crystal  solution-growth  surface  is  shown  in  Figure  5a,  along  with  the  recorded  B-B  image  in 
Figure  5b  showing  very  limited  .spatial  extent  of  the  cumulative  dislocation  strain  fields  at  the 
indentation  sites.  Such  indentations  produced  cracking  on  {001  }c  and  {241}  cleavage  surfaces,  as 
indicated  by  the  filled  points,  for  plane  normals,  and  particular  subscripted  identifications  in  Figure  5a. 
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Assessinciu  of  the  while-lo-black-to-background-gray  local  diffraclcd  intensities  at  the  indentation  sites, 
along  with  complementary  etch  pitting  results  [14|,  provided  llrst  experimental  evidence  of  very 
restricted  dislocation  slij)  occurring  in  RDX  and  related  energetic  crystals. 
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Figure  4:  An  (0330)  Lang  transmission  topograph  of  dislocation  bundles  spreading  from  nucleation 
sites  induced  by  growth  interruptions  for  a  sapphire  crystal  produced  by  chemical  vapor  decomposition. 


Figure  5a:  Stereographic  projection  for  ( 1 0  00)  R-B  image  of  hardness  indentations  in  an  RDX  crystal. 
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Kijiure  5b:  While-lo-hluck-lo-gray  diffraction  contrast  for  plasticity/cracking  at  hardness  indentations. 

A  significantly  greater  extent  of,  albeit  anisotropic,  dislocation  slip  and  cracking  behavior  is  observed 
for  two  similar  type  diamond  pyramid  indentations  put  into  a  clean  cleavage  facet  of  an  (ionic)  Al' 
oxidizer  crystal,  as  shown  in  the  matched  optical  and  (523)  H-H  images  of  figures  6a,b  —  with 
interesting  interpreted  consequence  for  the  role  of  dislocation  slip  being  responsible  for  the  observed 
cleavage  cracking  [  1 5|!  The  chosen  orientation  of  the  indenter  diagonals  for  the  pair  of  indentations  put 
into  the  (210)  surface  has  produced  (001 )  cleavage  cracks  emanating  only  from  one  side  of  the  diagonal 
edges  because  of  favorable  dislocation  slip  on  inclined  intersecting  systems  occurring  only  on  that  one 
side,  as  clearly  evidenced  also  by  the  larger  size  of  the  .same  half  of  the  residual  indentations,  riuis, 
otherwise  brittle-type  cleavage  cracking  has  occurred  in  the  indentation  strain  Held  where  greater 
dislocation  strain  is  evidenced  in  the  x-ray  topograph  and  is  associated  with  the  stress  concentration  ol' 
dislocation  pile-ups  against  reacted  dislocation  obstacles  at  the  slip  system  intersections. 


LM-AC'T  for  perfect  .silicon  crystal  devices 

A  different  research  effort  from  that  of  employing  rellection  XRDT  for  monitoring  dislocation  structures 
came  from  the  consideration  that  such  topographic  imaging  could  be  applied,  if  the  resolution  were 
improved,  to  monitoring  elastic  strains  accomitanying  layer  depositions  and/or  implantation  .strains  in 
electronic  devices  fabricated  on  e.ssentially  perfect  silicon  crystal  surfaces  [lb).  I'or  the  purpose, 
advantageous  u.se  was  made  both  of;  ( 1 )  the  dynamical  theory  prediction  for  reducing  the  angular  spread 
for  the  total  rellection  of  x-rays,  as  originally  proposed  to  play  an  important  role  in  observing  dislocation 
straiti  I'telds  1 17];  ainl.  (2)  further  reduction  of  adverse  geometrical  vertical  divergence  determined  at  the 
x-ray  source  by  employing  a  horizontal  line  geometry  in  a  so-called  line  modified  -  asymmetric  crystal 
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I'igiiri;  6a:  Optical  micrograph  oflwo  microindentations  on  a  clean  (210)  AP  crystal  surface;  and, 

61):  (?23)  li-U  topograph  of  dislocation  strains  producing  cleavage  cracking  at  slip  system  intersections. 

topography  (l.M-ACf)  design  |18|.  Predictions  from  such  dynamical  theory  parameters  account  for 
reduction  of  ohlitiue  "rocking  curve  wiilths”  in  synchrotron  radiation  applications  |19J  and  are  being 
carried  forward  currently  for  evaluation  of  energetic  crystal  powder  diffraction  experiments  |2()|. 

I'igure  7  shows,  on  top  and  bottom,  the  physical  experiment  and  combined  dispersion  surface 
construction,  respectively,  employed  for  application  to  investigating  diffraction  contrast  within,  and 
above,  the  substrate  of  an  otherwise  perfect  silicon  crystal  parametric  test  chip  [t8|.  fhe 
asymmetrically-cut  silicon  crystal  beam  conditioner,  in  reconstruction  of  the  beam  probe,  serves  the 
purpose  of  reducing  the  spatial  width  of  the  incident  horizontal  beam  while  maintaining  a  small  vertical 
divergence  sighted  at  the  x-ray  source,  fhe  theoretical  angular  rellecting  width  exiting  the  beam 
conditioner  is  just  greater  than  the  atigular  range  of  acceptance  for  total  rellectiou  frotn  the  perfect 
silicoti  crystal  sitbstrate  device,  tiow  with  a  further  dytiamical  theory  ba.sed  reductioti  of  the  perfect 
ci  jstal  ret]ected  beatn  width  directed  onto  the  emulsioti  ftltn  record. 

1  he  described  l.M-.ACT  set-up  has  produced  an  image  resolution  at  1-2  tnicrometers  involvitig  oOOX 
magnilication  of  exposed  etuulsion  images  in  an  optical  transmission  tuicroscope  |2I  |.  An  example  is 
showti  iti  I'igure  8a. b  of  a  fitie  line  strueture  re.solved  near  to  the  edge  of  a  depositctl  silica  graticttle 
tnttrkitig  etilarged  iti  part  in  8b  with  two  It)  micrometer  scale  marker  separatiotis.  fhe  graticule  tnarkers 
are  revealed  because  of  elastic  straitis  produced  within  the  silicon  crystal  substrate,  fhe  iticidcnt  x-ray 
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beam  direction  in  the  Figure  is  from  left-to-righl  along  the  horizontal  lengths  of  the  fine  line  structure. 

I  he  imaged  line  structure  extends  past  the  true  vertical  edge  of  the  line  pattern  because  of  the  beam 
probe  having  traveled  at  a  small  angle  tlirough  a  depth  of  ~2  micrometers  in  the  silicon  crystal  substrate. 
The  graininess  in  the  Figure  corresponds  to  the  developed  grain  size  of  the  nuclear  emulsion.  Other 
demonstrated  application  of  the  high  resolution  LM-ACT  method  has  been  to  monitor  sequential 
implantation  deposition  steps  for  fabrication  of  such  parametric  [22]  and  related  devices  [23]. 


TOPOGRAPHIC 
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ASYMMETRICALLY 
CUT  SIUCON  BEAM 
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Figure  7:  The  l.M-ACT  system  applied,  via  favorable  dispersion  surface  geometries,  to  reducing  the 
angular  width  for  total  reflection  of  x-rays  so  as  to  obtain  improved  imaging  of  a  parametric  test  chip. 
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Figure  8a, b:  High  resolution  LM-AC'F  images  obtained  of  detailed  structure  in  a  parametric  test  chip. 


Suiiiiiiary 

With  relation  to  pioneering  XRDT,  FIM,  and,  especially,  ThM  defect  observations,  such  as  reported  in 
the  latter  case  by  Professor  Gareth  fhomas  and  colleagues  and  students,  a  number  of 
experimental/modcl  XRD  T  results  are  reported  here  for  a  variety  of  relatively  perfect  crystals  of  quite 
different  materials.  I'he  presented  results  cover  observations  made  on:  (1)  individual  dislocation  strain 
fields;  (2)  dislocation  subgrain  boundary  geometries/misorientations;  (3)  cumulative  dislocation 
interactions/cracking  at  residual  hardness  indentations;  and,  (4)  elastic  strains  in  an  electronic  device 
fabricated  on  an  essentially  dislocation-free  silicon  crystal. 
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The  seeds  of  tornado  prevention 

J.  Gregory  Glenn  and  Ronald  W.  Armstrong 

AFRL/MNME,  Eglin  Air  Force  Base  EL  USA 

Abstract;  The  role  of  eleotrifieation  in  produeing  tornadoes  is  assessed  in  order  to 

present  a  new  inventive  "bottom-up"  method  of  fly-ash-like  lighter-than-air 
partiele  eloud  seeding  for  tornado  prevention. 

Keywords:  Cloud  eleotrifieation,  lightning,  ionization,  oondensation,  updraft,  vortex 
development,  fly-ash-like  nano-particles. 


1.  Introduction 

Tornadoes  generally  appear  after  severe  frontal  thunderstorms  involving  copious  amounts 
of  lightning.  Vonnegut  (1960)  points  out  the  insightful  association  made  in  a  referenced 
1 840  comment  by  Peltier  that  “Everything  proves  that  the  tornado  is  nothing  else  than  a 
conductor  formed  of  the  clouds,  which  serves  as  a  passage  for  a  continual  discharge  of 
electricity  ...”.  Figure  1  shows  the  designated  oldest  known  photograph  of  a  tornado 
(NOAA,  2003)  obtained  near  to  one -half  century  after  Peltier's  concluding  comment.  As 
will  be  seen  in  subsequent  description  given  here,  Peltier's  quoted  one  sentence 
elucidation  is  not  quite  correct  but  modern  observations  do  show  that  tornado 
development  does  trail  the  peak  of  lightning  activity  (Goodman  1999;  Buechler  et  ah, 
2000).  And,  despite  the  association  of  thunderstorm-generated-lightning  and  tornado 
occurrences,  the  fundamental  origins/forming  conditions  of  tornadoes  are  not  yet  fully 
deciphered.  For  example,  recent  description  of  electrical  features  associated  with 
tornadoes  (Winn  et  ah,  2000)  presumably  once  developed,  allows  for  conflicting  reports 
of  "remarkably  intense  lightning  ..."  and  the  "absence  of  lightning  in  the  storm". 

Recent  lightning-tornado  association  (Goodman  1999,  Buechler  et  ah,  2000)  is  to  be 
compared  with  previous  description  (Davies-Jones  and  Golden,  1975)  of  rare  electrical 
activity  near  or  within  tornadoes.  Current  emphasis  is  on  the  fundamental  role  of  all 
possible  parameters  on  determining  a  critical  updraft  feature  even  of  the  thunderstorm 
origin  while  leading  possibly  to  any  tornadic  activity  (Dye  et  ah,  1986).  Such  basic 
differences  in  attribution  of  “on  or  off’  electriflcation  in  tomadogenesis  accrues, 
undoubtedly,  in  part,  to  the  need  for  coming  together  in  not  quite  chaotic  circumstance  of 
a  combination  of  tornado  determining  parameters/conditions.  The  situation  bears 
similarity  with  the  anecdote  of  deflning  an  elephant  from  different  touches  on  disparate 
external  parts  —  despite  one  association  (Brooks  1951)  of  the  tornado  funnel  with  an 
elephant’s  trunk. 

Here  we  flrst  synthesize  a  series  of  events,  based  on  previously  reported  tornadic 
studies,  by  which  a  tornado  might  be  envisioned  to  originate  —  no  small  feat  because  of 
the  complexity  of  involved  parameters—  and,  then,  propose  a  novel  method  of  defeating 
such  tornadic  development  at  embryonic  operation. 
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Figure  1.  Oldest  known  photograph  of  a  tornado,  28  August  1884,  at  loeation  22  miles 
southwest  of  Howard  SD,  as  obtained  from  the  NOAA  Photo  Library; 
wea00206.  Historic  NWS  Collection. 


2,  Thunderstorm/tornado  elements 

Figure  2  illustrates  the  normal  thunderstorm  situation  of  a  build-up  of  positive  upper  and 
negative  lower  charge  in  the  thunderstorm  cloud  structure.  The  situation  leads  eventually 
to  lightning  strikes,  in  this  chosen  schematic  case,  being  shown  to  initiate  within  the 
cloud  (Uman  1994)  that,  at  thunderstorm  stage,  is  in  a  low-pressure,  multi-phase, 
(ionized)  plasma/gaseous/particulate  state.  The  vertical  charge  separation  is  thought  to 
involve  mostly  hydrogen/hydroxyl/nitrous  ion-associated  water/ice  molecule  clusters  or 
particles,  proposed  to  generate  the  charges  through  strongly-convective  ice  particle 
collisions  (Dye  et  ah,  1986).  In  such  charge  generation  process  for  establishment  of  the 
electric  field  within  the  thunderstorm  cloud,  a  fundamental  role  is  assigned  to  existence 
of  a  strong  intra-cloud  updraft,  whose  collisional  charge  generation  process  over-rides,  at 
least,  at  first,  the  reverse  flow  of  negative  charge  to  the  lower  cloud  level.  In  any  event, 
the  vertical  charge  separation  is  thought  to  develop  in  more  complicated  convective 
churning  of  the  “vaporous  soup”  than  simply  being  the  result  of  an  eventual 
gravitationally-determined  lower  heavier  hydroxyl  ion  attachments  as  compared  to  higher 
lighter  hydrogen  ion  ones.  Strong  friction  forces  are  proposed  to  be  associated  with 
minute  hail  nuclei  interactions.  The  buoyant  force  of  the  updraft  is  sufficient  to  support 
quite  large  hailstones  or  other  graupel  prior  to  their  release.  Cosmic  radiation  apparently 
contributes  little  to  the  process  (Tammet  1994). 

Figure  3a  is  a  more  detailed  elaboration  of  the  charge  build-up  occurring  within  a 
thunderstorm  “supercell”  structure,  also  showing  the  importance  of  wind  shear  properties 
as  part  of  the  updraft  in  contributing  to  tornado  development  (Klemp  1999).  The  build¬ 
up  of  vertically  separated  charge  densities  within  the  potential  thunderstorm  proceeds 
until  dielectric  breakdown  of  the  cloud  resistance  occurs  and  lightning  strikes  are 
produced,  as  indicated  in  Figure  3b.  Intra-cloud  strikes  are  of  special  interest  in  modern 
researches  that  show  essentially  no  correlation  with  cloud-to-ground  strikes  (Goodman 
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Figure  2.  Partial  sequence  of  timed  steps  in  a  cloud-to-ground  lightning  strike,  giving 
emphasis  here  to  intra-cloud  initiation  of  the  lightning  strike;  time  is  in  milli¬ 
seconds;  see  Davies-Jones  and  Golden,  1975. 


1999;  Buechler  et  ah,  2000).  Flere,  we  focus  on  the  intra-cloud  strikes  for  the  additional 
reason  that  such  strikes  produce  localized  cylinders  of  ionized  particles,  particularly,  of 
hydrogen  and  hydroxyl  ions,  which  then  are  accelerated  as  charge  carriers  within  the 
otherwise  global  electric  field  of  the  thundercloud  structure. 

A  role  for  consideration  of  the  smaller,  “unseen”  population  of  molecular  elements  has 
been  attributed  to  Rathbun  in  a  later  comprehensive  analysis  of  electrical  influences  on 
tornadic  behavior  (Wilkins  1964).  The  hydrogen  ions  are  accelerated  downward  and  the 
hydroxyl  ions,  up,  even  if  attached  to  water  or  other  molecular  clusters.  The  easier  flow 
of  hydrogen  ions  downward  produces  re-combination  to  FI2O  molecules  and,  thence,  to 
ice  or  water,  with  an  approximately  twenty  times  greater  exothermic ity  for  the  ionic  re¬ 
combination  as  compared  with  the  molecular  condensation.  The  liberated  heat  at  the 
lower  cloud  level  should  contribute  significantly  to  added  strengthening  of  the  updraft, 
otherwise  found  reasonably  difficult  to  account  for,  temperature-wise,  only  on  a 
condensation  basis  (Vonnegut  1960;  Wilkins  1964).  Greater  influence  on  the  mixed 
vapor  flow,  molecular-cohision-wise,  is  produced  by  the  heavier  hydroxyl  ions  that  are 
driven  upwards  to  effect  their  re-combination  at  the  upper  cloud  level  in  the  global 
electric  field.  Thus,  an  electrically-driven  addition  to  the  updraft  is  produced  with 
important  temperature  and  pressure  components.  The  encompassing  cloud  and  charge 
carriers  operate  as  an  electrostatic  motor  (Watkins  et  ah,  1978). 

An  essentially  vertical  lightning  strike  produces  a  radial  distribution  of  ions  traveling 
upwards  and  downwards  on  the  cylinder  circumferential  surface  surrounding  the  post¬ 
strike  “centerline”,  very  probably,  having  peak  distribution  within  a  particular  radius 
increment.  Perhaps  there  is  a  vortex-type  screw  character  to  the  lightning  strike,  though 
as  yet  only  referenced  (Wilkins  1964)  to  an  electrodynamic  theory  proposed  by  Rathbun. 
Subsequent  pioneering  laboratory  experiments  to  study  tornado  characteristics  have  been 
performed  by  way  of  demonstrating  the  character  of  electric  discharges  occurring  inside 
of  vortex-stabilized  vertical  arcs  (Wilkins  1964;  Watkins  et  ah,  1978)  that  have  provided 
indication  of  updraft  enhancement  of  lightning  strikes.  We  would  term  these 
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Figure  3. a.  J.  Klemp  three-dimensional  schematic  illustration  of  horizontal  wind  velocity 
influence  on  tornado  development,  including  charge  separation,  as  utilized  by 
H.B.  Bluestein  (1999). 


Figure  3.b.  Adaptation  of  Figure  3. a.  to  show  an  intra-cloud  lightning  flash  producing 
charge  carriers  to  activate  the  thunderstorm  “motor”  in  pre-tomadic 
development. 
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experiments,  along  with  most  of  the  refereneed  information  obtained  until  the  present 
time,  as  post-tornado  eharaeterization.  Aireraft  measurements  of  lower  eloud  positive 
eharges  have  been  attributed  to  the  “deposition  of  positive  eharge  by  lightning”  (Mo  et 
ah,  2002). 

Now,  the  deseription  given  above  applies  only  for  a  single  intra-eloud  lightning  strike 
among  an  inereasing  number  assoeiated  with  greater  thunderstorm  strength.  The  single 
vertieal  “motor  eell”  is  among  many  distributed  side-by-side  in  the  super-sized 
thundereloud  —  and  these  vortex-type  eells  repel  eaeh  other  with  a  foree  proportional  to 
whatever  vortex  eomponent  of  strength  may  have  been  aequired.  Again,  beeause  of 
previous  interest  in  lightning-assoeiated  tomadie  aetivity,  researeh  study  of  the 
assoeiation  has  indieated  no  reason  for  repeated  strikes  to  oeeur  within  a  onee  struek 
vortex  (Wilkins  and  MeConnell,  1968).  However,  the  meehanism  that  is  being  proposed 
here  for  an  individual  eell  to  aehieve  tornado  status  is  one  of  repeated  lightning  strikes 
having  to  oeeur  within  the  same  eellular  zone  of  an  earlier  strike.  Sueh  improbability 
relates  to  the  observation  of  tornado  development  being  partieularly  assoeiated  with 
extensively  repeated  lightning  aetivity  within  a  given  region  (Goodman  1999;  Bueehler  et 
ah,  2000).  And,  along  the  way  towards  tornado  development  is  the  tell-tale  beginning  of 
intensifieation  of  loealized  updraft  experieneed  at  the  lower  eloud  surfaee  of  the 
thunderstorm,  thermally-assisted  by  the  (downward)  arrival  of  hydrogen  ions  for  re- 
eombination/eondensation. 

Thus,  in  summary,  intra-eloud  lightning  is  a  neeessary  but  not  suffieient  requirement  for 
tornado  initiation.  The  need  for  multiple  strike  aetivity  within  an  individual  pre-tomadie 
eell  relates  to  the  relative  rarity  of  tornado  oeeurrenees.  Tornado  generation  is  attributed 
to:  [1]  suffieient  eharge  separation  in  the  thunderstorm,  with  aeeompanying  internal 
atmosphere,  to  generate  eopious  lightning  strikes;  [2]  repeated  intra-eloud  strikes  within  a 
developing  tornadie  eylinder  for  thermally-assisted  build-up  of  vortex  strength;  and,  [3] 
re-eombination  of  deposited  positive  eharge  at  the  lower  eloud  surfaee  for  temperature 
rise  plus  heavier  negative  ion  eollisional  flow  upwards  to  eontribute  vertieal  pressure  — 
all  aided  by  lateral  wind  eonsiderations  (Klemp  1999). 


3,  Charge  recombination/condensation 

A  tornado,  however  formed,  eontains  mueh  less  energy  than  the  thunderstorm  that 
produees  it.  Modern  evidenee  for  tornado  oeeurrenee  following  signffieant  lightning 
aetivity  (Goodman  1999;  Bueehler  et  ah,  2000)  is  shown  in  Figure  4.  The  strong 
eorrelation  has  led  to  suggested  orbital  lightning  monitoring  as  a  possible  tornado 
warning  system.  However,  while  elimination  of  the  global-type  eharge  separation  that 
oeeurs  in  thunderstorm  development  is  "out-of-the-question",  there  is  opportunity  of 
interfering  with  tornado  initiation  eaused  by  lightning  "eharging"  of  potential  vortiees. 
The  basie  idea  is  to  effeet  both  lightning  suppression  and/or  reeombination  of  lightning- 
initiated  eharges  by  supplied  nanometrie  fly-ash-like  partieles  (Akesson  et  ah,  2002)  put 
into  the  embryonie  updraft  zone,  loeated  just  below  the  negatively  eharged  lower  eloud 
level.  As  a  start,  the  insulating  eharaeter  of  the  eeramie  eonstituents  of  fly-ash  partieles 
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intracloud  and  cloud— to— ground  lightning  detected  by  OTD 
as  it  passed  over  the  tornadic  storm  cell. 

Cloud— to— ground  lightning  associated  with  the  tornadic  storm  cell 
that  was  detected  by  the  National  Lightning  Detection  Network. 


Figure  4.  Optical  transient  detector  (OTD)  plot  reported  by  NASA  Space  Science  News, 
as  obtained  from  orbital  observations  made  over  a  storm  in  Oklahoma.  "OTD 
registered  a  sharp  increase  in  lightning  flashes  as  the  storm  system  built,  and  as 
clouds  rose,  then  a  sudden  drop  as  the  system  collapsed  and  spawned  a 
tornado."  (Goodman  1999;  Buechler  et  ah,  2000). 


should  contribute  to  increase  of  the  dielectric  breakdown  strength  of  the  intra-cloud 
atmosphere.  Also,  the  low  density  and  small  size  of  the  suggested  particles  should 
promote  their  having  an  important  significant  lifetime  in  the  intra-cloud  chamber,  as 
presumed  useful  to  effecting  ion  recombination  and  subsequent  phase  transformation. 

The  proposed  method  of  tornado  prevention  by  a  cloud  seeding  procedure  is  novel  for 
several  reasons;  [1]  essentially  lighter-than-air  nanometrically-sized  particles  are 
involved;  and  [2]  the  particles  are  to  be  injected  into  the  early  updraft  development  of  the 
lower  cloud  level.  Of  course,  it  has  been  known  for  the  last  ten  years  that  cloud  seeding 
has  been  accomplished  quite  effectively  for  hail  suppression  (Sheremata  2001).  Hail  is 
known  to  be  a  by-product  of  the  storm  systems  that  produce  tornadoes  and  such  seeding 
has  been  utilized  to  effectively  reduce  ground  level  damage  through  reduction  of  hail 
stone  sizes.  Such  alteration  of  hail  characteristics  that  has  been  achieved  on  a  large  scale, 
effectively  altering  cloud  energetics  even  at  a  scale  near  to  tornado  spawning 
thunderstorms,  provides  encouragement  for  the  current  proposal.  In  contrast  to  other 
rain-making  cloud  seeding  activities,  both  hail  suppression,  as  referenced  here,  and  the 
proposed  charge  recombination/condensation  mechanism  involves  super-saturated 
atmospheric  conditions  that  need  only  relatively  much  smaller  particles,  even  suggested 
nanometrically-dimensioned  ones,  to  effect  charge  re-combination  and  thence  molecular 
condensation.  The  smaller  particle  sizes  relate  to  the  important  consideration  of  being 
able  to  disperse  a  sufficient  number  of  particles  to  influence  the  properties  of  a  sizeable 
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cloud  structure.  Under  such  severe  atmospherie  conditions,  the  finer  aspects  of  effective 
advantage  through  top-down  eloud  seeding  also  to  involve  epitaxial  aspeets  of  utilizing 
(heavier)  silver  iodide  particles  seem  relatively  unimportant,  although  sueh  eonsideration 
is  not  discounted  among  potential  low-density  candidate  particle  materials. 

The  basic  idea  is  to  foster  reeombination  of  separated  ionized  hydrogen  and  hydroxyl 
constituents  from  water,  along  with  gaseous  water  molecules,  into  water/iee  partieles  in 
the  heavily  supersaturated  eloud  formations  presaging  possible  full  tornadic  activity,  thus, 
producing  significant  energy  dissipation  in  the  formation  of  condensed  products  and 
consequent,  undoubtedly  heavy  rain  or  hail,  instead  of  lightning  strikes  and  destructive 
wind  power.  The  main  consideration  for  accomplishing  the  task  involves  a  form  of 
"eloud  seeding"  —  done  in  a  new  way  with  new  partiele  considerations,  as  follow: 

(i)  To  make  important  contact,  first,  with  the  partieular  constituent  of  a  heavy 
coneentration  of  hydroxyl  ions  at  the  lowest  eloud  surfaees,  the  cloud  seeding  is  to  be 
initiated  by  liberation  of  appropriate  nanometric,  or  near-nanometric,  sized  partieles 
below  the  lower  eloud  levels. 

(ii)  Something-like  nanometric-sized  fly-ash  partieles  (that  are  hollow  ceramic 
cenospheres,  that  unfortunately  are  normally  of  larger  30-300  micrometer  size)  would 
seem  to  be  ideal  from  a  desired  "  lighter- than-air"  property  in  that  an  important  coneept  in 
this  proposed  eonsideration  is  that,  if  tornadic  development  is  imminent,  the  particles 
would  be  dragged  or  pulled  upwards,  in  the  eonvectively-stirred  up-draft  presaging 
tornado  oecurrenee,  thus  going  through  the  clouds  from  the  bottom-up  while  serving  as 
particle  nucleation/growth  sites,  first,  for  hydroxyl  ion  absorption  and  then,  further  up  in 
altitude,  hydrogen  ion  reaction/nucleation  to  form  hail  or  rain. 

(iii)  Perhaps,  a  weak  positive  charge  state  could  be  added  to  such  particles,  first,  to  help 
prevent  their  own  agglomeration  and  then,  when  appropriately  dispersed,  to  help  in 
promotion  of  first  stage  hydroxyl  ion  adsorption.  As  an  aside,  if  only  positively  charged 
particles,  near  to  being  lighter  than  air,  were  liberated  simply  to  combine  with  the 
hydroxyl  ions,  the  advantage  of  gaining  dissipation  of  the  very  appreeiable  ionization 
potential  energy  in  forming  water  or  hail  particles  would  be  lost. 

(iv)  The  proposal  of  employing  nanometric  particles  in  order  to  aecomplish  tornado 
defeat  is  important  for  at  least  two  size-dependent  reasons:  (a)  sueh  small  size  is  in  the 
presumed  smaller,  or  smallest,  partiele  regime  theoretically  capable  of  nucleating  water 
droplets  under  large  supersaturation  conditions;  and,  (b)  there  is  the  more  practical 
consideration  of  being  able  to  carry/disperse  "billions  and  billions"  of  partieles  during 
eonveying  flight. 
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Abstract.  The  involvement  of  the  combined  constituents  RDX  (cyclotrimethylenetrinitramine)  and 
aluminum  particles  in  energetic  material  formulations  provides  a  reason  for  comparing  the  theoretical  x-ray 
diffraction  properties  of  the  two  very  different  materials,  as  determined,  for  example,  in  powder  diffraction 
or  topographical  imaging  experiments.  Also,  of  interest  is  x-ray  determination  of  the  internal  dislocation 
strain  fields  in  the  two  materials  because  of  their  very  different  deformation  and  fracturing  behaviors. 
Thus,  a  comparison  is  made:  (1)  of  the  limiting  theoretical  perfect  and  imperfect  crystal  integrated  x-ray 
intensities  computed  for  symmetrical  diffraction,  on  a  primary  and  secondary  extinction  basis,  respectively; 
and,  (2)  of  the  modeled  dislocation  strain  field  “widths”  affecting  the  x-ray  intensities  for  the  two  materials. 
Beyond  the  complications  presented  both  by  the  lower  diffracting  power  of  the  complex  orthorhombic 
RDX  crystal  lattice  and  its  larger  cell  dimensions,  the  advantage  of  lower  energetic  crystal  dislocation 
densities  is  offset  by  the  greater  extent,  and  influence,  of  the  dislocation  strain  fields. 


INTRODUCTION 

Compared  on  a  mechanical  properties  basis, 
RDX  is  elastically  compliant,  plastically  hard, 
and  brittle  (1),  whereas  aluminum  is  elastically 
stiffer,  plastically  softer,  and  considerably  more 
ductile. 


Figure  1.  An  RDX  scanning  electron  micrograph 


The  difference  in  properties  of  RDX  and 
aluminum  can  be  traced  to  the  bonding  of 
molecules  and  atoms  in  the  respective  lattice 
structures  and,  particularly,  traced  to  the 
properties  of  dislocations  in  the  two  materials. 
Also,  the  deformation  properties  of  RDX  are 
relatively  complicated  (1)  whereas  the 
mechanical  properties  of  aluminum  are  relatively 
simple,  especially,  in  terms  of  relating  those 
deformation  properties  of  aluminum  to  its 
consequent  dislocation  behavior. 


X-RAY  DIFFRACTION 

The  role  of  x-ray  diffraction  in  characterizing  the 
internal  perfection  of  metallic,  ionic,  covalent, 
and  molecular  lattice  structures  is  well- 
established,  even  spanning  early  powder 
diffraction  peak  profile  analyses  and  more  recent 
lattice  imaging  of  dislocation  strain  fields  akin  to 
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the  technique  of  transmission  electron 
microscopy  (2).  For  the  x-ray  case,  the 
theoretical  (relative)  intensity  associated  with 
reflection  of  x-rays  from  an  ideally  perfect 
crystal  is  expressed  as  (3) 


with  a  very  much  reduced  penetration  depth,  4, 
known  as  the  x-ray  extinction  depth,  correspond¬ 
ing  to  the  primary  extinction  of  x-rays 
undergoing  interference  through  multiple 
reflections. 


Pp  =  (8/37i)(l^/sin[20B])('Ftki'A^)(e^/mc^) 

{ 'sin[0B-l]  '/sin[0B+x] }  {P} , 

in  which  X  is  the  x-ray  wave  length,  0b  is  the 
Bragg  angle,  Fhu  is  the  cell  structure  factor,  V  is 
the  cell  volume,  (e^/mc^)  =  2.82  x  lO  '^  cm,  %  is 
the  (negative)  angle  between  the  Bragg  plane  and 
crystal  surface  normal  directions,  and  (^{P}  = 
{(1+cos[20b])/2}  is  the  (average)  polarization 
factor.  Pp  is  small  because  of  the  influence  of  the 
phase  interference  between  multiply-reflected 
rays  in  the  perfect  crystal  lattice. 

The  opposite  limiting  theoretical  diffracted 
intensity,  defined  x-ray-wise  for  an  ideally 
imperfect  crystal,  ps,  applies  for  the  model  of  the 
rays  satisfying  the  Bragg  condition  for  a  single¬ 
time  reflection  only,  as  given  by  (3) 

Ps  =  (l/2p)(XVsin[20B])(Fhu/V)2(eVc")' 

(1  -cot[0B]tan[x]).(;{P}^ 

where  p  is  the  normal  x-ray  absorption 
coefficient  and  all  other  symbols  have  their  same 
value  as  for  pp.  The  greater  penetration  of  x- 
rays  associated  with  Ps  is  compared,  generally 


The  computed  values  of  the  two  limiting  x-ray 
intensities  are  shown  in  Figure  2  for  both  RDX 
and  aluminum  symmetrical  plane  reflections 
obtained  with  copper  Ka  radiation  (4).  Richard 
Gilardi  (5)  supplied  the  values  of  Fhki  for  RDX. 
The  (111)  lowest  angle  aluminum  reflection 
occurs  at  20b  ~  38.5  degrees.  Comparable 
values  of  {hkl}  occur  for  RDX  at  much  smaller 
values  of  20b.  For  the  individual  materials,  of 
course,  ps  >  Pp,  although  more  so  for  aluminum, 
and  its  intensities  are  greater  than  those  of  RDX, 
for  either  case,  principally  because  of  the  larger 
values  of  F^ki.  For  RDX,  the  ratio  ps/pp  shows 
the  particularly  complex  dependence  on  20b  that 
is  indicated  because  of  the  complicated  form  of 
Fhki.  Alternatively,  the  substantial  range  in 
intensity  between  pp  and  ps  for  aluminum, 
particularly  for  the  (111)  reflection,  provides 
ample  room  for  characterization  of  its  (mosaic) 
dislocation  structure. 

The  indicated  situation  for  RDX  of  smaller 
Bragg  angles,  greater  x-ray  penetration  depths, 
lower  diffracted  intensities,  and  multiple 
reflections,  all  combine  to  make  more 
complicated  its  x-ray  diffraction  characteristics. 


•Al  perfect  crystal  diffracted  integrated  irlensity  A  Al  imperfect  crystal  diffracted  integrated  irlensity 

•  RDX  perfect  crystal  diffracted  integrated  intensity  •  RDX  imperfed  crystal  diffracted  integrated  intensity 
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Figure  2.  Integrated  x-ray  intensities,  pp  and  ps,  for  RDX  and  aluminum,  versus  their  diffracting  angles  20b 
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RDX  (unstrained)  AI3 


Figure  3.  Combined  diffraction  traces  for  RDX  and  aluminum  taken  under  otherwise  identical  conditions 

reasonably  compared  with  those  in  Figure  2;  see, 
for  example,  the  comparison  of  the  RDX  (102) 
and  the  aluminum  (111). 

DISLOCATION  INFLUENCES 


The  influence  of  crystal  straining,  further  to 
whatever  level  of  mosaic  imperfection  already 
exists  in  the  crystals  being  examined  by  x-ray 
diffraction,  is  to  lower  the  diffracted  peak 
intensities  for  particular  reflections  and  to 
broaden  the  peak  widths.  Thus,  Figure  4  shows 
the  lowered  intensities,  at  closed  circle  positions, 
put  onto  a  number  of  RDX  peak  intensities 
recorded  before  the  powder  specimen  was 
subjected  to  a  gentle  “crushing”  deformation. 
The  result  is  similar  to  that  obtained  in  a  prior 
experiment  (6)  where  the  focus  was  on  the 
lowering  of  the  recorded  intensity,  and  spread,  of 
the  otherwise  prominent  (400)  diffraction  peak  — 
because  of  identification  of  the  (040)  slip  plane 
and  b  =  a[100]  Burgers  vector  dislocation  in  the 
orthorhombic  RDX  crystal  lattice. 

In  x-ray  imaging,  or  topographical,  experiments, 
whether  at  low  Bragg  angles  in  transmission  or 
large  angles,  0b  ~  45  degrees,  in  reflection,  the 
emphasis  is  on  spatially  resolving  dislocations 
through  their  strain  field  influences  on  locally 


An  example  of  combined  RDX  and  aluminum 
powder  diffraction  traces  obtained  in  separate 
experiments  accomplished  under  otherwise  very 
similar  conditions  is  shown  in  Figure  3.  The 
ordinate  intensity  counts  should  be  considered  an 
arbitrary  scale  but,  for  relative,  order-of- 
magnitude,  comparison  basis,  the  results  may  be 


Figure  4.  Not-crushed/crushed  RDX  intensities. 
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■  Al,  g.b=0 
AAI,  g.bM 
XAI,  g.b=2 
XAI,  g.b=3 
♦  RDX,  b=a[100] 


Figure  5.  Theoretical  dislocation  image  widths  in  x-ray  diffraction  topographs  of  RDX  and  aluminum 


enhancing  the  diffracted  x-ray  intensity  in  the 
surrounding  vicinity  of  the  dislocation  line 
positions.  In  transmission  x-ray  topography 
experiments,  such  reflections  as  the  (200),  (102), 
and  (220)  ones  shown  in  Figures  3  and  4  have 
been  favorably  employed  (7)  whereas  in  surface 
reflection,  at  larger  Bragg  angles,  the  (632)  and 
(10  00)  have  been  necessarily  employed  (8). 

The  theoretical  dislocation  widths  seen  in  x-ray 
topographical  images  have  been  associated  (9) 
with  the  primary  extinction  distance,  4,  specified 
on  a  dynamical  theory  basis  as  for  pp: 

4=(mcV)(7iV/^l(;{P}"Fhki') 

['sin(0B  -  x)’][sin(0  +X)], 

where  all  parameters  have  been  previously 
identified.  The  dislocation  image  width  has  been 
given  by 

w~(^27i)  ghki*b 

with  ghki  being  the  reciprocal  lattice  vector  for 
the  (hkl)  plane. 

Figure  5  shows  the  computed  images  obtained 
for  <110>  type  dislocations  in  aluminum  as 
compared  with  the  a[100]  dislocation  in  RDX. 
In  general,  at  the  relevant  small  angles  shown  for 


aluminum,  the  image  widths  are  much  smaller 
than  those  shown  for  the  various  prominent 
reflections  of  RDX.  The  results  are  interpreted 
to  reflect  the  condition  of  the  dislocation  strain 
field  extending  to  relatively  larger  distances  in 
RDX,  and  to  relatively  greater  purpose  in 
connection  with  the  limited  deformation  and 
cracking-prone  behavior  of  RDX. 
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ABSTRACT: 

Nanometric  fuel-oxidizer  material  systems  that  release  energy  comparable  to  that  of 
HCNO  explosives  are  of  high  current  interest  as  replacements  or  enhancements  for  high 
explosive  materials.  Reaction  and  energy  release  rates  in  such  materials  may  be 
accelerated  by  reducing  particle  size  to  nanometric  scale.  This  increases  surface  area 
for  a  given  mass,  and  reduces  average  diffusion  path  lengths.  One  approach  is  to 
directly  coat  nanometric  fuel  particles  with  a  nanometric-thick  layer  of  an  organic 
oxidizer  or  explosive  material.  In  this  research,  coating  and  crystallization  techniques 
are  being  developed  for  direct  encapsulation  of  metal  fuel  nano-particles.  There  is  also 
high  interest  in  development  of  better  understanding  of  the  physics  controlling  the 
formation  and  growth  of  ultrafme  (including  nanometric)  crystals  of  explosive 
materials,  and  of  composite  nano-structures  composed  of  organic  explosives  and  metal 
clusters,  so  that  techniques  for  the  direct  crystallization  of  materials  as  coatings  on  metal 
fuel  particles  can  be  optimized. 

At  the  Air  Force  Research  Lab,  Energetic  Materials  Branch,  two  methods  are  currently 
being  used  to  recrystallize  RDX,  CL-20,  and  other  organic  materials  onto  spherical 
Aluminum  particles  ranging  in  size  from  30  nm  to  1  |J,m  diameter.  In  each  method,  the 
organic  material  is  first  dissolved  in  a  suitable  solvent.  The  Aluminum  particles  are 
then  added,  with  continuous  agitation.  In  the  first  method,  a  suitable  non-solvent  is  then 
added.  After  precipitation  of  the  organic  material,  the  solvent  and  non-solvent  are  then 
removed  by  vacuum  distillation.  In  the  second  method,  non-solvent  is  not  used;  instead, 
the  solvent  is  slowly  removed  by  distillation  with  constant  agitation.  Materials 
produced  in  these  ways  are  generally  about  three  times  the  diameter  of  the  starting 
nanometric  Aluminum  powder  used.  Powder  x-ray  diffraction  data  indicates  that  the 
materials  recovered  are  of  the  desired  composite  configuration  rather  than  a  mixture  of 
distinct  Aluminum  and  organic  crystal  particles. 

At  TNO  Prins  Maurits  Laboratory,  similar  results  are  being  pursued,  with  some 
differences  in  the  techniques  being  developed.  The  crystallization  of  the  energetic 
material  RDX  on  nanometric  aluminium  particles  is  studied  by  using  a  drowning-out 
process.  It  is  based  on  the  addition  of  an  RDX  solution  in  an  organic  solvent  to  a  non¬ 
solvent,  in  which  the  aluminium  particles  are  suspended.  The  results  obtained  with  this 
technique  under  different  process  conditions  and  the  characterization  of  the  obtained 
products  (by  means  of  e.g.  scanning  electron  microscopy  and  thermal  analysis)  will  be 
presented  in  this  paper. 
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The  Power  and  Strength  of  Energetic/Reactive  Nano-Materials 
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ABSTRACT 


A  number  of  relatively  new  energetic/reactive  material  formulations  employing  nanometric 
particle  size  ingredients  might  be  thought  to  be  somewhat  novel  because  of  the  consideration  of 
bringing  their  energy  release  rate  and  decomposition  times  closer  together.  In  addition,  there  is 
the  advantage  of  such  ultrafme  particle  solids  being  stronger  generally  and,  consequently,  less 
sensitive  to  mechanical  forces  and  deformations,  for  example,  as  experienced  in  drop-weight 
impact  and  shock  wave  tests. 

Thus,  nanometric-sized  energetic  materials,  as  well  as  counterpart  combustible  metals,  and 
reactive  materials,  are  all  of  current  research  interest,  for  one  reason,  relating  to  surface-control  of 
their  decomposition  properties  providing  for  the  possibility  of  orders  of  magnitude  increases  in 
energy  release  rates  at  ultrafme  particle  sizes.  An  indication  of  this  benefit  is  shown  in  Figure  1 
by  extrapolation  to  nanometric  particle  sizes  of  recent  bum  rate  results  obtained  for  aluminum 
particle  combustion  properties  in  both  laboratory  and  gun  propellant  experiments  [1]. 

From  a  research  point-of  view,  the  surface  energy  property  in  the  nano-particle  regime  is  of 
paramount  importance,  as  evidenced  practically  by  nano-particle  processing  difficulties  with 
agglomeration.  Particle  coating  can  alleviate  the  problem,  especially  relating  to  coating  particles 
with  materials  that  have  relatively  lower  surface  energies.  In  this  regard,  a  comparison  of  solid 
surface  energies  (or  stresses)  of  materials  with  their  volume  free  energies  provides  interesting 
insight  into  the  differences  between  the  fundamental  material  lattice  stmctures  and  their 
physicaFmechanical  property  consequences;  see  reference  [2]  for  a  reported  listing. 

And,  then,  dislocation  defects  are  presumed  to  control  the  mechanical  strength  and  fracturing 
behaviors,  even,  of  nanometric-scale  material  particles  and  stmctures  —  thus  presenting 
challenging  areas  of  research  so  far  as  concerns  a  dearth  of  material  testing  methods  and  also 
prediction  of  size-dependent  material  deformation/decomposition  properties.  Recent  atomic 
force  microscopy  (AFM)  observations  on  nano-indented  and  otherwise  deformed  RDX  crystals 
have  indicated  the  presence  of  shock-generated  nanoscale  shear  cracks  with  displacement  vectors 
smaller  than  the  length  of  individual  molecules  [3]. 

1.  R.W.  Armstrong,  B.  Baschung,  D.W.  Booth,  and  M.  Samirant,  “Enhanced  Propellant 
Combustion  with  Nanoparticles”,  Nano  Eetters,  3,  253  (2003) 

2.  W.E.  Elban,  R.W.  Armstrong,  and  T.P.  Russell,  “Plasticity/Interfacial  Energy  Influences 
on  Combustion-Driven  Cracking  of  RDX  Energetic  Crystals”,  Philosophical  Magazine 
A,  78,  907  (1998). 

3.  J.  Sharma,  C.S.  Coffey,  R.W.  Armstrong,  W.E.  Elban,  and  S.M.  Floover,  “Sub- 
Molecular  Fracture  Steps  in  Shock-Shattered  RDX  Crystals  and  Follow-on  Nano- 
Indentation  Evaluation  of  Early  Stage  Plasticity”,  in  Shock  Compression  of  Condensed 
Matter  -  2001,  M.D.  Furnish,  N.N.  Thadhani  and  Y.  Florie,  eds.  (AIP  Press,  N.Y.,  2002) 
Conference  Proceedings  620,  Part  2,  p.  837. 
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Figure  1.  Burn  rate  as  a  function  of  aluminum  particle  size,  also  showing  indicated  pressure 
ranges  for  the  vertical  lines,  and  providing  indication  at  the  smallest  particle  sizes  of  an  inverse 
square  dependence  of  bum  rate  on  particle  diameter  [1],  The  Technanogy  results  are  for  strands 
3.175  mm  diameter  and  6.35-8.47  mm  length  employing  a  standard  ammonium  perchlorate 
formulation  in  88%  solid  condition  after  pressing  to  80-85%  density.  The  French-German 
Research  Institute  of  Saint  Louis  results  apply  for  a  new  double-base  gun  propellant  containing 
15%  either  of  Alex  or  micrometer-sized  aluminum  particles  in  the  same  Pechiney  YX  76 
propellant  formulation  of  55%  nitrocellulose  and  30%  dioxyethylnitramine  dinitrate  (DINA),  as 
measured  at  different  pressures  in  a  high-pressure  closed  vessel  of  50  cm^  volume. 
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Dislocation  properties: 

Pile-up/crack  analogy;  grain  volume/boundary-controlled  dynamics;  and, 

energetic  crystal  interpretations  ^ 
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ABSTRACT 


The  analogous  properties  of  dislocation  pile-ups  and  sharp  cracks  [  1  ]  are  employed 
recently  to  explain  the  increased  indentation  fracture  mechanics  stress  intensities 
measured  at  small  crack  sizes  for  alumina  polycrystals  [2].  For  Armco  iron  employed  as 
a  reference  structural  material,  the  twinning-stress-determined  Hugoniot  elastic  limit 
(HEL)  is  shown  to  follow  an  Hall-Petch  inverse  square  root  of  grain  diameter 
dependence,  in  agreement  with  conventionally  determined  twinning  stress  measurements 
obtained  at  low  temperatures  and/or  high  strain  rates  [3].  Special  dislocation  properties 
for  RDX  and  related  energetic  materials  offer  a  rationale  for  the  crystals  being  plastically 
hard  and  brittle,  as  compared  with  their  otherwise  “polymer-like”  elastic  compliances  [4]. 

References 

1 .  R.W.  Armstrong,  “Plasticity:  Grain  Size  Effects”,  in  Encyclopedia  of  Materials: 

Science  and  Technology  .  (Elsevier  Science  Publishers,  Oxford,  2001),  pp  7103- 
7114.  „ 

2.  R.W.  Armstrong,  “Grain  size  dependent  alumina  fracture  mechanics  stress 
intensity”.  Intern.  Journal  of  Refractory  Metals  and  Hard  Materials,  vol.  19 
(2001),  pp.  251-255. 

3.  R.W.  Armstrong,  W.  Arnold,  and  F.J.  Zerilli,  “Grain  size  dependence  of  shock- 
induced  twinning  stresses  in  Armco  iron”.  Intern-.  Workshdp:  “New  Models  and 
Hydrocodes  for  Shock  Wave  Processes  in  Condensed  Matter”.  Chair:  V. 

Klimenko,  C.S.  Coffey,-A.  Rajendran,  and  D.  Swift,  Edinburgh,  Scotland,  19-24 
May  2002. 

4.  R.W.  Armstrong  and  W.L.  Elban,  “Dislocations  in  energetic  crystals”,  submitted 
for  Dislocations  in  Solids.  Editor:  F.R.N.  Nabarro,  Vol.  12  (Elsevier  Science 
Publishers,  Oxford,  2002). 


176 


W.C.  Unwin,  "On  the  Yield  Point  of  Iron  and  Steel  and  the  Effect  of 
Repeated  Straining  and  Annealing",  Proceedings  of  the  Royal 
Society  (London)  58,  178  (1895). 


The  LOEiD  KELYIIT,  D.O.L.,  LL.D.,  President,  in  the  Chair, 

A  List  of  the  Presents  received  was  laid  on  the  table,  and  thanka 
ordered  for  them. 

'The  following  Papers  were  read  : — 

I.  ‘-On  Measurements  of  small- -Strains  in  the  Testing  of 
Materials  and  Structures.”  By  J.  A.  EwiNG,  M.A.,  F.R.S.,. 
Professor  of  Mechanism  and  Applied  Mechanics  in  the 
University  of  Cambvirlo’p.  P.pnoLrpri  Am-u  9/i  mot; 


177 


Stress  fMPa] 


Length  fmm] 


178 


DISLOCATION  PILE-UP  DESCRIPTIONS  OF  SHEAR  BANDING 


I.  Static  pile-up  results 
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Q.  Gao  and  H.W.  Liu,  "Characterization  of  the  Tip 
Field  of  a  Discrete  Dislocation  Pile-up  for  the • 
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The  fracture  mechanics  stress  intensity  at  small 
grain  sizes  follows  prediction  of  an  H-P 
dependence:  ' 

1 .  Lower  single  crystal  stress  intensity  measurements 
occur  because  of  orientation  influence;  and, 

'  ’  t 

2.  Thermal  strain  influences  apparently  reverse  the  grain 
size  dependence  at  large  grain  sizes. 
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Hardness  stress,  kgf/mm 


RDXi  h,  nrn 


« 

Strain,  (Indentation  dia./ball  dia.) 


=  (V3rt)[(l  +  _  u^)/E,]-'((i/D),  ,  (1) 

wheie  E^,  and  t'j,  £5  are  Poisson’s  ratio  and  Young’s  modulus  for  the  ball  and 
specimen,  respectively,,  r/  is  the  elastic  contact  diameter  and  D  is  the  ball  d.iameter 
1  he  ratio  {d/D)  is  the  effective  hardness  strain.  The  theoretical  hardness  cracking 
stress  is  given  on  an  indentation  fracture  mechanics  basis  as: 

>  [4E,'y/nI){\ -  +  4)]'/^(d/Z))"'/^  (2) 
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